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accepted papers are retained by the Society for one year after publication; unless their 
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latitude being permissible, however, in the case of an author referring to his own work). 
The following examples indicate the style of reference appropriate for a paper and a book, 
respectively :— 

A. Corlin, Zs. f. Astrophys., 1§, 239, 1938. 
H. Jeffreys, Theory of Probability, 2nd edn., section 5.45, p. 258, Oxford, 1948. 
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recommendations of Commission 3 of the International Astronomical Union (Trans. 
1.A.U., Vol. VI, p. 345, 1938). Council has decided to adopt the I.A.U. 3-letter abbrevi- 
ations for constellations where contraction is desirable (Vol. IV, P.3 221, 1932). In general 
matters, authors should follow the recommendations in Symbols, Signs and Abbreviations 
(Loados : Royal Society, 1951) except ieass these conflict with I.A.U. practice. 
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MEETING OF 1954 MAY 14 
Dr J. Jackson, President, in the Chair 


The election by the Council of the following Fellows was duly confirmed :— 


Helmut Arthur Abt, Yerkes Observatory, Williams Bay, Wisconsin, U.S.A. 
(proposed by S. Chandrasekhar); 

Peter John Allender, B.A., B.Sc., Ph.D., 104 Somerset Road, Erdington, 
Birmingham, 23 (proposed by F. J. Wright); 

Raymond Cooper, Royal Fort, The University, Bristol (proposed by 
C. R. Burch); 

Fazli Hussain, Pakistan Military Academy, Kakul, N.W.F.P., Pakistan 
(proposed by A. Hunter); 

Gerald Edward Kron, M.S., Ph.D., Lick Observatory, Mount Hamilton, 
California, U.S.A. (proposed by J. Stebbins); 

Edward Ronald Niblett, M.A., Department of Geodesy and Geophysics, 
Downing Place, Cambridge (proposed by S. K. Runcorn); 

John Beverley Oke, Ph.D., David Dunlap Observatory, Richmond Hill, 
Ontario, Canada (proposed by L. Spitzer); 

Harry Seddon, B.Sc., Royal Observatory, Edinburgh (proposed by W. M. H. 
Greaves); 

Charles Donald Shane, Ph.D., Lick Observatory, Mount Hamilton, Cali- 
fornia (proposed by H. M. Jeffers); 

David Colin Smith, Royal Greenwich Observatory, Abinger Common, 
Dorking, Surrey (proposed by H. Spencer Jones); and 

Henry Tompkins, 6 Crescent Road, Crouch End, London, N.8 (proposed 
by P. A. Moore). 


Ninety-four presents were announced as having been received since the last 
meeting, including :— 
H. Alan Lloyd, Horology and meteorology (presented by the author); and 
V. A. Hughes and E. E. Treadwell, A form of equatorial drive for controlling 
an altazimuth mounted radio-telescope (presented by the Central Radio 
Bureau, Ministry of Supply). 


The President gave a brief address on the award of the Gold Medial, after 
which the Medal was handed to Mr Chester E. Beaman, Chief of Protocol, 
United States Embassy, for transmission to Dr Walter Baade, to whom it had 


34 
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been awarded for his observational work on galactic and extra-galactic objects. 
The President said : 


‘‘ Before proceeding to the reading of papers it is my pleasure to welcome 
here today Mr Chester Beaman from the United States Embassy. He has 
come to receive on behalf of Dr Walter Baade, of the Mount Wilson and Palomar 
Observatories, the Gold Medal of the Society awarded to him. At the Anni- 
versary Meeting in February an Address was given on Dr Baade’s work (see 
p. 370), and so it is unnecessary today for me to say more than that the award 
was for outstanding researches on galactic and extra-galactic objects. But for 
Mr Beaman’s information I would like to say that the Gold Medal is the highest 
award of the Society. It is usually bestowed once a year and in selecting its 
medallist the Council is not restricted by nationality or in any other way. It is 
no new thing for the Medal to be awarded to an astronomer working in California. 
There you have excellent atmospheric conditions for making observations, and 
there you have erected magnificent telescopes. But the Medal is not awarded 
to Dr Baade because of the superb conditions under which he works but for the 
superb use he has made of his opportunities—opportunities which were offered 
to him as a result of his earlier work. Mr Beaman, in forwarding this Medal 
to Dr Baade we would like you to convey to him not only our admiration for 
the work he has already done but our hope that he shall be able to continue and 
extend it with the great 200-inch Hale telescope on Mount Palomar.” 


MEETING OF 1954 OCTOBER 8 
Dr J. Jackson, President, in the Chair 


The election by the Council of the following Fellows was duly confirmed :— 


William Ross Hossack, M.A., Ph.D., David Dunlap Observatory, Richmond 
Hill, Ontario, Canada (proposed by C. A..Chant); and 

Patrick Aloysius O’Brien, M.Sc., Ph.D., Department of Physics, University 
of Natal, Durban, South Africa (proposed by G. Burbidge). 


Two hundred and fifty-one presents were announced as having been 
received since the last meeting, including :— 


H. Alfvén, On the origin of the solar system (presented by the Association of 
Special Libraries and Information Bureaux); 

G. van den Bergh, Eclipses in the second millennium B.C. and how to compute 
them in a few minutes (presented by the author); 

N. O. Bergquist, The Moon puzzle (presented by C. Hamburgers Bogtryk- 
keri); 

F. Chemla-Lamech, Carte topographique de la Lune (presented by the 
executors of the late Mr M. E. J. Gheury de Bray); 

S. G. Hacker, Arithmetical view-points—an introduction to mathematical 
thinking (presented by the author); 

S. G. Hacker, Theoretical astrophysics. Parts I and II (presented by the 
author) ; 

International Union of Geodesy and Geophysics, Photographic atlas of 
auroral forms (bequeathed by the late R. L. T. Clarkson); and 

L. Rudaux, Sur les autres mondes (presented by the executors of the late 
Mr M. E. J. Gheury de Bray). 





PHOTOELECTRIC INVESTIGATIONS OF SOLAR CORPUSCULAR 
RADIATION, II 


M. F. Smyth 


(Communicated by the Astronomer Royal for Scotland) 
(Received 1954 August 13) 


Summary 


Photoelectric measurements of the solar K-line in 1951 and 1953, with 
improved equipment, confirm the conclusion of a previous paper that 
absorption in the violet wing, corresponding to particles approaching with 
speeds between 750-1150 km/sec, does not exceed one-half per cent at any 
time following moderately intense flares. There is also little evidence of 
absorption of this magnitude associated with sequences of recurrent geo- 
magnetic storms. 





1. Introduction.—A previous paper (1)* gave an account of some observations 
in 1949-50 intended to detect or set an upper limit to anomalous absorption in the 
violet wing of the solar K-line due to approaching corpuscles associated with flares. 
These observations failed to reveal in most cases any absorption exceeding 0-5 per 
cent corresponding to particles with velocities in the range 650-1500 km/sec. 

Further series of observations have been made at the Royal Observatory, 
Edinburgh, during 1951 and again during 1953, with an improved apparatus. 
As flare activity declined, 27-day sequences of geomagnetic disturbance became 
prominent, and it was hoped that improved instrumental stability would make it 
possible to investigate also any absorption effect associated with recurrent magnetic 
sequences. The solar source of these sequences is still obscure, despite many 
inconclusive attempts to correlate them with various solar features (2). In the 
absence of evidence to the contrary, however, persistent corpuscular streams from 
unidentified ‘‘ M-regions”’ are usually regarded as responsible for recurrent 
geomagnetic activity. It might then conceivably be possible to detect a recurrent 
absorption effect. 

2. Apparatus.—A six-inch object-glass by Steinheil provided a solar image 
21mm in diameter. For the 1951 observations a coelostat was used, so that the 
spectrograph slit, restricted to 4mm in length, fell centrally and roughly in a 
North-South orientation on the non-rotating 21 mm solar image. In 1953, the 
coelostat being required for other purposes, the equipment was transferred to a site 
using a siderostat, so that the image now rotated slowly during the course of the 
day. This was of no consequence in looking for a corpuscular cloud or stream 
which might be expected to cover the Sun’s disk. 

With a solar image of this size, conduction by the slit jaws was sufficient to 
prevent overheating, and no heat filter was required. 


* Referred to as Paper I in the text, 
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The Cooke spectrograph was of conventional design, with separate collimator 
and camera, and an 80 mm x 50 mm aluminized plane grating, forminga first-order 
spectrum with dispersion 9-7 angstroms/mm. In the focal plane was a set of 
fixed exit-slits, as indicated in Fig. 1 and described later. ‘The slits were all cut 
on a Dumoulin dividing engine in the emulsion of a single fully fogged and 
developed process plate, so that stability of their dimensions and relative positions 
was assured. Light from any one of the slits V, X, Y, R, chosen by means of 
shutters, passed through a field lens to a selected 931-A photomultiplier. A weak 
prism deflected the light from the ‘“‘ Hg” slit, after passing through the field lens, 
to a second 931-A. 


KH Hg 44047 





| | | | SPECTRUM 
VXYR 





Fic. 1.—Evit-slits in focal plane of spectrograph. X and Y lie in the wings of the K-line, 
V and R determine the continuum, while Hg is centred on the mercury line X 4047 as wave-length 
reference and also serves for compensation of seeing in the solar spectrum. 

The photomultipliers were fed from a divider network in the usual manner, 
voltage being obtained initially from dry batteries and later from a stabilized mains 
power unit. The photocurrents went to two calibrated taut-suspension galvano- 
meters, side by side and indicating on the same scale, so that they could be read 
simultaneously. Thus with each reading for one of the V, X, Y, R slits a simul- 
taneous reading of ‘‘Hg”’ (in the continuous spectrum) was obtained. This 
enabled the intensity ratios V: X ; Y; R to be determined with great reduction in 
the error arising from changes in sky transparency, and an estimated accuracy of 
0™-003 (0-3 per cent) for a single ratio was usually attained. In practice, each 
‘‘ observation” involved the average of three or four readings for each ratio. 

3. Wave-length adjustment and reduction of observations.—In the Cambridge 
apparatus (Paper I) the appropriate wave-lengths were brought on to the exit-slits 
by using the solar A-line itself as reference, adjusting to equality a further pair of 
exit-slits placed symmetrically close to the centre of the line. It was found that 
this method of wave-length adjustment produced a small systematicshift depending 
on the “slope” of the spectrum and the consequent distortion of the K-line 
profile. So for the present observations a laboratory wave-length reference was 
used. The slits were arranged so that V, X, Y, R fell in their correct positions 
when the line A 4046-7 of a mercury lamp was centred on the double “ Hg” slit. 
This was really a single slit with a central opaque bar corresponding to the width 
(0-15 mm) of the primary slit. ‘The line was correctly centred when its image fell 
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on the opaque bar and the only light passing was due to ghosts, scattered light, etc. 
Very sharp “tuning” was possible by obtaining the minimum deflection on the 
second galvanometer by means of the wave-length adjustment. 

This adjustment was made for each set of observations, so that the main 
exit-slits were always set on their correct wave-lengths, determined approximately 
by measurement as follows :— 


V: 43870°6 to A 3880-7 
X: A3g918-8 to A3924:0 
Y: A3942°8 to A3948-0 
R: A3gg0°1 to A3993°1 


‘The primary slit-width of about 0-15 mm corresponded to 1-9 angstroms in the 
focal plane. 





Fic. 2.—Positions of exit-slits relative to the H and K-lines, showing effect of variation of “‘ slope’’ 
of spectrum with solar altitude. 


Fig. 2 shows the positions of these slits relative to the 1 and K profiles. V and 
R, placed symmetrically about the K-line, served to determine the continuum; 
X corresponded to absorption byapproaching Ca1i(A)ions with velocities between 
750-1150km/sec, the expected range of particle velocities, while Y, similarly 
placed in the red wing, acted asa control. It will be noted that Y corresponded 
also to Cat1(#) ions with approach velocities between 1500-1900 km/sec, but we 
do not expect such high particle velocities to be associated with minor flares or 
M-regions. The different widths of the slits were chosen to give approximately 
equal photomultiplier responses. 

In estimating possible extra absorption at X or Y it was now necessary to allow 
for the effect of changing atmospheric extinction. Fig. 2 illustrates schematically 
the changes in intensity AV, AX, AY, AR corresponding to a change of solar 
altitude. For each observation, using a logarithmic calibration table for the 
galvanometer readings, the quantities* } log(X/R), $log(Y/R), $log(V/R) 
were formed. In Paper I it was shown that as the solar altitude varies the changes 
in these observed quantities should be proportional, provided that the solar 
spectrum itself is undisturbed. Therefore linear formulae were assumed : 


3 log (X/R)= A + B{Z log(V/R)}, 

$ log (Y/R)=C + D{$ log(V/R)}. 
The values of the constants A, B, C, D were determined by the method of least 
squares, using a whole year’s observations in each case. This procedure was 
preferred to using theoretical values for B, D derived from the Rayleigh A~4 


* All logarithms are to base 10. 
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formula as in Paper I, as it involved no assumption about the “ effective’? wave- 
lengths of the wide slits. ‘‘ Undisturbed”’ values of $ log (X/R) and $ log (Y/R) 
were calculated from the above formulae for each observed value of } log(V/R), 
and subtracted from the observed values of 3 log(X/R) and $log (Y/R). Any 
significant negative X- or Y-residual would then indicate excess absorption 
(in magnitudes) in the corresponding wing of the K-line. 


sec z 
1.0 3.0 5-O 7-0 
T q 














Fic. 3.—Correlation of atmospheric reddening with secz for Sun, 1953 observations. 
V and R are continuous-spectrum intensities at two points about 60 angstroms to the violet and 
red respectively of the K-line, so that lower values of $ log (V/R) indicate increased reddening. 


Open circles : observations before Fune 21. 
Filled circles : observations after Fune 21. 

As a matter of interest, a plot was made (Fig. 3) of the relation of § log (V/R), 
which is a measure of the local reddening of the continuous spectrum, to the 
secant of the Sun’s zenith distance. It will be seen that there is considerable 
scatter, showing that zenith distance is not a sufficient criterion for atmospheric 
reddening. Further, the scatter is greater and the reddening generally stronger 
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for the observations before midsummer (open circles, Fig. 3). This seems to be an 
indication of the greater water droplet content of the atmosphere during spring and 
early summer which is a familiar feature of the Edinburgh climate. 

Before discussing the results. comment should be made on the use of a 
laboratory source as wave-length standard. The geocentric radial velocity of the 
Sun has the maximum absolute value e(1 — e?)-!* times the mean orbital velocity, 
or o:50km/sec, attained in April and October. Also the Earth’s horizontal 
rotation velocity at the latitude of Edinburgh is 0-26 km/sec. Hence the maximum 
relative velocity between laboratory source and Sun, at sunset in April or sunrise in 
October, is +0°76km/sec, corresponding to a Doppler shift of o-o1 angstrom, or 
I micron in the focal plane. This would scarcely be detectable, and is not 
revealed by any systematic trend in the readings. Thus the use of a terrestrial 
wave-length standard is justified in this application. 

4. Observations following flares, 1951.—After declining in the latter half of 1950, 
solar activity rose again during 1951, but the number of flares was of course much 
lower than during the maximum years 1947-9. ‘There was considerable geo- 
magnetic activity, in some cases associated with flares, but mainly of the recurrent 
type with indefinite onset (3). Few opportunities therefore arose for detecting 
corpuscular emission from intense flares during the period 1951 June-December 
covered by these observations. 

TABLE | 
K-line observations following flares, 1951 
Delayed 

Date Class r/R Interval geomagnetic Observed 

effect absorption 

1951 , 

184 98h 28m 0:22(E) ? SC June 18923"14™ No absorption at 14". 

194 132 o5m 0°31(W) e: ae Absorption <o°6 per 
cent 5™ after flare, 
zero after 1" (see 
Fig. 4). 

224 yo 4om 0:27(E) i uke Absorption <0°3 per 

13% 23m 0:89(W) sis se cent 35™ after first 
flare, <(0°5 per cent 
2 after second flare. 

39 132 15m 0°62(E) 55" SC Sept. 5%20% 45" No absorption on 
Sept. 4. 

74 roh 54m 0°36(W) we os No absorption 535 
after flare. 

147 132 43m 0°55(E) a idl Slightly negative 
readings preceding 
flare ! 

r54rshyom 2- o:12(W) eis ie Absorption <o-6 per 
cent 35™ after flare, 
zero after 2" (see 
Fig. 4). 

94 08h 4gm dy sae Absorption <o-8 per 
cent 50™ after flare, 
zero within 2" (see 
Fig. 4). 

174 yh 26m o-78(W) —e sia No absorption _ be- 
tween 22 and 4h 
after flare. 





M. }. Smyth, Photoelectric investigations 





' oasis C 
flere 2° 95! September |5 


—_- 


’ 
- Spaces sso + ages 


195! October 9 











Fic. 4.—K-line wing residuals following flares, 1951 (negative indicates absorption). 


TABLE 


27-day analysis of geomagnetism and K-line residuals for six solar rotations, 1951. For 
the second gives the daily mean violet-wing residual X and the third the daily mean red-wing 


1951 June 15 II 8S 4 6 143 4 6 3 2 2 6 
eee eee +2 +3 +23 wn FE. oF “— TS — 
+S +3 +1 avs o +2 .. +I 


July 12 10 II 6 6 12 


Aug. 


Sept. 


Averages 
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A few flares* followed by K-line observations are listed in Table I, together 
with the distance r/R of the flare from the disk centre, any subsequent geomagnetic 
disturbance}, and the results of the K-line observations. Only one flare of 
importance 3 occurred during the period, and observations at Edinburgh were 
unfortunately not possible until the following day. Flares of importance 1 are 
included in the table only when confirmed by ionospheric effects. 


Geomacnetic _'c' -_ a 
tsps ge yt ijn ST Pond Ree 
weaat v ' 


os 





K- LINE: REO WING 


2 Og POA gee - -0-0- - SFEHR- 











7 
Fic. 5.—Superposed-epoch analysis of geomagnetic character-figures and K-line wing residuals 
(negative indicating absorption) for six solar rotations, 1951 (see also Table II). 


II 
each rotation the first line gives the magnetic character-figures (with decimal point omitted), 
residual Y (in units of o™-o01, negative figure denoting absorption). 


8 8 + f . mw ee 6 ” 4 5 
+r +1 gic ae os we 


. 8 


Q°'2 I1'O 14°2 14°2 12°3 II‘2 97 7:0 80 
+3 eo -§ -F FE =e ~5 o -3 -!I 
+1 +1 -—-2 —2 +1 -—2 o |.-!I o -—2 


* Flare data from the Quarterly Bulletin of Solar Activity. 
t+ Geomagnetic data from the Journal of Geophysical Research. 
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It will be seen that there is little evidence of absorption. In Fig. 4, the detailed 
graphs for the three most important cases where observations were obtained for 
some time following flares, the violet-wing residuals X do perhaps show some 
tendency towards low values for an hour or so after the flares, but the effect is not 
strongly marked and is certainly less than one percent. The red-wing residuals Y 
show somewhat less variability. None of these flares preceded a magnetic storm. 

The year’s greatest magnetic storm, that of 1951 October 28 with range of 5457 
in H at Abinger (3), had the characteristic SC onset associated with solar flare 
origin, and stands completely isolated on a Bartels diagram. Nevertheless no 
flare or ionospheric disturbance had been recorded for nine days preceding the 
storm. <K-line observations on October 26 showed no trace of absorption. 

In 1953 solar activity reached minimum values, with long sequences of spotless 
days. Flares were so infrequent and small that no purpose would be served by 


TABLE 

27-day analysis of geomagnetism and K-line residuals for nine solar rotations, 1953. For 

the second gives the daily mean violet-wing residual X and the third the daily mean red- 
1953 Apr. 21 a 10 10 5 2 5 6 2 I 


4 5 $e — eee 4 2 eee I 
—6 —s5 ae ha. eck ae ca 


4 
+8 


“=—9 
8 


Ta 
Tr 


Average 
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attempting to correlate them with the present observations. Magnetic storms 
were small and of the recurrent variety (4). 

5. Analysis of recurrence tendencies, 1951.—Since much of the magnetic distur- 
bance during 1951 was of the recurrent type, but still presumably of corpuscular 
origin, it was important to look for any recurrence tendencies in the observed 
K-line residuals. It had been noted that the X (violet-wing) residuals showed 
rather greater variability, both during the day and from day to day, than did those 
for the red wing. This was most marked during 1951 September and October, a 
period of numerous prolonged and overlapping magnetic storms with indefinite 
onset, but there was no clear relationship between the magnetic activity and these 
variations. 

The daily means of the X and Y residuals were used in a superposed-epoch 
analysis in order to detect any 27-day recurrence tendency. Instrumental 


III 
each rotation the first line gives the magnetic character-figures (with decimal point omitted), 
wing residual Y (in units of o™-001, negative figure denoting absorption). 

re Yi II fe) 6 3 ° ° I sw 


+3 2 2-—!I ee Gale mie as Dae 
—IrI -—-9 -—7 8 ae 


16 % 
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stability was important in assessing long-term variations ; in the present method of 
spectrophotometry the principal cause of drift is gradual changes in the trans- 
mission of the exit-slits due to the accumulation of dust. It was satisfactory that 
the daily means of the Y (red-wing) residuals in 1951 were constant within 0-4 per 
cent for five months. The X (violet-wing) residuals showed greater variability 
and a systematic downward trend not exceeding one per cent over the same period. 
This would not be sufficient to obscure any genuine 27-day periodicity. 

Table II gives details of daily mean residuals and magnetic character-figures for 
six solar rotations in the latter half of 1951. The averages, represented graphically 
in Fig. 5, show a broad peak in the character-figures, indicating a rather poorly 
defined M-sequence. The average red-wing residuals do not depart significantly 
from zero. ‘Those for the violet wing show a slight drop a day or two after the 
peak of the magnetic disturbances. This cannot be regarded as significant. The 
smoothing of the averages for the K-line residuals is evidently poorer than for the 
magnetic character-figures, since the observations represent only 48 out of the 
possible 162 days. ‘Thus an upper limit of 0-5 per cent can certainly be set for the 
average absorption, associated with this M-sequence, by Cali ions with approach 
velocities between 750 and 1150 km/sec. 








L 











Fic. 6.—Superposed-epoch analysis of geomagnetic character-figures and K-line wing residuals 
(negative indicating absorption) for nine solar rotations, 1953 (see also Table III). 


6. Analysis of recurrence tendencies, 1953.—In 1953 the very low sunspot and 
flare activity were accompanied by two pronounced M-sequences. ‘The second of 
these, persisting for nine rotations from May until December, was well covered by 
K-line observations. A similar superposed-epoch analysis, of which details 
appear in Table III, shows a pronounced peak in the magnetic character-figures. 
It will be seen that in 1953 the red-wing residuals showed considerably greater 
variability than those for the violet wing. This can scarcely be regarded as other 
than accidental, but is nevertheless surprising. The averaged curves (Fig. 6) 
show no relation to the magnetic curve. Again the K-line residuals are less fully 
smoothed since they represent only 76 out of the possible 243 days. There is 
evidently no trace of absorption in the violet wing connected with the M-sequence. 
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7. Conclusion.—Flare observations presented in this paper support the con- 
clusion of Paper I, that there is no detectable absorption of as much as one-half per 
cent in the violet wing of the A-line following moderately intense flares. This 
result, as was pointed out, is not in conflict with the minimum particle density 
required for the generation of a magnetic storm. 

The most favourable combination of circumstances, a long series of 
observations following an intense central flare producing a great magnetic 
storm, has however not arisen during the course of the investigations described in 
these papers. A spectrophotometric technique has been developed which may be 
capable of detecting a borderline effect in such favourable circumstances, and 
should be worth applying at the next sunspot maximum. It may be that corpus- 
cular absorption rises above the threshold of observability only for a limited time 
and over a limited area of the Sun’s disk, and some experiments have been made 
with a disk-scanning device for detecting localized A-line asymmetry. 

It now appears that there is little hope of detecting absorption associated with 
M-regions while spectrophotometricerrors remain of the order of a fewthousandths 
of a magnitude. This is not surprising, since it is well known that recurrent 
magnetic storms are, asa class, of lower intensity than those associated with flares. 
It would be desirable inany further investigation to check the exit-slit transmissions 
frequently, and to attempt the most precise wave-length setting, probably involving 
temperature control of the spectrograph, in order to ensure optical and mechanical 
stability consistent with the high precision of the photomultipliers. In addition, 
the present observations have investigated absorption only by particles in the 
velocity range 750-1150 km/sec. It has been suggested (5) that “‘ slow” corpus- 
cular radiation in the range 350-600 km/sec is responsible for M-sequences, and 
this would involve absorption closer to the centre of the A-line. In any case a 
more complete survey of variations in the profile will eventually be desirable. 
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RELAXATION OF STELLAR VELOCITIES 
R. v. d. R. Woolley 


(Received 1954 June 17) 


Summary 
The stability of truncated Maxwellian velocity distributions is examined by 
calculating the fraction of stellar encounters in which the exchange of energy 
lies within a given interval. 





1. In a previous paper* the writer examined the distribution of stars in a 
globular cluster, supposing the velocity distribution to be “truncated Max- 
wellian’’, that is to say a Maxwellian distribution with all velocities in excess of 
a certain value absent. This paper attempts to enquire whether such a distri- 
bution, if set up, would have a long life without substantial change. 

The exchange of energy between two stars in a gravitational encounter 
may be calculated in a well-known manner. ‘The calculation is repeated here 
because the result is wanted in a particular form. There is no exchange of energy 
when the velocities are measured relative to the centre of gravity of the star pair, 
but when that centre of gravity moves relative to an origin (such as the centre 
of a globular cluster) there is, in general, an exchange of energy. 

Consider an encounter between two particles of which the first has mass 
m, and scalar velocity c, while the second (the field star) has mass m, and 
scalar velocity x. Let the velocity of the centre of gravity be g and the initial 
relative velocity between the stars be v. Further, let »,=m,/(m,+m,) and 
/,g=m,/(m,+m,). Fig. 1 represents the initial velocities. If the angle XOC is 0 
and the angle OGC is 8 we have 


(1.1) 


2gv sin B= 2xc sin@, 
2gv cos P=x*(1+p,—p,)-C(1+m,—p,) + 2xe(u,—p,) cosO. 





Fic. 1. 





* M.N., 114, 191, 1954. 
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The plane of Fig. 1 is the plane containing the two vector motions relative to the 
origin and is not in general the “orbital” plane which is the plane containing 
the two stars and their relative motion. Fig. 2 represents the celestial sphere 
on which the great circle C’XBC is the orbital plane and the great circle XOC is 
the plane of Fig. 1. The two points C and X, 180° apart, are the initial velocities 
of the stars relative to the motion of the centre of gravity and the point O is the 
direction of motion of the centre of the cluster relative to the centre of gravity 
of the star pair. The effect of the encounter is to move the velocity of the C star 


Fic. 2. 


in the orbital plane from C to C’, where XC’=2y, the angle between the 
asymptotes of the hyperbolic orbit. 
We have OC=f, as in Fig. 1. Let OC’=’. Now 


c2=g" + ,20? —2u gv cos B, 
and if c’ is the scalar velocity after the encounter, 
c’2 = 92 + y? wv? — 2 ,gvcos Pf’, 


since the scalar values of g and wv are the same after the encounter as before it. 
Hence if the gain of energy by the C star is AE, 


AE = }m(c'? —c*) =m,u,gv (cos B’ — cos f). (1.2) 


Now let © be the angle OCB in Fig. 2 and drop OB perpendicular to XBC. 
Let OB=i and BC=¢. Then 


cos B=cos¢ cos, 
cos B’ = cos (7 + 2s — d) cost, (1.3) 


sin 8 cos © =sin¢ cosi. 
Combining (1.1), (1.2) and (1.3) we get 
AE = m,p1,[Cos? p{x(1 + pg — fhe) — (1 + pe — Mg) + 2KO(M, — fg) COS B} 
+ cos siny . 2xc sin @ cos O}. (1.4) 
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Now it can be shown* that if D is the distance at which the stars would pass if 
there was no deflection, and if I is the constant of gravitation, 


cos? ys = {1 + D*v4/T2(m, + m,)*}", 


sin ys cosy = {Dv?/T'(m,+m,)}/{1 + D®v4/T2(m, + m,)*}. 
Write 
D*v4/T?(m, + m,)? =r'v4. 
Then from (1.4) 
( I+ rv*)AE, Mf = x°(1 +he— Me) _ c*( I+p.— Me) + 2xc(u, ra Ma) cos # 


+ ryt. 2xc sin 8 cos @. 


Again, v? =x? +c?—2xc cos 6, and using this and writing AE/u.,m,=2, we get 


a= {1+r4(x? +c? —2xc cos 6)*} x1 +,—p,.)-—C(1 +p, —bMe) 
+ 2x¢(,—p,) Cos A + 72(x? + c?—2xccos@).2xcsinOcos@]. (1.5) 


Equation (1.5) shows the energy exchange z in a form which is useful because 
the velocities which appear, x and c, are scalar velocities relative to the origin, 
while the angles @ and © have simple distributions. It is assumed that at any 
point in the cluster (or other aggregation of stars with which we wish to deal) 
the distribution of velocities is isotropic. ‘Then the fraction of the encounters 
in which @ lies between @ and @+ dé@ is } sin@ dé and the fraction in which @ lies 
between © and ©+d0 is dO/z. We also suppose that the fraction of the 
encounters in which D lies between D and D+dD is 2DdD/D,?. The nature 
of the limit D, is discussed at length by Chandrasekhar in his Principles of Stellar 
Dynamics, pp. 55-57. We shall follow Chandrasekhar and suppose that D, is 
the average distance between stars. 
2. If we write 
r(x? + c? — 2xe cos 0) =y, 


x(1+p,—p,)-C(1+p.—p,) + 2xe(u, —,) cosd=a, 


2xc sin#=5, 
then from (1.5) 


at+bcos®.y 

z= ~~ (2.1) 
and, for a given x, c and 6, the quantities a and 5 are constants. We now enquire 
in what fraction of the encounter does z lie between k and k + dk when x, c and 0 
have fixed values, while all values of © between o and a are equally likely. 
The method of calculating this is illustrated in Fig. 3. Let the solution of 
f(x, y)=h be y=4(x, h). In the figure the curve ABCD represents f(x, y)=h+ dh 
and the curve EFGH represents f(x, y)=h. The variables x and y may have 
any values between zero and X and between zero and Y respectively, and the 
fraction of cases in which x lies between x and x + dx is supposed to be p,(x) dx/X, 
with a similar supposition about y. Now the fraction of cases in which f(x, y) 
lies between h and h+dh is the sum of the fractions of cases in which it lies in 

areas of which BCGF is an example, i.e. 


¥FG.p,(OL). BFp,(OP)/XY. 





* e.g. S, Chandrasekhar, Principles of Stellar Dynamics, Appendix 1, Chicago, 1942. 
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But if FG =dx, then 
BF = ar |(A” * an) 
and the required ere is 


i f pal). py(b(x, A). dx (he ) ; 














Fic. 3. 


We apply this technique to calculate p,(k)dk, the fraction of the encounters 

(having a given @) in which z lies between k and dk. All values of © between 

o and = are equally likely, while y lies between o and Y, the density being 
Py(¥) = 29/¥ 

Before applying (2.2) we examine the limits involved. If z= then 


cos © = {k(1 + y") — a}/by (2.3) 


2ky = bcos © + {b? cos? © — 4k(k —a)}*". 
No real solution for © exists unless b?>4k(k—a), and for @=o or 7, 
2ky= +b + {b®—4k(k—a)}*?, 
so that the curve (2.3) cuts the boundaries cos@= +1 and cos@= —1 at four 
points, in two of which y is positive and two negative (unless k=o). When 
k(k—a) is positive the two positive points are 


I }' 
Ya=|b/2k| + | (04h a)}9|, 


I 
yi=|b/2k| — [5 (0 4k(k— ay], 


and both of these points lie on cos @ = +1 or both on cos@= —1. When k(k—a) 
is negative the two positive points are 


Ya=  |b/2k|+| = (08 4h(k—a)}?? |, 


yi= —[b/ak| +] (0 4e(k—a)}*® |, 
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and one point lies on cos@= +1 and the other on cos@=-—1. In all cases 
Vey: 
Equation (2.2) gives 
‘Ll 2y dy 
p(k) = »,7Y® |(6z/00), | 
and the upper limit L is y, or Y, whichever is smaller. Since 
by sin © 


| (dz/6@),, | = . 


we have 


2 


~ (14+?) dy 
p{k)= =>] Sot 


y, [osinO] ° 
When L=yz, i.e. when 
b I 
i food nie 
" ahk| * | 2k 
the integration leads simply to 


{6*— 4k(k—a)}*|, (2.4) 


5b? + 2ak 
p(k) = 2y2 R83 | ’ 


whereas, if (2.4) is not satisfied, 


(2.5) 


1 |b%+2ak 2/Y?+n I 
ee -1 —S 72 4 1/2 
p{k)= | > 6S (mi? — gn) pls +mY* +n) . (2.6) 
in which /= —k?, m=b?—2k(k—a) and n= —(a—k)*. The inequality (2.4) 
may also be stated in the form 





a+by 
1+ Y? 


|k| > . (2. 

3. The number of encounters which occur in time dt between a star with 
velocity c and n(x) dx stars per unit volume with velocities between x and x +dx 
is 


“Dy “71 
n(x) dx dt | 27D dD | (x? +c? — 2xc cos 0)!” } sin @ dé 
~0 ~0 


x? +c? ae 
= n(x) dx dt .rD,? x 3 ~— ” ae (3-11) 
3: 
x? + 2c? os 
= n(x) dx dt .1D,? x — fw #<c. (3.12) 
3C 
Similarly the number of encounters N(k, x, c)dxdkdt with energy exchanges 
between k and k+dk is given by 


N(k, x, c) =n(x).7D,? | i sin@.(x®+c?—2xccos@)"*.p(k).d0. (3-13) 
~ 0 
Since a and 4 are functions of 0, p,(k) is a function of @ (and k). For simplicity 
we take, from now on, the case where the stars have equal mass, »,=y,. Then 
a=x*—¢?, b=2xcsiné@, 
y=r(x2+c?—2xccos0), Y=R*(x*+c?—2xccos6), 
R* = D,/20m. 
— at+by 
The condition k= Ta yz becomes 


(x? —c®) + axc sin @. R(x? +c? — 2xc cos @) 
1 + R4(x? + c? — axe cos 6)? 





k= 
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When R is large and x <c this is 
k 


+ 2xc sind 
~ 
R*(x* +c? — axc cos 6) 





2XC 


; , tS Ee ee. ey 
which has maximal values where cos @=2xc/(x*+c*) and k= + Rt) 


, SO 


that if |k| > Ras equation (2.5) applies. When x=c and sin@=o 


the energy exchange is zero. 
Using 2I'mR* = D, equation (3.13) becomes 
N(k, x, c) = 47T*m?R'n(x) . | } sin (x? + c? — 2xc cos 6)". p(k) dé. 
0 
Let 


(x) 2xc 
x F(k,x,c) when Vel >| paar (3.21) 


27T?m*n 
[k* | 
27*T'm?n(x) 


| ®* | 


N(k, x, €) = 


2XxC 
Rx2—2)|" (3.22) 


x G(k,x,c) when |k] <| 


The values of F(k, x,c) may be given simply. When |k| > RE c)|’ 
bi I 4x°c? sin? 6 — 2(x? —c?)k 
p(R) =; 2R4 ks | = (x? +c? — 2xc cos 6)? 


=o when sin? @<sin?4, 


when sin? @>sin?4, 





and xc? sin?6,=k?—k(x*—c?), with the conventions 0<@,<7/2; 6,=o if 
k? — k(x? —c?) is negative; and 6,=7/2 if x?>k and k is positive or if c?>|k| 
and k is negative. ‘These limitations arise through (2.7). Carrying out the 
integration we have 
k negative.—|k | must be less than c? as the ¢ star cannot give up more energy 
than it has. 
From x=o to x*=c?—|k] 
8x” — 6k 
F(k, x, c)= oo (3.31) 
With x? >c?—|k] 
: 8c? + 2k 
F(k, x, c) = (c2 + k)¥?. (3-32) 
3xc 
k positive.—F(k, x,c) is zero when k>x* as the x star cannot give up more 
energy than it has. 
From x?=o to x2=k 


F(k, x, c)=0. (3-33) 


From x?=k to x27=c?+k 
8x? — 2k ; 
id 2__ py 
F(k, x50) = = (at) (3-34) 
With x2 >c?+k 
a 8c? + 6k 
F(k, x, c)= oe (3-35) 
Detailed balancing.—Let the distribution functions for c and x be c®f(c) and 
x°f(x), that is to say, let n(x)=.x*f(x) and let the c stars be the normal field stars. 


35* 
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Then the number of collisions which destroy two velocities c and x and set up 
two velocities c’ and x’ (with energy changes k= c’? —c? = x? — x”*) is 
I 

[| 
per unit area per unit time, when |k| > |2xc/R%(x*—c?)|. These arise out of 
a set of orbits all of which can be described in the opposite sense, destroying 
c’ and x’ and setting up c and x, with the energy change kh’ =c?—c’?=x’*—x?. 
The number is 


c'2x'2f(c’) dc'f(x') dx’ . 


c?xf(x) def(x) dx x . F(k, x, c) dk x 2nT*m? 


I 


Lane 





F(k’, x’, c’) dk’ x 2nT?m?. 


The system is conservative if 
c®x®f(c) f(x) . de dx . F(k, x, c) dk =c'*x"*f(c’) f(x’) de’ dx’ . F( —k, x’, c’) | dk’ | 
with dk=c' dc’ = —x' dx’, dk'= —cdc=xdx and dk= —dk’. Hence 
xe dc dx dk =x'c' dc’ dx’ | dk’ | 

and the system is conservative if 

wef(c)f(x)F(k, x, c)=x'c'flc’)flx’) . F(—k, x’, ¢’). (3.4) 
Consider for example k positive and x?>c?+k. Then F(k, x, c)=(8c? + 6k)/3x. 
by (3.35). In the reverse exchange k is negative and x’*>c’®—k so that 


2 
F(—hy2', c= et by (3:52) 
3x'c 
8c? + 6k 
= 77 xc 
3x Cc 





XC 
= F(k, x,c) x WO 


Hence (3.4) gives KOfl*) =fle lx’), 
which is of course satisfied by Maxwell’s distribution f(y)=e"*””, since 
2+ x2 =¢'2 + x2, 
Detailed balancing is therefore satisfied. It may be verified in the other cases 


that 
xcF(k, x, c)=x'c' F(—k, x’, c’), 


provided c’?=c?+k and x"*=x?—k. It is clear that we must also have 
xeG(k, x, c)=x'c’G(—k, x’, c’). 
4. Now suppose that the velocity distribution is truncated Maxwellian, or 


n(x)=N’x’e-?=* when x<xq 
=o when x>xy (4-1) 


and n= /{™ xe?" dy, 
0 


Let jx)=J. We suppose that e~’" is a small number, so that N’ = 4j37"!2N, 

The number of stars removed from the class c, de by collisions x, dx : k, dk 
in unit volume in time dt is 
dk p i 
n(c) de dt n(x) dx TF | x 2al*m?H(k, x, c), (4.2) 
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H(k,x,c)=F(k,x,c) if k>k, 
=G(k,x,c) if k<k, 


and ky = | 2xc/ R*(x? —c?) |. 
In full equilibrium this removal is exactly balanced by the gain to ¢, de 


n(c’) de'n(x’) dx’ Fy 2nl*m?H(—k, x’, c’), (4.3) 
where c’dc'=cdc, x'dx'=xdx, c*=c?+k and x'*=x*—k and k may be 
positive or negative. 

Under condition (4.1) certain unbalances occur. The net effect on the 
population n(c) can be calculated by counting those transitions which do not 
occur as a result of the cut-off. These result in a gain, or no change, or a loss, 
to the population n(c), as the case may be, as follows :— 

If all four states in the transition c, c’ : x, x’ are of lower energy than the 
cut-off x), then the c and c’ transitions both occur, and balance, giving no net 
change in n(c). If two or more states are of greater energy than x, the process 
and its reversal are both cut off, giving no net change in n(c). 

If however one only of the four states c, c’ : x, x’ is of greater energy than 
the cut-off, a net change does occur. Since c?>x,? is forbidden ex hypothesi 
there are three cases to be considered, namely (1) x? >x,?, (2) x’? >x,? and 
(3) c’?>x92. In case (1) the transition but not its reversal is forbidden, so that 
there is a net gain to n(c). In other cases the transition is permitted but the 
reversals are suppressed so that there are net losses to n(c). 

(1) When x? > x,? but both c’?<x,? and x’*<x,?, k must be positive and 


Xp? - C2 >k>x? —x,?. 
Let k=+/*. Then the loss to n(c) dc is 


— anP%mn(e') de’ at | " n(x’) H(-P,x',c’)dx’, (4.41) 
(x*,—I*) 


the negative sign in H(—/?,x’,c’) occurring because the unbalanced transitions 
are losses to c’.. Equation (4.41) holds from /?=0 to 2 =x,?—c*. 

(2) When x’?>x,? but both c’?<x,? and x?<x,?, k is negative, equal to 
—/*, say. Then 


Xp? > x? > x2 -—P? 
and the loss to n(c) dc is 


2ldl 7 
2m1*m?n(c) de dt ae | n(x)H( —I?, x, c) dx (4.42) 
(a*,—l*)* 
valid from /?=0 to /?=c?. 
(3) When c’?>x,?, k is positive and greater than x)?—c?. Then x? >x,.?—¢? 
as the x star cannot give up more energy than it has. The loss to n(c)dc is 
2m1m?n(c) de dt = n(x)H(I?, x, c) dx (4.43) 
~ (a*—c*)* 
valid from /? =x,?—c? to ?=.x,?. 
Case (3) is unimportant except when c* is nearly equal to x,’. 
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5. Accurate evaluation of the function G(R, x, c) is a matter of great complexity. 
We shall assume that (when x >c) 


O-Aisne 8c? (1+a3)/ 


3x + o3/,°’ (5-1) 
a form which makes G= F when /=/, and w aie keeps GI-* finite when/=o. The 
value of « may be chosen to make the total seis - collisions 


ex? 


aT %m?n(x) 2 H(k, x, €) a5 : _ (5.2) 





equal to the number of collisions given (accurately) by (3.11) or, 

n(x) .7Do*(3x* + c)/3x. (5-3) 
Substituting (5.1) in (5.2) and making approximations suitable to x)? >/,? and 
c?>/,? the number of collisions given by (5.2) is 


* 8c? (1 bk (r+a) | 4 He 
or 1 Resi ade an Nes 
271 ?m?n(x) 3x {7 (1+ [ 4m I+a3 V3 (tan a3 ag =) |) 


and equating this to (5.3) with pe Re, and I 2xc/ R(x? —c*) we find 


1+a3 (1+«)8 2 2-a 7 3x27+c2 = 8xc? 
i+— } In - + —- (tar *—- +7) /= wa ere 
on? I+% V3 a3 6 8c (x* — c*)* 


2 3 
This value of « makes the total number of collisions given by the approximate 
formula (5.1) correct. 
The case in which we are most interested is x>c, and then 


a= (4x%/7e2)"" 
approximately : « therefore takes large values. ‘This has been established to 
show that the terms involving « in equation (6.3) can be neglected. 
6. From equations (4.41), (4.42) and (4.43) 



































ri dn(c) _ en(c’) (vm aldl a 
2ml@m? n(c) dt c'n(c) Jo _ | ee n(x’) H( —1?, x’, c’) dx 
ais {" n(x)H(—I?, x, c) dx 
i V(ae*—I*) 
ts 2ldl (7 
f | (a_2—c?) | e,2—c2) n(x)H(I?, x, c) dx. (6.1) 


It is now necessary to approximate. If x)?—c? is not small, the third integral 
may be neglected in comparison with the other two, since /* occurs in all 
denominators. We further take x )?—c?>c* and neglect the range from /=c 
to 1= \/(x,?—c?) in the first integral, again because small values of / contribute 
most to the integrals. We are left with 


I 1 dn(c)_ (¢2ldl cn(c dy : 
eae Ses. eto ahi _ 
(6.2) 


Now en(c’)/c'n(c) =c’ /c exp {j*(c? — c’?)} =(1 + /c?)"e-*"", and we can approximate 
to H(—I?, x, c) by 





H(-F*, x, c)= fll), 
where f()=1 when />h, 


[6 
A= 5 e ee . if 1<Jh, following Section 5. 
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Then 


_1 dic) 167 129,27’ 2 x a as —7t f —i*a* 
eae ewe [et a) el a, ee 


32 al ?msNe'(1 — § °c?) [mos 


rn *9n7_ J? e: Cc " I+a3 
= 164/nl*mj3Ne* (1 — §j%c?) | In q2t $(1 +3) In = _ ae 
The terms involving « are unimportant, with « large, since c?/J,? is large. This 
justifies the approximation in Section 5. In what follows, these terms will be 
neglected. 


We now define the time of relaxation 7(c, x») by 
1 dn(c) I 


n(c) dt ~~ T(c,%)° 
Then, substituting }R*x yc for c?//,? (supposing x) >c) we have finally 
[T(c, xp)} 1 = 16 y/aP2m43N( 3%? — 1) In (J RExyc)e” (6.4) 
This may be compared with Chandrasekhar’s time of relaxation* 7, which is 
given, in our notation, by 
a 321 7) pemj2N In (3R*/2/%), (6.5) 
7. For small values of c, namely j*c?< 3, the time of relaxation due to the 
Maxwellian truncation is negative, i.e. the number of stars with small velocities 
increases. Now if n(c)=427"?PceF, 
S OM Aeneas 
“ee cee 


so that equation (6.4) gives 


a= 
showing that the collisions in the truncated distribution bring about an increase 
in the quantity /, i.e. a decrease in the “temperature”. But the main difference 


; 6 
f= OX ATA In (IR Ce 


between equations (6.4) and (6.5) is the appearance of the term e~”’ in the former. 


Accordingly if e~’* is a small number, that is to say if the cut occurs far out in 
the Maxwellian distribution, the inverse time of relaxation is much smaller, 
and therefore the time of relaxation much larger, than the quantities given by 
Chandrasekhar. If then his time of relaxation is of the order of magnitude of 
the time required to set up, substantially, a Maxwellian velocity distribution, 
this distribution will be stable over a much larger period of time if the escape 
of fast stars cuts off the high velocities. ‘This is the point which the present paper 
attempts to establish, and the work is offered as a justification of the “ model” 
velocity distribution taken in the writer’s earlier paper. 

8. Acknowledgment.—I am indebted to Mr J. Burr who has checked nearly 
all the algebra shown in this paper. 
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* §. Chandrasekhar, Principles of Stellar Dynamics, equation (2.379), p. 67, Chicago, 1942. 








TEMPERATURE AND TURBULENCE IN THE CHROMOSPHERE 


R. O. Redman and Z. Suemoto 
(Received 1954 September 4) 


Summary 


Measures of chromospheric line profiles in high resolution spectrograms 
obtained at the total solar eclipse of 1952 February 25 have been used to 
investigate further the question of temperature and turbulence in the solar 
chromosphere. In trying to deduce a distribution of line-of-sight velocities, 
arising from either temperature motions or turbulence, two problems are 
encountered ;: (a) to estimate the chromospheric heights to which the 
spectra refer, and (5) to correct the measures for broadening by self-absorption , 
and Stark effect. Heights have been estimated partly from a comparison of 
line intensity ratios with those measured by Dr J. Houtgast on slitless spectro- 
grams obtained at the same eclipse, and partly from the time of exposures, 
using the motion of the Moon relative to the Sun. Broadening by self- 
absorption has been estimated in the case of metallic lines by an empirical 
method. The corrected line widths, together with others of He 1, suggest an 
approximately linear increase of line-of-sight velocity with height, from 
&)~2°5 km/sec at the base of the chromosphere, to &)~16 km/sec at height 
2600 km. ‘Taking into account the uncertainty in effective height, these 
results agree with those obtained in 1940. 

The hydrogen Balmer lines have been examined in more detail. The 
earlier members of the series are broadened by self-absorption, while later 
members are broadened by Stark effect. Both kinds of broadening were 
first estimated by an empirical method, which led to the conclusion that the 
temperature of the lower chromosphere does not exceed 10 000 deg. K. The 
hydrogen line widths were then computed more elaborately by a theoretical 
method, which showed that the observations could be fitted almost as weil 
by a temperature 6 000 deg. as by 10000 deg. It seems unlikely that the 
present data will yield any closer estimate of temperature than this. 





At the total solar eclipse of 1940 October 1 line widths in the chromospheric 
spectrum were measured which seemed to indicate two important results: 
(a) the turbulence in the low chromosphere is small, the average (horizontal) 
line-of-sight velocity being about 1-8 km/sec, (6) the kinetic temperature of the 
low chromosphere is about 30000 deg. K.* While no objection has so far been 
raised to the first of these conclusions, the second has provoked a great deal of 
discussion and the general trend of opinion today is that such a high temperature 
is most unlikely. Thus there has been a call for more work, first to re-measure 
the line profiles to ensure that the 1940 measures contain no substantial error, 
and then, if they turn out to be correct, to seek some alternative explanation 
consistent with a temperature not far from that of the photosphere. 

With this in mind, observations were made at the total eclipse of 1952 
February 25, in Khartoum, with a large grating spectrograph. A preliminary 


* R. O. Redman, M.N. 102, 140, 1942 (Paper i). 
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description of this work has already been given elsewhere* and further details 
are now in the press.t| The more essential points may be repeated here. The 
spectrograph used the second order of a 6-inch Rowland concave grating, having 
15000 lines per inch, in a Wadsworth mounting. The collimator mirror had 
a focal length 12 feet, the nominal grating radius was 21 feet. For thermal 
protection the instrument was installed in a pit in the ground. The slit was set 
tangential to the solar image near the point either of second or third contact. 
The solar image had a diameter 52 mm; the slit width was 0-03 mm; scintillation 
was estimated to give an image spread roughly 2 seconds of arc. Eight photographs 
were taken at second contact covering the region 5550-6040 A, with dispersion 
2'2A;mm, and seven at third contact covering 3430-4110A, with dispersion 
2-4A/mm. The average exposure time was about 1-5 seconds, with an interval 
between the centres of successive exposures about 2-3 seconds. The plates were 
calibrated with a Hilger rhodium step wedge, and the instrumental profile was 
determined by means of a water-cooled *Hg isotope lamp. The spectrograph 
gives very satisfactory definition. The sum of the Rowland ghosts in the second 
order is about one per cent of the parent line. The general photometric procedure 
was in principle the same as in 1940, although almost all the instrumental details 
were different. 

Results from narrow lines are given in Paper III and confirm quite closely 
the 1940 measures. The assumption that for these lines the observed width, 
after correction for instrumental broadening, is chiefly due to line-of-sight motion, 
is supported by the fact that the two spectrum regions yield nearly the same mean 
velocities. Following Unséld’s notation], the velocities € are taken to have a 
distribution given by exp —(€/&,)*, so that the results may be expressed in terms 
of the parameter &. ‘The 1940 measures (Paper I) were given in terms of 
the arithmetical mean velocity 6. The two parameters are connected by 
£) = 54/7. 

In this paper we shall merely quote the results for narrow lines and then 
concern ourselves mainly with the problem of wider lines, very many of which 
are known to be broadened by self-absorption. ‘Two results will emerge: 
(a) evidence is found which suggests a considerable increase of velocity with 
chromospheric height, (b) the kinetic temperature in the low chromosphere 
does not exceed 10000 deg. K, while 6000 deg. K may almost equally well be 
reconciled with the line widths. At present it appears doubtful whether a more 
accurate kinetic temperature can be deduced from the available line profiles. 

Heights.—The practical problem of deducing a velocity distribution from 
chromospheric line profiles leads to two other questions: (1) the determination 
of effective heights, (2) the correction of the measured line widths for various 
distortions, by the instrument, by self-absorption, and in the case of hydrogen 
by Stark effect. 

The assignment of a chromospheric height to observations made with a slit 
spectrograph is a difficult question to which no entirely satisfactory answer has 
been found. As explained in Paper III, the only tolerably satisfactory way of 
estimating the heights at which our ultra-violet (third contact) observations 
were made is to compare line intensity ratios with those from slitless spectra. 

*R. O. Redman, Accademia Nazionale dei Lincei, Convegno Volta, 11, 72, 1953 (Paper II). 


ft R. O, Redman, Vistas in Astronomy, in press, 1955 (Paper III). 
A. Unséld, Z. f. Naturforschung, Ja, 121, 1952. 
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What we aim at getting is only an apparent height, measured from the solar 
limb, and we try to express the physical properties of the chromosphere in terms 
of what would be seen, integrated along the line of sight, in ideal conditions. 
What ultimately we require are the physical properties in terms of true heights 
within the chromosphere, but the derivation of these can be regarded as a separate 
problem, which is not tackled here. 

The apparent heights deduced from line intensity ratios are not very precise 
for at least two reasons: (a) What the s/it spectrograph records is an indeterminate 
blend of spectra from an appreciable range in chromospheric height, the degree of 
mixing depending chiefly on scintillation, guiding and length of exposure (the 
coelostat having been rated to follow the Moon), (5) the structure of many of the 
spectrum lines changes greatly with height, and what slitless spectrographs 
record is the sum of all the chromosphere not covered by the Moon at the time of 
exposure. Ifthe seeing is very good and the instrumental resolution high enough, 
these various line structures tend to be recorded side by side; if the scintillation 
is great they are all blurred together. 

However, despite various difficulties, this method of height determination, 
i.e. by comparison of line intensity ratios between slit and slitless spectrograms, 
appears to be the least unsatisfactory for the present case. Using data kindly 
supplied in advance of publication by Dr J. Houtgast, apparent heights have been 
estimated as follows: 


3rd contact exposure I 1200: km 
2 1100 
3 level 1 < 50 
2 <100 
3 600 


The chief cause of the small differences of effective height between exposures 1 
and 2 is probably poor guiding. This arose froma slight mishap with the coelostat 
drive, which slipped out of gear about 30 seconds before the end of totality. The 
solar image was recovered only just in time for third contact. The three heights on 
exposure 3 refer to measures made from different points of the spectrograph slit. 

For the yellow observations at second contact, heights have been estimated from 
the motion of the Moon across the Sun. The last traces of Fraunhofer spectrum 
are seen on exposure 3 and the narrow lines measured in Paper III are from 
exposure 4, ata height <sookm. Subsequent exposures show nothing but NaI, 
D, and D,, and Her Ds. For these the estimated heights were based on the 
assumption that second contact took place between exposures 3 and 4, and that 
the guiding was good. In some respects this method is more satisfactory than the 
use of relative line intensities, but when a rather slow succession of exposures is 
made, as here, it can be applied only to strong lines which extend to comparatively 
great heights. ‘The adopted heights are as follows: 


2nd contact exposure 4 200 km 
5 800 
6 1400 
7 2000 
8 2600 


Metallic lines.—In the violet and ultra-violet spectrum from low chromo- 
spheric levels, say less than 600 km, there is an adequate selection of lines, notably 
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of the singly ionized rare earths, which appear to be little broadened by self- 
reversal, and these give the line-of-sight velocities already reported and reproduced 
here in Table I. 
TaBLe I 
Line-of-sight velocities from narrow lines 
Height v & Number of lines 
ultra-violet << 50km 1°5 km/sec 2°7 km/sec 38 
<100 1°8 3°'2 54 
~ 600 2°2 3°9 27 
yellow <500 12 21 18 
1940 violet <500 18 3°1 30 
At these same comparatively low levels the strong metallic lines, especially those 
which are also strong in the Fraunhofer spectrum, are generally considerably 
broader than the narrowest lines, and often their profiles show clear indication of 
of self-reversal.* At greater heights, unfortunately, practically the only lines which 
persist in sufficient strength to be measurable are these same lines which at low 
levels are known to be affected strongly by self-reversal. True, the optical depth 
in the line of sight decreases with increasing height, so that the self-absorption 
decreases and somewhere becomes negligible, but we cannot be sure without 
further examination that this happens before the line becomes too weak to measure. 

A fairly certain answer may be obtained if there are enough lines to construct a 
curve of growth, for on the first linear part of the curve there should be negligible 
self-absorption. Construction of an empirical curve of growth from Fe I multiplets 
showed that most of the easily measurable iron lines on our u.v. spectrograms fall 
near or on the second section of the curve, so that their profiles are unlikely to give 
satisfactory line-of-sight velocity distributions without correction. It is possible 
to compute the necessary corrections on the basis of some model chromosphere, 
but here we have preferred to use an empirical method, from extrapolation of the 
observed data. 

As many lines as possible of a given element are measured at some given height, 
and their half-widths (ordinates) plotted against line intensities (abscissae). As 
expected, the widths increase with growing intensity, as the broadening by self- 
absorption increases. Extrapolation of the observed curve back to zero intensity 
should, however, give a ‘‘ true” line width, free of the effect of self-absorption. 
The extrapolated curve is assumed to start from the ordinate axis with zero slope. 
This procedure was first tried for u.v. exposure 3, at the 600 km level, where there 
are also sufficient rare-earth and other narrow lines to give a “‘true”’ half-width 
without correction. ‘The extrapolated half-widths (corrected for instrumental 
broadening) for Fe 1 and Ti 1 give reasonably accordant results : 


600 km rare earths, etc. 5 = 2:2 km/sec €) =3'9 km/sec 
extrapolated Fe 1 22 39 
extrapolated Ti 11 2-7 48 


At higher levels there were rather few lines measurable, even of Fert and Tin, 
and since exposures 1 and 2 refer to approximately the same height they were 
combined, with the following result : 


1200 km extrapolated Fe 1 5=5 km/sec €)=9 km/sec 
extrapolated Ti 1 5 9 


* Not merely self-absorption, but an actual self-reversal at the centre of the line profile. 
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The exact coincidence of the Fe and Ti velocities here should be regarded as 
fortuitous, and not a sign of great accuracy in the result. 

In the yellow the only metallic lines visible above 500 km were the D lines of 
sodium. These are very heavily self-reversed in the low chromosphere, in the 
manner typical of all low-excitation lines which are also strong in the Fraunhofer 
spectrum. However, at the highest level at which they are measurable, 1400 km, 
both D lines have the same width and we have supposed them to be free of 
serious self-absorption broadening at this height. (We note, however, that their 
intensities here are in the ratio 1: 1-4): 

1400 km D lines 5=5-3 km/sec £)=9'4km/sec 

As far as these spectrograms are concerned there are no further data to be 
obtained from metallic lines at greater heights; for such levels a faster spectro- 
graph or longer exposures would be required. 

Helium.—The only Het lines measurable on our spectrograms are 5876 A (Ds), 
4026A, and on one exposure 3819A. A few other helium lines could be expected 
in the 3500-4000 A region, but at the fairly low levels where they might otherwise be 
measurable they are concealed by blends with other lines. ‘Throughout the violet 
and ultra-violet spectrum, blending is generally severe at low levels, especially below 
100km, and owing to the intrinsic width of the lines an increase of instrumental 
resolving power can bring only partial amelioration. 

All three of the measured Het lines are multiple. D, consists of a close unre- 
solved doublet with a partially resolved fainter companion, the others contain a 
single line with a partially resolved fainter companion. In each case the fainter 
companion has, as far as possible, been removed from the measures. After this 
correction, and after correction for instrumental broadening, we get: 


Height Line 
600 km 4026 A v=6-2 km/sec o=11'0 km/sec 
1100 4026 6-6 11°7 
1100 3819 6°6 I1°7 
1400 5876 6-0 10°7 
2000 5876 7°5 13°3 
2600 5876 8-9 i ty 


On the 1940 spectrograms only the 4471 A Het line was measured, without 
correction for the faint companion, giving = 8, = 14 km/sec at some uncertain 
height probably less than 1500km. It was remarked at the time that this line 
width was probably a little too great. 

Variation of &, with height.—If now we plot all our results for €) against height, 
including those of Paper III, we find an almost linear increase from about 
§9= 22 km/sec, at the base of the chromosphere, to 16 km/sec at 2500 km (Fig. 1). 
No definite explanation can be offered for the two seriously discordant points, 
from He! 4026 and 38194, but the figures for all the He! lines have to be accepted 
with some reserve. In the first place no examination has been made of possible 
widening by self-absorption, but if comparable with self-absorption in hydrogen 
lines of the same intensity, the effect is quite weak for 4026 and 3819 A at 1100 km. 
In the case of D; nothing can be said except that the line is weakening from 
exposure 6 to 8, but at the same time the measured line width actually increases ; 
that is, if there is really an appreciable amount of self-absorption broadening (which 
must diminish with height), the velocity increase with height is still greater than 
that just deduced. 








No. 5, 1954 Temperature and turbulence in the chromosphere 529 


There is another complication in the case of helium. As already explained, the 
heights we have tried to assign to the measures are those at which the spectrograph 
slit would have to be placed to give the observed line widths in perfect observing 
conditions (no scintillation, etc). ‘These apparent heights are not the same as 
heights in the chromosphere itself. For instance, were there an exponential 
intensity distribution with height, with no self-absorption, half the scale height 
would need to be added to the apparent height to give the true mean height in the 
chromosphere. With helium the relation between intensity and height is certainly 
more complicated than this. Many observers have noted that at low levels Het 
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lines tend to behave differently from most others, in that their intensities at first 
increase with apparent height, a behaviour which is presumably closely associated 
with the general increase of spectrum excitation with height. A simple exponen- 
tial law does not represent the facts adequately. The conversion of behaviour with 
regard to apparent height to behaviour in terms of true height can be made 
sufficiently well only with slitless spectra and we do not attempt it here. Fig. 1 
ignores the correction from apparent height to true height, and we have to 
remember that the relation between these heights is probably less simple for Her 
than for any other element. 
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Thus although the evidence for an increase of £, with height is fairly convincing, 
the linear relation of Fig. 1 should not be taken too readily at face value. The 
points have been taken from a variety of lines, and the height question in particular 
cannot yet be regarded as properly settled. The evidence, taken together with 
what follows on the hydrogen lines, suggests that the velocity increase is chiefly 
one of turbulence, rather than of temperature, although more observations are 
certainly needed before there can be a really satisfactory separation of turbulent 
from thermal motions. ‘The idea of a turbulence increasing with height has been 
favoured by some theorists, especially Unséld.* 

Hydrogen.—There are no hydrogen lines on the plate exposed at second contact, 
but the third contact photographs cover the whole Balmer series from Hé to the 
series limit. H6é itself was heavily overexposed, and not measurable on any 
exposure; others of the earlier and stronger series members could be measured 
only at the higher levels. ‘The measured half-widths, i.e. the total widths at half 
the maximum intensity, are in Table II. Gaps inthe measures are due to blends, 
which are particularly troublesome at low levels, to overexposure of lower members 
of the series, or underexposure of higher members. In cases marked * the 
whole line could not be measured because of blends and the half-width was 
obtained by assuming a symmetrical line. The heights were estimated as before. 

Three specimen hydrogen profiles are given in Fig. 2, the points being single 
readings from one exposure. Lines near m= 30 are affected by the Balmer 
continuum, which is just measurable on exposure 3, at height 600km. At lower 
levels the Balmer continuum is concealed by a general continuum; at higher 
levels it is too weak to appear on these spectrograms. When finding the half- 
widths the Balmer continuum has not been removed from the measures, but has 
been regarded as a part of the lines themselves. (It becomes important only for 
lines with m>25.) ‘The behaviour of line width with 7 is more regular when the 
Balmer continuum is left in than when it is removed. 

A first inspection of the hydrogen line widths shows a tendency to increase 
as one approaches the series limit, which qualitatively is what is to be expected 
from Stark effect. In passing we may note that the lines do not overlap near the 
series limit, as far as these measures go. In the opposite direction, as one goes 
to stronger lines there is ultimately also a systematic widening, which is fairly 
certainly due to self-absorption. All lines now measured are substantially 
narrower than the widths for HS, Hy and Hé (m=4, 5, 6) on the 1940 spectro- 
grams, which tends to support the suggestion that all the earlier and stronger 
members of the Balmer series are considerably broadened by self-absorption 
when observed at fairly low chromospheric levels. 

We suppose now that the measured hydrogen line widths are due to 
(a) instrumental broadening, (+) Stark effect, (c) self-absorption, (d) Doppler effect. 
Of these (a) has already been removed and we wish chiefly to find (d). There 
is a weakness in the raw material, inevitable in the present case, and not easily 
avoided altogether in eclipse spectroscopy of this kind, that because of over- 
exposure the accuracy of measurement falls off with the strongest lines, on 
which we mostly rely for estimating the effect of self-absorption, and because of 
underexposure it falls off also with the weakest lines, from which we deduce the 
Stark effect. Ideally, the measures should cover a considerably greater range 
than is possible on a single photographic exposure. 

* A. Unsdld, loc. cit. 
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TaBLe II 
Half-widths of chromospheric hydrogen lines 
Blended lines marked* have been assumed symmetrical; other blended lines have been 
omitted. Uncertain measures are marked with a colon. The half-width is the total line 
width at half the maximum intensity. 


Total Measured h.w. corrected for 
intensity half-width instrumental distortion 


Height 100 km 1360 0°374A 0349 A 
14* 1085 ‘345 ‘321 
15 865 328 304 
16 765 333 -308 
548 ‘333 -308 
18 458 ‘347 323 
19* 402 ‘313 287 
20* 316 ‘308 282 
21 265 : 318: 2904: 
Height 600 km 1375 342 318 
1190 335 ‘311 
955 318 ‘204 
815 328 "304 
725 "340 -316 
615 “364 340 
508 ‘342 318 


n 


345 "364 "340 


398 


2 
2 
2 
2 
2 


398 
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‘419 


uw 
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Height 1100 km : 421 
°333 
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“23% 
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‘289 
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Empirical reduction of hydrogen line widths—The measures cover n =7 to 31, 
but not all on one exposure or at one level. There are not enough observations 
to determine the variation of Stark effect with height, so that we have taken all 
the measures together, assuming they refer to some average height between 
600 and 1200km. ‘The Stark effect has been taken to give a line width pro- 
portional to n(m—1). This is a good approximation for large m, and is amply 
accurate for our present needs. We shall find that for »<15 the Stark effect 
makes a very small contribution to the total line width; for n>25 it dominates 
the line profile. ‘The self-absorption broadening has been taken to vary only 
with the measured total line intensity, irrespective of height. What is left after 
Stark effect and self-absorption broadening have been removed is taken to be 
the Doppler profile, which is of course proportional to wave-length, although 
the range of variation with the lines measured here is less than ten per cent. 
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Fic. 2.—Measured profiles (single readings) of H16, H21 and H26 at height ~ 600 km. 1 mm=2'4 A. 


With these considerations in mind, empirical values for the Stark broadening 
and self-absorption broadening were estimated by graphical means and the 
results are in ‘Table IJ]. It was assumed that for this preliminary examination 
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TABLE IITA 
Doppler half-widths of hydrogen lines, empirical reduction 


100 km 600 km 1100 km 1200 km 
0°253A 
‘282 
308 
‘289 
289 
‘270 
‘253 
‘287 
268 
‘253 
294 
294 
282 


Tas_e IIIB 
Stark half-widths from empirical reduction 


n h.w. 

7 oo17A 
fe) 036 
15 085 
20 "154 
25 ‘243 
30 “352 


TasBie IIIc 
Corrections for self-absorption from empirical reduction 
(to be subtracted from the observed half-width) 

I Correction 

250 0-002 A 

500 ‘007 

750 ‘018 
1000 036 
1250 068 
1500 ‘121 
1750 +187 
1905 °233 
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Stark and Doppler widths could be combined as Gaussian errors. The self- 
absorption was estimated as a correction to be subtracted from the measured 
half-width. 

After removing Stark and self-absorption broadening the mean Doppler 
half-width is o-276A at a mean wave-length 3699 A. In obtaining this value 
a few of the strongest and weakest lines were rejected. This width gives a 
kinetic temperature 11 000 deg., if no allowance is made for turbulence. The 
correction for turbulence, following the earlier part of the paper, will depend 
on the assumed height. If we take this to be 850 km, with a &, (turbulence) =7, 
the corrected kinetic temperature is about 10000 deg. K. 

Theoretical calculation of hydrogen line widths.—This treatment assumes that 
for present purposes we can regard the chromosphere as a layer of gas with 
a given constant temperature and pressure at all points along any particular line 
of sight. The Stark profile of a hydrogen Balmer line has been taken from 
Verweij, as quoted by Miss van Dien.* ‘The profile is U(8), where 8 is the 
ratio of the field strength to a standard field F,=2-61N*%e. WN is the number 
of charged particles per cm® and e the charge on the electron. £ has a certain 
probability distribution, a property of the inter-particle fields in the gas. 

For convenience U(f) is represented by a Voigt profile, V(a,, AA,), where 
the parameters a,, AA, can be regarded as a fictitious damping width and Doppler 
width respectively. a, is expressed in units of A\,. By trial and error a 
satisfactory fit for U(8) was found with a,=o-4, the quality of the representation 
being shown in Table IV. Here w is the distance from the centre of the line, 
expressed as a fraction of the (fictitious) Doppler width; V(v) is the Voigt profile, 
taken from Hjerting’s computations (his Table 4).+ The ratio of 8 to wv has 
been chosen so that V(v) fits U(8) atv=1. The fifth column in Table IV gives 
the Voigt profile, and the sixth gives the Stark profile, each in units of its central 
value. 

The AA,,, =s,,.8F, of Miss van Dien is identified with AA, when v=1, B=1°-35. 
Hence AA,=1°35s,9, where Fy = 2-61 N®%e = 1-253 x 10°* x N*% e.s.u. and, with 
AA, measured in angstroms, 
n® 


-) {n(n —1)+2}. 


n2 — 


S, = 0°00256 ( 


Some typical values of AA, are given in Table V. 

We must now combine a real Doppler broadening with the Stark profile, 
the latter being now expressed as a Voigt profile with a pseudo-damping width a, 
and a pseudo-Doppler width A\,._ The real Doppler width, following Hjerting’s 
convention and in the customary notation is 


A 
a= /(=)2. 
My fe 
The resulting profile is another Voigt profile, V(a, AA), where AA, = 1/(AA?2 + AA,?). 
If we continue to express the (pseudo-) damping parameter in terms of the 
(pseudo + real) Doppler width, now AA,, it becomes a=0'47, . Some figures 


for the particular case 7'= 6000 deg., log N=11°'5, are given in Table VI. 


* Elsa van Dien, Ap. F., 109, 452, Table I, 1948. 
+ F. Hjerting, Ap. F., 88, 514, 1938. 
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TaB_e IV 

V(v) B U(6) V(v)/V(0) U(8)/U(e) 
0°67 0°00 0287 1'00 1°00 

57 0°68 243 0°85 0°85 

"36 1°35 "154 ‘54 54 

‘19 2°02 ‘082 28 ‘29 

“090 2°70 ‘O41 134 "143 

"049 3°38 ‘022 ‘073 ‘077 

0300 4°05 ‘O14 "045 "049 
0'021 4°72 0008 0'031 0'028 


TABLE V 
Values of AA, 
log N=12:0 II's II‘O 
5 o'0135A 0:00627A 0:00290 A 
10 0432 "0201 00931 
15 0951 0442 0205 
20 ‘169 0783 0363 
25 264 123 0569 
0°381 0177 00820 


TaBLe VI 
T=6000 deg. log N=11°5 

n AA, Ary; AA, A),/AA,q a 

5 0:00627A 0'1446A o145A 00432 0'O17 
10 ‘0201 1266 128 ‘157 063 
15 "0442 ‘1237 ‘131 338 135 
20 0783 ‘1227 146 537 ‘215 
25 123 1223 ‘173 oy 284 
30 0°177 0'1220 O'215 0823 0°329 


Having obtained the Stark plus Doppler profile for various n, the next step 
is to compute the effect of self-absorption. For this it is necessary to know the 
absorption coefficient «, at each point within the hydrogen line. ‘The first stage 
is to determine «,, the value of «, at the line centre, and this we have obtained 
by way of B,,hv, the total strength of the Balmer line. B,, is the Einstein 
probability of absorption, and may be computed from Menzel and Pekeris’ 
values of A,, by the relation 


g, 
By,= — 3 
£, 2hv' 


Ay. 
Table VII gives an outline of the results. 


TasLe VII 

n A log A log B log hy log Bhy 

5 4340 A 6-401 Q'512 12°660 2°172 

10 3798 4°850 8-389 12°718 3°107 

15 3712 3°962 7°823 12728 4°551 
20 3683 3°334 7°425 12°73 4°157 
25 3669 2849 7°139 12: 5872 
30 3662 2°452 6°897 12°734 5°631 

In order to find x, the values of B,,hv must be divided by a factor depending 


on the line width. First we remark that Hjerting, whose computations we use, 


36* 
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tabulated «,/«* + as a function of v, the distance from the line centre in units 
of the Doppler width, and normalized his profiles so that 


_* 
Jo Ko 2 


We now define an effective line width v,,, so that 


Ky “lk, /t 
wy age { aie M. 
Ko J v=0 Jo \Ko - 


(x,/K*),-9 is the actual absorption coefficient at the line centre, as tabulated 
by Hjerting. Then the relative values of the absorption coefficients at the 
Balmer line centres are B,,hv/2Av._, where 


a 
| Ky YT 


c c 

2AVeg = = 2A ge = y2 
a“ A 

A summary of the results of this section of the work is in Table VIII. 


~« 2Uer - AA,. 


TasBie VIII 
T =6000 deg. log N=11°'5 


B, oh 
AA, a (y/K5 v=o Verf log 2Averr log —2 “4 
~_ 2Averr 


o'145 / O'O17 0-981 0'903 10°619 13°553 
128 063 034 0948 10°701 14406 
‘131 "135 ‘866 1'024 10°765 15°786 
146 ‘215 797 I‘112 10°856 15°301 
‘173 284 748 1°184 10:960 16-912 
‘215 "320 714 1‘241 11076 16°555 


Suppose now that the optical depth of the chromosphere in the line of sight 
at the centre of a certain line is 7», and is 7)., at the point where the observed 
intensity is one-half the central intensity. ‘The intensity distribution within 
the line, including the effect of self-absorption, is taken to be 


I «(1-—exp[-—7]), 
so that \ = exp ( oc 70-5) 
2= Texp(—7) * 
The relation between log 7, and log 7,.;/7) is given in Table IX. 7, is proportional 
to By, hv/2Av., for the particular Balmer line which is in question. We 





TABLE IX 


log 7 log T9.5/To 
00 1699 

2000 1°699 

1-000 1°681 

1699 I 

0000 

0°301 

0°477 

0 '602 

0699 

1°000 

2°000 

3000 

4°000 

+ «%, is Hjerting’s «xo, the fictitious absorption coefficient at the centre of the line when the 
damping is vanishingly small. 
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proceed by assuming 7, for one convenient line of the series, and choose n= 30. 
Table IX then gives 7).;/7), and from Hjerting’s Table 4, with a=o-4.(AA,)/(AA,), 
values of v are chosen so that 


(=) (3) ia TH-5 
* * Hh SS 
Ko | KO J v=0 T? 


v is the distance from the line centre in units of the Doppler width, so that 
2u x AA, =2AA,,, should be the half-width of the line as observed. 


TABLE X 


Values of 2AXgai 
r 6000 deg. 6000 deg. 6booodeg. 6o000deg. 6o0o deg. 8000 deg. 
log N 12'0 12°0 II'5 II's II‘o II's 
To(n = 30) ov! 0°05 or! 0°05 ov! 0°05 





t= § 0'844A 0'722A 0682 A 0'624A 0:626A 0-702 A 
fe) 528 446 "412 346 358 388 
15 436 372 ‘310 272 262 "304 
20 478 446 "304 "290 246 314 
25 628 624 364 "352 256 378 
30 0872 0860 0°452 0°448 0'292 0466 


T 8000deg. r1oooodeg. 10000deg. 10000deg. 10000deg. 12000deg. 
log N II'5 II°5 I1°5 Ir‘O II‘o II‘o 
To(2 = 30) 0°03 0°05 0°02 0°05 0°02 0-02 


n= § 0'654A 0'776A 0°674A 0'734A o628A 0682A 
10 "346 ‘422 ‘344 ‘374 "318 ‘344 
15 ‘290 336 +306 "304 ‘284 308 
20 "314 "342 "340 "204 ‘204 "318 
25 ‘374 "400 "396 ‘312 ‘308 "332 
30 0°460 0°484 0°478 0°342 0°334 0°356 





Table X contains half-widths calculated in this manner for various com- 
binations of parameters within the ranges 


a] 


6000 deg. < 7 
II‘o 


12.000 deg., 


< 
< logN < 12°0, 
0°02 < T)(m=30)< 


orl. 
For a number of reasons, some of which will be briefly discussed in a moment, 
a perfect fit with the observations cannot be expected. Fig. 3 shows the 
observations with two computed curves: 
(a) T= 6000deg. log N=11°5 T) (w= 30) =O'1 
(d) 10 000 deg. 11'S 0°02 
Both curves lie fairly near the observations. Of the two, (6) appears to be 
somewhat the better, but it is doubtful whether much significance can be attached 
to this. The results suggest that an accurate determination of the kinetic 
temperature by this method is difficult, and could be made only from more 
comprehensive observations covering most of the Balmer series to m=30 at 
a number of different heights. 
The difficulties may be summarized as follows: 
(1) Owing to the large range of light intensity to be covered the present 
observational accuracy falls off at each end of the measured range of x. 
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(2) Measures at different heights 600-1200 km have been combined, although 
it is to be expected that both self-absorption broadening and Stark broadening 
will decrease with increasing height. The computed half-widths for the earlier 
members of the series are rather sensitive to the assumed “‘ Balmer optical depth” 
(expressed as 7, at m= 30), so that they are strongly dependent on height. There 
are two places at which this fact needs particularly to be borne in mind. In the 
measures quoted in Table II there is a strong suggestion that those from effective 
height 1200 km correspond to a substantially smaller 7) (n= 30) than the others, 
so that to combine them all together is not entirely satisfactory. We should also 
refer to the 1940 measures, in which the half-width for n=5 was 0-53 A. ‘This 
is smaller than one would expect either from extrapolation of the 1952 measures, 
or from the computed cases (a) or (b) above. Satisfactory data are not available 
for the height at which the 1940 measures were made, although 1500 km was 
tentatively suggested at the time. No great reliance can be placed on this figure, 
but it suggests a t, lower than for the 1952 measures, and qualitatively this 
agrees with the line widths. 


Half-width 
inA 








1 1 1 A 


lO 15 20 Z3 





Fic. 3.—Half-width of hydrogen Balmer lines plotted against n. 
—— Curve (a), for T=6 000 deg. Curve (b), for T=10 000 deg. 


(3) We have made no allowance for turbulence in our computations. If 
the increase with height suggested by Fig. 1 is substantially correct, turbulence 
is becoming important, even for hydrogen, at a height about 2500km. The 
turbulence makes a kinetic temperature determination from metallic lines virtually 
impossible by present techniques, and not very practicable from helium. 

(4) The computations have neglected all variations in the properties of the 
chromosphere along the line of sight. Even ina static atmosphere, what is observed 
is an integration over a wide range of physical parameters. If the chromosphere 
has to be treated as a field of filaments or jets the question becomes forbiddingly 
complicated. 
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Since this paper was completed, work by Kawaguchi has appeared* in which 
he fits theoretically computed hydrogen profiles to the 1940 measures made by 
one of us. His conclusion with regard to these lines is the same as ours, viz. that 
it is possible to interpret the observed profiles with a kinetic temperature of 
10000 deg. K or less, if self-absorption is taken into account. 
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Summary 


The apparatus function of the Oxford solar spectroscope has been deter- 
mined in the wave-length region A 5461-A 5945. Comparison with Plaskett’s 
determination of the core of the profile at A 4320 shows that Redman’s 
suggestion that the linear dimensions of the profile are independent of 
wave-length requires modification. Differences between the curves at the 
two regions can be accounted for by the existence of large-scale inhomo- 
geneities of the order of 3 x 10~° in the refractive index of the prisms. The 
wings of the apparatus function show a number of sharp subsidiary maxima, 
whose spacing is six to seven times the spacing of the subsidiary maxima 
of the theoretical diffraction pattern. They are attributed to diffraction by 
striae, known to be present in the prisms. 





Introduction.—In his well-known investigation of the central intensities of 
certain Fraunhofer lines (8,9) Redman made a very accurate and extensive 
determination of the apparatus function of the f/77-6 McClean prismatic spectro- 
scope. He found two noteworthy departures from the predictions of diffraction 
theory. Firstly the core was much wider and the wings stronger than the 
diffraction-theory curve. Secondly the linear dimensions of the core of the 
apparatus function over the spectral range A 4200-A 6400 did not vary linearly 
with the wave-length but were approximately constant—a slightly greater 
observed linear width at A 4200 being attributed by Redman to errors caused by 
temperature changes in the spectroscope. 

As a step towards finding the profiles of certain Na lines, the apparatus function 
of the Oxford prismatic spectroscope (7) has recently been determined in the 
region A5461-A5945. ‘The technique of this determination, essentially the 
same as that used by Redman and based on his experience, is briefly described 
in the final section of this paper. This apparatus-function determination, 
supplemented by the apparatus-function core determined at A 4320 by Plaskett 
(7), shows that some of Redman’s results cannot be applied to the Oxford spectro- 
scope without slight modification. The wings, moreover, show a number of 
subsidiary maxima which appear to be a newly observed phenomenon with 
prism instruments. ‘The new results and conclusions, together with possible 
theoretical explanations, are set forth in the first two sections. 

1. The core.—The continuous curve of Fig. 1 shows the core of the apparatus 
function as determined from the profiles of the emission lines A 5461, A 5770 and 
A 5791 given by a mercury-isotope lamp and from A 5945 given by a neon discharge 
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tube. Ordinates are intensities and abscissae are millimetres in the focal plane 
of the spectroscope. ‘The filled circles are the mean observed points and the 
vertical lines through them show the r.m.s. errors of these means found from the 
residuals. ‘The dashed curve is that found by Plaskett from the Kr line A 4320, 
the abscissae again being millimetres in the focal plane. Both the full line and 
dashed curves were obtained with field-lens illumination of the spectroscope slit. 
For comparison the dash-dot curve gives the profile of the mercury-isotope line 
A 5461 for a mode of slit illumination mimicking that given by the solar telescope. 











| 





0 +0.04 +008 +021 


Fic. 1.—Apparatus-function core. 
Abscissae: millimetres in the focal plane. 
Ordinates: intensity. 


The continuous curve refers to a mean wave-length of 45730, the dots being observed points 
and the vertical lines r.m.s. errors; the dashed curve refers to 44320. The dash-dot curve was 
obtained at 5461 with a mode of slit illumination mimicking that produced by the solar telescope. 


Comparison of the full-line curve for a mean wave-length of A 5730 and the 
dashed curve for a wave-length of \ 4320 shows fairly close agreement on the 
violet side. On the red side, however, there is a marked difference between the 
two profiles, that for A 4320 showing a strong satellite which is almost completely 
absent at 15730. ‘This difference is also well shown in Plate 14 where on the 
right appears a fifteen-minute exposure of 5461 and on the left a three-hour 
exposure of A 4320—the latter showing the satellite. 

A possible explanation of the appearance of the satellite at A 4320 can readily 
be offered. In accordance with Hilger’s standard practice the prisms of the 
spectroscope have been locally figured to compensate for inhomogeneity of 
refractive index in the glass. ‘This compensation, however, can only be effective 
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at the wave-length for which the prisms have been tested in Hilger’s Twyman- 
Green interferometer. At any other wave-length, the re-figured portion of the 
prisms will divert light from the central part of the line image. A simple numerical 
example shows that a very small variation in refractive index is sufficient to account 
for the satellite observed at 14320. Assume that we can replace the actual 
prism by two prisms placed one on top of the other. Let the upper one have 
refracting angle A and refractive index m and the lower an angle 4+5A and 
index n+6n, where for simplicity we shall assume 5n to be independent of 
wave-length. Then since at minimum deviation 
_ A+O/. A 
n=sin— /sin= ; 


the two prisms will produce at some wave-length, say A 5730, two lines separated 


by 
sp= 2 sin A/2 int { n cos A/2 -1}e4. (1) 
V 


V1—n? sin? A/2 ‘1 —n? sin? A/2 








To bring these two lines into coincidence at this wave-length we must make the 
right-hand side of (1) zero, so that for Chance glass 613369 of the Oxford prisms 
and A = 60° we must make 

5A = —1:23765n. (2) 
This corresponds to the case of a prism locally figured at wave-length A 5730. 
At A 4320 this composite prism will give two lines separated by 5@ according 
to (1). Inserting the appropriate refractive index for this wave-length for 
Chance glass 613369 and the value of 5A from (2), we find 

50 = —0:0700 6n. (3) 
Now for the equivalent of five prisms we find from Fig. 1 a satellite with an angular 
separation of o-1 mm/F, where F=8958 mm for the Oxford spectroscope. 
For a single prism the angular separation would be a fifth of this; hence from 
(3) the difference in refractive index between the top and bottom halves of the 
prism necessary to produce the observed satellite is only 
- : (: x sess) 32% 10° 

5 8958 





in = — =— 
0:0700 0:0700 
a very small amount. 

The strength of the satellite at 44320 suggests that the Oxford prisms are 
composed of a localized mixture in roughly comparable amounts of two glasses 
with refractive indices differing by 3:2 x 10-5. Equally the greater width of the 
violet wing at this wave-length suggests the presence of other irregularly dis- 
tributed inhomogeneities with differences of refractive index still smaller by a 
factor of the order of ten. It seems unlikely that large blocks of optical glass 
can at present be made without large-scale refractive index variations of the order 
of 10-° or 10-**, so we can conclude that in most cases local figuring will be 
necessary even to yield moderate performance in one wave-length region. For 
such locally figured prisms the width of the core of the apparatus function will 
be a minimum for that wave-length at which the re-figuring has been carried out, 

* F. Twyman (11) quotes W. M. Hampton as follows: ‘If a piece of glass is uniform in 
temperature while being held prior to the annealing operation, the only index changes introduced 
into it will be the plus or minus variations due to the temperature gradient during cooling, and 
these will probably be small, probably not more than a few units in the fifth decimal place....” 
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and for any other wave-length will depend on the difference between the refractive 
index at this wave-length and that at the wave-length of re-figuring, and also on 
the extent of the re-figuring. It therefore seems that Redman’s tentative sugges- 
tion that the width of the core is independent of wave-length should be modified 
in the light of the above facts. 

A word should be added about the dash-dot curve of Fig. 1. As already 
mentioned, this curve was found for a slit illumination mimicking that given by 
the solar telescope. In so far as this type of illumination is non-coherent and the 
field-lens illumination used for the other two curves is coherent, the greater 
width of the dash-dot curve is in accordance with van Cittert’s prediction (2). 
It should, however, be emphasized that this curve is based on long exposures 
during which atmospheric pressure changes occurred, and though a first-order 
correction has been applied for these changes (cf. Section 3) the resulting apparatus 





| ! | ' 











Fic. 2.—Apparatus-function wings. 


Abscissae: millimetres in the focal plane. 
Ordinates: logarithm of intensity. 
Continuous and dashed curves are red and violet wings respectively. 


function is subject to some uncertainty. The difference between the results 
obtained with the two modes of illumination may therefore be more apparent 
than real. 

2. The wings.—The red and violet wings of the apparatus function are shown 
in Fig. 2 as continuous and dashed curves respectively. Ordinates are proportional 
to the logarithm of the intensity, abscissae are millimetres in the focal plane, and 
the curves may be taken to refer to the same mean wave-length, A 5730, as the 
core. Because of the importance of the far wings of the apparatus function in 
fixing the observed central intensities of Fraunhofer lines, a fact originally pointed 
out by Allen and Redman (8, 10), the present determination has been carried 
down to g x 10° of the central intensity. 

In agreement with Redman’s results, the wings are considerably stronger 
than the maxima of the side bands of the theoretical diffraction pattern. They 
differ, however, from those of the McClean solar spectroscope in showing a number 
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of irregularly spaced subsidiary maxima. ‘These maxima are well shown in 
Plate 1a—an over-exposed spectrum of the Na D lines given by a Philips’ discharge 
lamp taken through a step wedge. Such maxima, however, tend to be washed out 
by atmospheric pressure changes, and their absence from the wings of theapparatus 
function of the McClean spectroscope may possibly be due to the occurrence of 
such pressure changes during the long exposures which Redman required to 
reach out into the wings (twenty-five hours as opposed to only four hours in the 
present work). 

The secondary maxima in Plate 1a bear some resemblance to diffraction 
maxima, and this resemblance is heightened by their rough symmetry about 
the core. Nevertheless, they cannot be due to diffraction over the whole aperture 
of the spectroscope, since this would lead to some eighty maxima on either side 
of the central peak instead of the twelve to fifteen observed. ‘The maxima may, 
however, arise from diffraction by striae in the glass of the prisms. In common 
with most large blocks of optical glass, the Oxford prisms contain a number of 
striae, well shown for example in Foucault tests. In so far as these striae are 
opaque, occur more or less at random throughout the prisms and, much more 
dubiously, may be regarded as having the same size, shape and orientation, 
they will give rise from Babinet’s theorem to diffraction side bands identical with 
those given by an aperture of the same size, shape and orientation as the striae. 
The irregularity in spacing and intensity of the secondary maxima in Fig. 2, 
somewhat resembling beats, is presumably a consequence of the striae not all 
having the same dimensions and orientation. If, however, the striae are such 
that they can be placed in one or other of several families, the members of each 
family being closely the same in the pertinent characteristics, then each family will 
give rise to its ownset of diffraction maxima and the superposition of these maxima 
will produce oscillations which wax and wane as do those in Fig. 2. Superposed 
on these hypothetical diffraction effects are of course effects due to inhomogeneity 
of the glass (for example, the satellite extending from 0-3 to o-8 mm on the red 
wing), together with a general background of scattered light, some considerable 
part of which is probably produced by the striae. In so far as this interpretation 
is correct and in so far as striae are common to all large blocks of optical glass, 
secondary maxima similar to those shown in Plate 1a may be expected to occur 
in the wings of the apparatus function of any large prismatic spectroscope. 

3. Experimental procedure.—In general the experimental procedure followed 
that of Redman, differing to some extent, however, in the method of slit illu- 
mination, sources of light, and photometric procedure. A field-lens was used 
to produce coherent illumination. The lens was illuminated by a diaphragmed 
light source and produced an extra-focal patch covering the full length of the slit, 
where a step wedge was situate:!, the image being in focus on the collimator and 
filling it with light. ‘The second arrangement was designed to mimic the non- 
coherent illumination produced by the solar telescope. This was done by means 
of a good-quality lens stopped down to the appropriate focal ratio and provided 
with a central stop to represent the obscuration caused by the Cassegrain mirror. 
In this case the image of the diaphragmed light source was formed on the slit. 
This second arrangement is uneconomical of light, and for reasons described 
later produced results of very low weight. For both arrangements the slit 
width was set at 0:036mm, the value used for solar work in the region 
A §500-A 6000. 
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The three light sources used were a water-cooled '*Hg isotope lamp made 
by G.E.C., a B.T.-H. neon discharge lamp, and a Philips’ sodium lamp. The 
isotope lamp has been described by Barrell (1), who states that at 17 deg. C the half- 
width is 0-008 A, being chiefly due to, Doppler effect and collisional broadening, 
a value less than one-tenth of the apparatus-function half-width at the wave- 
lengths concerned. ‘The lines A 5461, 45770 and 45791 were used. The neon 
discharge lamp contains neon at a pressure of a few millimetres in a tube 9g mm in 
diameter; the operating current is 0-5 amp. Jackson (4) finds that the con- 
ditions of excitation necessary to obtain the lines free from self-reversal are not 
particularly stringent. From Jackson’s published interferograms the _half- 
width of the particular line used, A 5945, is <o-o2 A, which is only one-sixth the 
apparatus-function half-width at this wave-length. ‘The weak satellite 0-02 A 
to the violet of the main line, due to **Ne, present in atmospheric neon to the 
extent of one part in ten, can be neglected. Redman has also used neon lines 
(8). Finally, a Philips’ sodium lamp was used for the far wings. ‘The use of 
such a bright source enabled an intensity of g x 10-5 of the central intensity to be 
recorded in exposure times extending only up to four hours. Owing to the 
danger of self-reversal, this lamp was not used for the core of the profile. Table I 
gives a summary of the observational material obtained with the three light 
sources. 

TABLE I 


Observational material 

Minimum 
intensity 
recorded 
A 5461 coherent 6 15™—6h 0°02 

A 5770, A 5791 ‘a 3 2h—4}h ‘4 

A 5890, A 5896 - 308 -4h -00009 
A 5945 Ra 6 1om—4h 02 

A 5461 non-coherent 5 2h—gh > 


Mode of No.of Exposure 


Wave-length illumination _ plates times 


Plate calibration was performed by means of a step wedge consisting of six 
steps of rhodium evaporated on to glass together with a ‘“‘clear glass” step. The 
wedge, which was made by A. Hilger, Ltd according to Plaskett’s specifications, 
is designed for use on the slit of the solar spectroscope, a method of use which 
automatically corrects for scattered light. In order to avoid the formation of 
Edser-Butler fringes, which would occur if the glass base upon which the 
rhodium is deposited were plane parallel, the base is made in the form of a double 
prism, the two glasses being chosen so that the deviation of the light passing 
through the wedge is nearly zero for all wave-lengths in the range for which the 
step wedge is required. (Owing to the large focal length of the spectroscope, 
a single prism cannot be used, since the corresponding deflection of the beam would 
throw the light off the collimator.) ‘The wedge is so constructed that by shifting 
it parallel to itself across the slit the light can be made to pass through the glass 
base without passing through the rhodium steps. This makes it possible to 
correct the characteristic curves for any residual deviation produced by the glass 
base, and more important, for any non-uniformity of the slit illumination. 

The photometric procedure differed from that of Redman (g), who used 
either a separate calibrating spectroscope or a tube sensitometer. The method 
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adopted here was the method of self-calibration. Inthis method a single exposure, 
taken through the step wedge, exhibits both the profile under measurement and 
the calibration marks necessary for constructing the characteristic curve, the 
profile being obtained through the “clear glass” step. An additional exposure, 
with the wedge displaced to one side, allows one to correct the characteristic 
curve for any non-uniformity of slit illumination in the manner described above. 
The great advantage of the method is that both profile and calibration marks 
are obtained under identical photometric conditions. Provided the centre 
of the wedge exposure is not ‘burned out’’, the characteristic curve can be 
constructed in terms of the central intensity. The reduction then gives directly 
the intensity at any given point in the profile in units of the central intensity. 
In determining the outer parts of the profile, however, it is necessary to over- 
expose the centre of the line (cf. Plate 1a), and the characteristic curve must 
then be constructed from some reference point other than the centre. If the 
intensity at this point is known from the reduction of previous plates in the series, 
however, relative values may readily be converted to central intensity units. 
The subsidiary maxima described in Section 2 were chosen as reference points ; 
they served very conveniently for the measurement of distance from the centre 
of the profile as well as intensity. 

Backed Ilford Rapid Process Panchromatic plates were used throughout. 
The plates were ‘‘ brush” developed, using a felt pad, in order to reduce Eberhard 
effect, and were developed in Parkhurst two-solution developer, which produces 
extremely little background fog. Microphotometry was performed with the 
Observatory microphotometer (6), a Moll-type instrument, using a plate-tracing 
magnification of fifty. Since the width of the microphotometer-slit image on 
the plate is only about one-seventh the half-width of the profile being measured, 
the microphotometer reproduces very closely the true profile. 

The Oxford solar spectroscope is immune to temperature changes in its 
surroundings owing to thermostatic control (7). When the source is so faint 
that the core of the profile requires a long exposure time, however, there is a 
chance that a sufficiently large atmospheric pressure change (accompanied by 
a change of the refractive index of the air in the spectroscope) may occur to distort 
the profile seriously. ‘This condition is not met with in solar use, of course, 
nor was any difficulty experienced in any of the twenty-five plates taken with 
field-lens illumination, since the exposure times for the core of the profile were 
short in the case of A 5461 and 5945, while in the case of A5770 and A 5791, 
although the exposure times were relatively large (see Table I), the pressure 
changes happened to be moderate, never exceeding o-2mm. However, three 
of the five plates taken with non-coherent illumination, none of which recorded 
more than the core of the profile, were affected by pressure changes of 0-9, 0-9 
and 1:8 mm of mercury respectively. Calculation shows that the latter value 
causes a line shift on the plate comparable with the half-width of the profile. 
Rather than discard the three plates, which gave profiles whose half-widths 
were well correlated with the pressure changes, it was decided to correct the 
profiles for the known changes which had occurred. ‘This was done as follows: 
if we assume a linear variation of pressure with time which causes a total line 
shift 5A, then the observed profile is the true profile obliterated by a rectangular 
obliterating function of width 5A, and the true profile may be found from the 
observed profile by solving the equation of obliteration. In view of the magnitude 
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of the shifts, it is clear that low weight should be attached to the curve obtained 
with non-coherent illumination. 

Finally a few remarks are necessary concerning the step-wedge calibration. 
Calibration was originally performed using the solar spectroscope and four 
“‘neutral”’ filters which had been calibrated at six wave-lengths by the National 
Physical Laboratory. This calibration, however, gave a suspicious wave-length 
variation for the transmission of rhodium. A check calibration using a Zeiss 
platinum step wedge (3) showed that the Hilger step wedge was in fact neutral, 
so it was decided to carry out a full-scale re-calibration using a Hilger E.316 
spectrograph. ‘The photometric standards for this calibration were seven per- 
forated diaphragms (5). A continuous source of light was used, microphotometry 
of the plates being performed at the three wave-lengths at which the original 
calibration had been made. The results of this final calibration (see Table IT) 
show that over the wave-length range involved, and within the experimental 
error, the transmission of rhodium is independent of wave-length. 


TABLE II 
Transmission of Hilger rhodium-on-glass step wedge 


A 5461 A 5780 A 6143 Mean 
1000 1-000 

‘006 668 
003. “427 
259 
155 
‘102 


I 


: 


652+ 003 = “649 
‘431+ 007 = “433 
266 + -004 265 


"005 ‘656 
"430 
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003 
003 


“002 


+ I+ H+ I+ 


‘158+ ‘oor ‘155 156 
‘108+ ‘oor 103 


‘0737 


‘104 
"0004 ‘07315 * ‘0735 


‘OO! 


He Ie I+ + I I+ 


> 
+ 
+s 
7» 4 
+s 
+ 


+ I+ 


0003 


It was suspected that the error in the original wedge calibration was due to 
unexpectedly sharp variations in the transmissions of the four filters, making it 
unsafe to interpolate values from a smooth curve through the six points at which 
the filters had been calibrated. One of the filters was accordingly returned to 
the National Physical Laboratory for re-calibration at intermediate points. 
The resulting transmission curve confirmed this suspicion. Moreover, when the 
step-wedge values given in Table II were applied to the calculation of the trans- 
mission of the filter, using the original wedge calibration plates, the filter trans- 
missions found by the National Physical Laboratory were correctly reproduced. 
Confidence can therefore be placed in the values given in Table II. 


The writer wishes to express his deep sense of gratitude to Professor H. H. 
Plaskett for his unfailing encouragement and generous assistance, particularly 
in the theoretical interpretation of the results, and to thank the Master and Fellows 
of Balliol College and the Department of Scientific and Industrial Research for 
awards during the tenure of which the work was carried out. 
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FIELD REVERSAL IN MAGNETIC VARIABLE STARS 
P. A. Sweet 


(Communicated by the Director of the University of London Observatory) 


(Received 1954 August 24) 


Summary 


In connection with oscillation theories of magnetic variable stars the 
limitation in the variation of the observed magnetic field strength H, imposed 
by Cowling’s boundary condition is examined. 

Discussion is restricted to fields that remain symmetrical about a fixed 
axis during their variation. In spherical stars with a coefficient of limb 
darkening 0°45 it is found that, provided suitable gas motions can arise, 
reversal is possible only in the following cases :— 

(i) dipole-type fields with distortions that are not symmetrical about 
the equatorial plane, when the inclination i of the magnetic axis to the line 
of sight exceeds 79°; 

(ii) symmetrical quadrupole-type fields with symmetrical distortions, 
when 50°<i< 69°; 

(iii) quadrupole-type fields with unsymmetrical distortions, for all values 
of 7; 

(iv) certain fields of higher complexity even when symmetrical, for all 
values of 7. 

Statistical considerations indicate that it is improbable that dipole or 
symmetrical quadrupole-type fields can account for observations. 

The paper is not concerned with the dynamics of the gas motions. 





1. Introduction.—Cowling has shown (1) that in the so-called oscillation 
theories of magnetic variable stars, where the magnetic field variation is ascribed 
to distortion in the lines of force by moving material in the star, the field in the 
reversing layer would be equal to the surface value of the field outside the star. 
As this latter is a potential field a condition is imposed on the reversing layer field. 
This condition limits the possible variation in the resultant field strength H, 
derived from the longitudinal Zeeman effect in the star’s spectrum. In particular 
the question arises as to whether this resultant field could change sign as is observed 
in a considerable proportion of magnetic variables. By showing that no reversal 
occurred in a particular case of extreme distortion, Cowling illustrated the difficulty 
of finding a distortion which could produce a reversal. It is the purpose of the 
present paper to examine this question more fully. 

2. The invariance of field type.—As mentioned in the Summary the discussion 
will be restricted to axially symmetric fields. By this it is meant that the fields 
remain axially symmetric throughout the oscillation cycle. It is shown in 
Section 7 that it is the type of field, whether dipole or quadrupole, etc., that 
determines the possibility of reversal when the star is viewed from any particular 
direction in relation to the magnetic axis. 

Since the field outside the star is axially symmetric and potential it has no 
¢é-component (component perpendicular to meridian planes). Furthermore, 
because the movement of gas inside the star is supposed to be axially symmetric, 


37 
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the meridian-plane component of the field inside the star is independent of the 
¢-component.* Hence the meridian-plane component is all that needs to be 
considered in the present work. It will therefore be regarded, in what follows, as 
the whole field. 

The field is determined by the system of lines of force in a meridian half-plane 
bounded by the axis. Any field produced by finite electric currents in the star 
possesses neutral points in this plane; these are of two types. A point encircled 
by closed lines of force in its neighbourhood will be called a circular neutral 
pointt, while a point at which two lines of force intersect will be called an x-type 
neutral point. 

A dipole-type field is defined as any field, vanishing at infinity, possessing one 
circular neutral point. ‘There can be no x-type neutral points in such a field. 

A quadrupole-type field is defined as any field vanishing at infinity, possessing 
two circular neutral points of opposite polarities, while the x-type neutral points 
involved lie on the axis. In an unsymmetrical quadrupole-type field there are two 
of these points, the positions of which depend on the positions of the circular 
neutral points and the strengths of the electric current systems associated with 
them. If the field is symmetrical about the equatorial plane one of the points 
is at the centre of the star and the other is at infinity. 

In the neighbourhood of circular neutral points an electric current density is 
present, consequently all these points lie inside the star. ‘The typical times of 
variation in magnetic variable stars are sufficiently small, and the electrical con- 
ductivity sufficiently high, to use the familiar picture of lines of force frozen into 
the gas. Since the field referred to is an independent meridian-plane field all its 
lines of force, except a finite number of limiting ones that extend to infinity, are 
closed curves. Because gas cannot cross the lines of force any given closed tube 
of force contains a constant total mass of gas. The requirement of mass con- 
servation therefore shows that no tube of force can be created or destroyed. 
Consequently the number and polarity of the circular neutral points are unaltered 
by any gas motion in thestar. The field type is therefore invariant. 

3. The reversing layer field—The external field is potential, solenoidal and 
axially symmetric. It can therefore be expressed in the form 


H = — grad V, (3-1) 
where 
bead R n+1 
V=RE a,(=)" Py(u) (3.2) 
n=1 


The a,, are constants, the P,, are Legendre polynomials, R is the radius of the star, r 
is the distance from the centre of the star and px = cos 8, @ being the angular distance 
from the axis. Cowling’s boundary condition is that the field H, in the reversing 
layer is the surface value of the above field, hence H, is given by 


H,,= > (n+1)a, Py(H) (3-3) 
H,=sin@ ¥ a, P,'(). (3-4) 


n=1 





* This would not necessarily be true if the Hall effect had to be allowed for; this effect, however, 
is negligible in stars. 
+ This term also includes neutral areas encircled by lines of force. 
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The surface of the star has been assumed spherical; according to Cowling (2) the 
effect of ellipticity is not large, consequently the assumption is not considered to be 
unduly restrictive. 

4. The observed field strength H,.—Denoting the line-of-sight component of H, 
by H,, and the surface brightness by J at any point on the star’s disk, H, is given by 


H,= | | i Hllds / | ie 14S, (4.1) 


where dS is an element of area of the disk. Referring to Fig. 1 it is seen that at a 
point O on the surface of the star 


H,=H,,cos 6, — Hg cosy sin 6,. (4.2) 
From the theory of limb-darkening J is usually expressed in the form 
I=I,(1-—u+ucos6,) (4.3) 


when u is the coefficient of limb-darkening. An element of area of the surface of 
the star is sin 0, d0, dd,, thus an element dS of the disk is given by 


dS =cos 6, sin 6, d6, dd). (4.4) 


= magnetic pole 


on™ 


To Earth 





~ 
"me 


Fic. 1.—Geometry of star with magnetic axis inclined to line of sight. 
Using the last three expressions, together with equations (3.3) and (3.4), equation 
(4.1) can be written 
6 o me [ 


H,= La | 
: 7(3-U) nai : $,=0- u,=0 


[(n+1) P,(u) 4, — P,,'(u) sin @ cos 7 sin 4,] 
x [I — ut up|, dy, dd, (4.5) 
where ».,=cos@,. From the spherical triangle L PO 
sin 6 cos 7 sin 6, = fy — fy, (4.6) 
=o, + Sinzsin 6, cos dy, (4.7) 
where p,=cosi. Differentiating equation (4.7) partially with respect to p,, 
keeping ¢, constant, 


SiNZCOSh, — fg — HH 
sind, — sin?6, ° 





i =Ho~ Pa (4.8) 


sy" 
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Thus 
0p ’ eaee 
P,,'(#) sin 8 cosy sin 8, = P,,'() (Ho — wpa) = sin? A, in, [P,(#)]- (4-9) 


Using this relation the second term in the integral in equation (4.5) can be inte- 
grated by parts, and (4.5) further reduces to 


[ f. i P,(u) a, | [(1 —u)({n — 2} 42+ 1) 


+u({m—3}uy9+ 2p)]dey. (4-10) 
By integrating the addition formula (3) for Legendre polynomials with respect to 
¢, it is seen that 


ao 


= —— La | 
© WF—Bans  Jant 





J), Pale) de =m Py (Hs) Palio) (4-11) 
Using this formula in conjunction with the following relations 
[Pato du= 2), J. wPal)au= as; 
(n—2) [ut Py)de= — Gari, + (4.12) 
(n—3) [ us P,(u)du= — Preericar 





forn=1,2,... 0, where P_, (0) is defined by P_,(0)=1, H, can finally be written 
in the form 
2u(n +1) 


6 2 I—u 
= 2 E ay Palte)| oe Pr) eres Pree) |. (4-13) 


5. H, in terms of H,(y).—By transforming equation (3.3) the a, are given by 








+41 
ay = M8 Hu) Pal) de (5.1) 

Equation (4.13) can then be written 
HAuo)=| Half (4 #0) de, (5.2) 

where 
fata) = = (n+¥) 
I—u 2u 

x Fences P,,s(0) + dataiate P »-2(0) |P nH) Pr(Ho)- (5-3) 


Forlargen, P,,(0)—=cos = He ; also, | P,,(u) P,,(u@o)| <1 for u and py in the ranges 
nv 


|| <1, ]9] <1. Thus the terms in the above series tend to zero at least as fast 
asn-*, The series is therefore convergent, although rather slowly. f(j, 49) was 
computed numerically, with u=0-45, for a series of values of « and py. The 
function is shown graphically in Fig. 2a. For the purpose of calculation the first 
32 terms of the series in equation (5.3) sufficed, in general, for the order of accuracy 
required. ‘Tallquist’s tables (4) were used for the P,. In the region »=1, 
© <jg <0°2 it was found necessary to make a correction for the remainder of the 
series, ‘This was carried out by taking the asymptotic forms of P,,(o) and P,,() for 





‘asofins 4vjs yo xnyf JNauSou fo sida} ut fF] Sutats suoyguny—'e *D3.J 














eo= 


s9o0= 


DO Stee eee 





9O= 


sen" 
O=°w 


syuiod buiusny 
JNOUYIM PUD UJIM SaAuNd wineie tiles aia 
buryouodas sansn> buyjiwy --- Saree 


PUD UjIM SBAIND SajDIDGas “$1 OseW 





























554 P. A. Sweet Vol. 114 


large m and making the usual approximate transformation into integrals. The 
functions f,=f(u, 49) +f(—, Mo) and f,=f(u, 49) —f(—, wy) were also derived ; 
these are shown graphically in Fig. 24, c. 

6. The contribution of a single shell of force to H,.—Consider the magnetic field 
ina shell of revolution about the magnetic axis, defined by two neighbouring lines of 
force in a meridian-plane cross-section. Let the magnetic flux across a strip of a 
surface of revolution that intersects the shell be 2~AQ. AO will be called the 
strength of the shell; it will be reckoned as positive if the field is directed in the 
clockwise sense about the direction of increasing ¢, and negative when in the 
opposite sense. A point of exit is where, travelling in the clockwise sense along 
the line of force, the line leaves the star; likewise, a point of entry is where the line 
enters the star. At any point of intersection of the shell with the surface of the star 
let « and »+Ap, where Au is positive, define the intersections made by the two 
neighbouring lines of force. ‘Then 


+AQ=H,,(u) Ap, (6.1) 
the positive sign being taken for a point of exit and the negative sign for a point of 
entry. Equation (5.2) therefore shows that the contribution to H, made by the 
shell at this point is + AO f(y, 19), the positive and negative signs again being taken 


for points of exit and entry respectively. ‘The total contribution of a shell is 
therefore 


AQ LS (+15 40) —S (#2 Ho) ], (6. 
where =p, and ;, at the points of exit and entry, respectively. 
If f(4, #9) is monotonic in yu for any particular jy the expression (6.2) cannot 


change sign for any change in the positions of the points of intersection. ‘The 
argument holds for any number of re-entries of the shell into the star. ‘Thus for 
such a value of 1, no distortion of the shell can reverse the sign of its contribution to 
i. 

A change in sign of the contribution is only possible, given suitable distortions, 
when py has values for which f(y, 9), as a function of «, has at least one turning 
point. 

For a given , the maximum and minimum contributions that could be made by 
theshell, with AQ takenas positive for definiteness, are when f (1.1, 49) —f (42, My) has 
its maximum and minimum values, respectively, subject to 4, >j9. ‘To take an 
example, when the star is viewed pole-on (1) = 1) the graphs in Fig. 2a show that the 
maximum contribution is 0-77AQ, occurring when the points of exit and entry are 
convected to the two poles of the star. ‘The minimum contribution in this case is 
positive, but can be made arbitrarily small by convecting the points of exit and 
entry sufficiently close together. If the shell can be convected entirely below the 
surface its contribution of course vanishes. ‘The contribution therefore cannot 
change sign. ‘To take a second example, when the star is viewed in the equatorial 
plane (jg =0) the maximum contribution is o-10AQ, occurring when the point of 
exit is at the equator and the point of entry is at the south pole (uj=—1). The 
minimum contribution is — 0-10AQ, occurring when the point of exit is at the north 
pole (9 = 1) and the point of entry is at the equator. ‘The contribution can there- 
fore change sign. 

7. The behaviour of complete fields\—(a) Unrestricted dipole-type fields. The 
definition given in Section 2 shows that AQ has the same sign for all shells. ‘The 
results of Section 6, in conjunction with the graphs in Fig. 2a, therefore show that 
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H, cannot change sign if 0-19 <j») <1, but that reversal is possible, for suitable 
distortions, ifo <py<o-19. Interms of the inclination 7 of the magnetic axis to the 
line of sight, reversal is therefore impossible if i<79°. It could take place, how- 
ever, for inclinations greater than 79° if movements of gas of sufficient amplitude, 
alternately northwards and southwards, took place at the surface of the star. 

The maximum range in H, for cases in which reversal is possible arises when 
the star is viewed in the equatorial plane (u,=0). H, can then vary between 
+ 0-10Q, where Q is the total shell strength of the field when these extreme values 
occur. ‘The maximum value +0-10Q is effected when all the points of exit are 
convected into a small area near the equator while the points of entry are com- 
pressed into a small area near the south pole. The minimum value —o0:10Q0 
arises when the points of exit are all near the north pole and the points of entry are 
near the equator. ‘The shell strength Q, for a dipole-type field is given by inte- 
grating equation (6.1) thus: 

“1 
Qo=| Hep(u) dp, (7-1) 
P 
where y: =p’ is the point at which the limiting shell of force that just fails to leave 
the star touches the surface. For a simple dipole field, for example, 


H,, 1] “H,, 


pe’ -_ O, (7.2) 
where H, is the polar field strength; thus 


Qo= 3H, (7.3) 
The extreme values of H, attainable by distorting this field are therefore 
+005H,QQ». Ifthe movement of gas producing the distortion is confined toa 
thin layer near the surface of the star, such as might be expected from Cowling’s 
work (1) on horizontal oscillations, then relatively few new shells of force would 
be convected to the surface from inside the star, and few existing ones would be 
carried below the surface. According to a theorem of Bondi and Gold (§) no lines 
of force running entirely below the surface can appear above the surface whatever 
the motions. ‘This theorem is based on the fact that in hydrodynamic flow no two 
adjacent particles can be separated. However, hydrodynamic flow often breaks 
down at surfaces, as for example when a greased body in water floats to the top and 
breaks surface. Fluid motion is free along lines of force, hence the surface of shells 
of force behave like greased surfaces in this respect. The Bondi and Gold 
theorem therefore breaks down; the effective shell strength QO might therefore 
vary if suitable motions arose in the interior of the star. For example, if the neutral 
point were convected towards the surface it might be expected that Q would 
increase, and vice versa. But, as has just been mentioned, for surface motions 
O=Q, in any case, thus H, would vary between + 0:05 H,. 

Hence in general the variation in H, is within +0-05 H,. ‘To take an example, 
for the star HD 125248 Babcock (6) gives a variation of approximately + 2000 gauss 
in H,.* Therefore in the undistorted field H, >4 x 10 gauss. Moreover, since 
the case on which this limit is based probably represents a far more severe distortion 
than could be realized in any actual gas motion, the limit would probably have to be 
exceeded by a large margin in order to show the required range in H,. 


* The figures — 6000 to + 7000 actually quoted refer to an equivalent polar field strength equal 
to H,/0°303. 
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(6) Symmetrical dipole-type fields. It should be noted in the previous 
example that the gas motions required to reverse H, are not symmetrical about the 
equatorial plane. If the actual motions are symmetrical and the undistorted 
field is symmetrical, as seems more likely for mechanical reasons, the field remains 
symmetrical during the oscillation cycle. In this case H,, (—)= — H,,(u), and 
equation (5.2) can be written 


H&ivo)= | Haste) fas to) de (7-4) 


The contribution of a single shell of force intersecting the surface at +p is 
therefore AOf,(u,49). The graphs in Fig. 2b show that f,(, 49) 20 for all 
relevant . and py. Reversal of H, is therefore impossible in all cases. 

(c) Symmetrical quadrupole-type fields. Here H,,(—y)=H,,(u), thus 
equation (5.2) can be written 


1 

HH) = | Her (H) fa (Hs H0) de (7-5) 
Again only the field in the northern hemisphere (o<y <1) need be considered. 
This part of the field contains only one circular neutral point, and it can therefore 
be treated as an unrestricted dipole-type field. (The same argumentas in part (a) of 
this section can therefore be applied in conjunction with the graphs of f, (u, 449) in 
Fig. 2c. ‘These graphs show that reversal is allowed if 0-36<p,<0-65, corre- 
sponding to 49°-5 <i<68°-9. 

Using arguments similar to those employed in the second part of (a) it can be 
shown that for HD 125248 H,, >8 x 10° gauss, where H, is the polar strength of an 
undistorted quadrupole field. 

(d) Unrestricted quadrupole-type fields. If 0-19 <p <1, f(u, Mo) iS mono- 
tonic in x, thus the two systems of shells associated with the two circular neutral 
points make contributions to H, of opposite signs. By convecting all the lines of 
force of one of the systems, say the one with positive shell strength, assumed to be 
in the northern hemisphere for definiteness, into the same region, or by convecting 
them below the surface, the contribution of this system can be made arbitrarily 
small. ‘The negative system would then predominate and H, would be negative. 
By similarly suppressing the negative system instead H, could be made positive. 
H, can therefore be made to reverse. 

If 49 <0-19, then by sufficiently suppressing one of the systems the field can be 
made to behave like a dipole-type, and the results of part (a) of this section show 
that reversal is again possible. 

(e) Other fields. An example will be given of a symmetrical field that can 
change sign for all inclinations. 

Consider a symmetrical field with four circular neutral points, the two neutral 
points in each hemisphere having opposite polarities. ‘The contributions of the 
two systems in one hemisphere have opposite signs, hence by suppressing each 
system in turn H, can be made to change sign. Surface motions would be 
sufficient for producing this reversal. 

8. Statistical considerations.—There are two statistical reasons why it is 
unlikely that any axially symmetric dipole-type or quadrupole-type field with the 
additional restriction of symmetry about the equatorial plane can account for the 
fields actually observed. 
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In the first place, the ranges of i within which reversal is possible are rather 
narrow. If it is assumed that the axes of the stars are randomly orientated then, 
for unsymmetrical dipole-type fields the condition i<79° means that less than 
19 per cent of the stars could have reversing fields. For symmetrical quadrupole- 
type fields the range 49°:5 <i <68°-g means that less than 29 per cent of the stars 
could show reversal. Further, for the reasons mentioned in (a) of Section 7, the 
above limiting percentages are probably far in excess of what would actually occur. 

A recent table of magnetic variables given by Babcock (2), however, lists 35 
stars of which 37 per cent are definitely reversing, while it seems likely that some 
of the others will eventually turn out to be reversing. 

Secondly, it was shown in Section 7 that, for the unsymmetrical dipole types, 
the median strength of the field at the poles of a star in which reversal is observable 
is at least 20 times the extreme values observed in H,. The same star if viewed 
pole-on would therefore show an H, of at least 6 times the extreme reversing value. 
A similar argument holds for the symmetrical quadrupole types. It would be 
expected from this that, on the whole, those stars that do not show reversal would 
have a much higher H, than the reversing ones. Babcock’s table, on the other 
hand, lists a total of 13 reversing stars for which the mean range in 
H, ( =0°303 x tabulated H,) is about + 840 gauss, while 10 non-reversing stars in 
which the variation is given have a mean H, of about 810 gauss, a comparable 
value. 

g. Conclusions.—It has been shown in the present paper that, in oscillation 
theories, the simplest magnetic fields, i.e. axially symmetrical dipole and quadru- 
pole types are capable of showing a reversing H, when suitable gas motions can be 
maintained in the stars, provided the magnetic axes are suitably inclined to the line 
of sight. ‘The gas motions referred to need only be confined to the surface layers. 
It is unlikely, however, for statistical reasons, that, with the exception of the non- 
symmetrical quadrupole types, such fields could account for the observations. 
More complex or non-axially symmetric fields are required in oscillation theories. 

University of London Observatory, 

Mill Hill Park, 


London, N.W.7: 
1954 August 16. 
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ON THE DEVELOPMENT OF THE Fe I CHROMOSPHERIC LINES 
DURING THE 1951 ECLIPSE OF 31 CYGNI* 


Anne B. Underhill 
(Received 1954 May 3) 


Summary 

The development of the Fe1 lines AA 3920, 3922, 3927 and 3930 in the 
chromospheric spectrum of 31 Cygni during ingress is studied from high- 
dispersion spectrograms. Intensity tracings of the B-type spectrum and 
chromospheric features are presented and a light curve of the B star during 
ingress is derived. It is found that the light of the B star at AA 3915-3950 is 
reduced relative to that from the K star for several days prior to first contact. 
Evidence is presented that this dimming, which is probably caused by 
continuous absorption, and a sudden increase in the equivalent widths of the 
Fe 1 lines occur when the line of sight begins to pass through a dense region 
in the atmosphere of the K star. It is found that the Fe 1 lines are somewhat 
broadened by turbulence. 





Although chromospheric absorption lines of Ca 11 appear in the spectrum of 
31 Cygni some two to three months before the total eclipse of the B star by the 
K star, chromospheric lines due to Fei are not observed until just a few days 
before eclipse. ‘The present paper traces the development of some of these lines 
during the 1951 eclipse of 31 Cygni. A good series of high-dispersion plates was 
obtained at the Dominion Astrophysical Observatory during the last days of 
ingress, but because of unfavourable weather few plates were obtained during the 
first days of egress. Consequently it is possible to study in detail the development 
of the Fe 1 chromospheric lines during ingress only. All of the plates were obtained 
under the direction of A. McKellar and all are calibrated for intensity work. 

The observations—The chromospheric Fe 1 lines are barely perceptible on a 
plate taken 12 days before totality and are definitely absent 15 days after totality. 
The present study will be based on the series of third-order grating plates, 
dispersion 4-6 A mm at A 3925, obtained during ingress, together with plates of 
the K-type spectrum obtained during totality and plates of the composite 
spectrum obtained in June and July of 1952. All of the plates used give spectra 
of excellent definition and good density in the region 3900<A< 4000 A. 

In order to separate the chromospheric features from those due to the B-type 
star, it is convenient to study chromospheric lines which fall in regions free from 
B-type absorption lines. For this reason and because the lines appearing early in 
the chromospheric spectrum are intrinsically strong lines arising from the ground 
level of the carrier atom, the lines AA 3920-260, 3922-914, 3927°922 and 3930-299 of 
multiplet number 4 (1) of Fe 1 were chosen for study. The other lines in this 
multiplet were found to be unsuitable for detailed study either because they lay in 
the wing of a B-type hydrogen line, or because they were disturbed by strong 
features which could not be unambiguously defined in the K-type spectrum. 


* Contributions from the Dominion Astrophysical Observatory, No. 37. Published by permission 
of the Deputy Minister, Department of Mines and Technical Surveys, Ottawa, Canada. 
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Microphotometer tracings magnified 100 x in the direction of the dispersion 
were made of all the plates for the wave-length region 3700 <A <4100A, and these 
density records were transformed to intensity tracings magnified 200 x in the 
direction of the dispersion and rectified so that the apparent continuum lay at a 
height of 10 inches above the level of zero light for all wave-lengths. A correction 
of 5 per cent of the continuous spectrum was made to allow for the grating ghosts. 
In all the spectra peak intensities were found to occur at certain well-defined wave- 
lengths. ‘The apparent continuum was determined by putting a smooth line 
through these points of peak intensity. Particularly in the region 3900 <A <4000A 
it was found that this apparent continuum could be consistently and uniquely 
defined. All the measures of intensity on any one plate are made with reference to 
this apparent continuum, which hereafter will be referred to as “‘ the continuum” 
In the region 3915 <A <3950A, which is to be studied here, the continuum repre- 
sents the mean peak intensity of short spectral regions near A 3915 and A 3950. 
At all phases in the light variation of 31 Cygni, the monochromatic flux at these 
wave-lengths is greater than that at any neighbouring wave-lengths. It is not 
assumed that this mean peak intensity, called here the continuum, represents the 
true intensity of the light which would be emergent from both stars if there were no 
absorption lines present. ‘That level of light intensity represents strictly the sum 
of the continua of the two stars. However, the true continuum is an unobservable 
quantity for a K star in the wave-length region being studied and is of limited 
practical interest. 

In the study of narrow spectral lines, the distortion of the true shape of the 
absorption line by the instrumental profile of the spectrograph may be important, 
especially in the case of astronomical spectra, which are usually obtained with a 
relatively wide slit. The instrumental profile of the spectrograph (designated 
Litt GIII, BL84) was estimated by deriving the mean profile of the four lines 
AA 3920, 3922, 3927 and 3930 in the iron arc comparison spectrum. ‘The central 
part of the profile was obtained from lightly exposed spectra, while the wings were 
derived from heavily exposed spectra. ‘The mean instrumental profile is given in 
Table I. The linear dispersion on the tracings is about 0-586 A/inch in this wave- 
length region. For each individual line no measurements were made over 
wave-length regions where a blending line might be suspected to occur. The 
asymmetry of the wings of the instrumental profile reflects the asymmetrical 
distribution of the intensities of the grating ghosts, one of which lies very close to 
the main line. With this instrumental profile estimates were made of the instru- 
mental distortion that might be present in the observed profiles of the stellar 
features. It was found that the central parts of the observed profiles of narrow 
lines had been distorted by the spectrograph and that the true line profiles are 
deeper, by approximately 15 per cent of the continuum, than the observations 
indicate. In view of the large observational uncertainties in the final line profiles 
(see below), it hardly seems significant to attempt to correct the observations for 
instrumental broadening. However, instrumental broadening is an important 
source of error in the numerical value of any result dependent on the measured 
apparent residual intensities of the narrow spectral lines. 

The tracings described above give the variation with wave-length of the 
intensity of radiation of the K star alone (from spectra obtained during totality), 
of the K star plus the B star (from spectra obtained when the stars were far from 
eclipse), and of the K star plus B star plus chromospheric features (from spectra 
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obtained during atmospheric eclipse). In each case these intensities for the wave- 
length range 3915 <A<3950A are expressed relative to the peak intensities at 
A3915 and A3950. It is evident that if the mean brightnesses at these wave- 
lengths of the B star and of the K star are b and k, respectively, and that if b+k=1, 
then if the intensity tracing for the K star be reduced in vertical scale by a factor k 
and subtracted from the tracing for the K star plus B star, the result will be an 
intensity tracing giving the B-type spectrum. Since the spectrum of the B star 


TABLE I 
Instrumental profile for Litt GIII BL84 spectrograph 


Intensity 
25°8 


a* Intensity 
—13 foe) 
12 o'2 
II o'4 
o'8 
°'9 
13 
1°8 
1°5 
1‘o 
2°8 
10°2 
29°3 13 
75°6 14 
100°0 15 
71°! +16 


= 

° 
~~“ = ; 
HK OC CO'I ANUAwWDn & 


_ 
N 


9 
8 
7 
6 
5 
4 
3 
2 
I 
° 
I 


* Units of x are 0-1 inch on a tracing magnified 200 x from the plate in the direction of 
the dispersion. Negative displacements are to the violet; positive to the red. Tracings 
were made from a plate having a projected monochromatic slit-width of 0-0020 inches. 


is known to have constant intensity (no absorption lines except perhaps a narrow, 
weak K-line) over the wave-length interval 3915 <A<3950A, the result of the 
subtraction should be a straight line at intensity 1—k. Similarly if the K-type 
spectrum reduced in vertical scale by the factor k be subtracted from a tracing of 
the K star plus B star plus chromospheric features, the result will be a tracing 
giving the B-type spectrum as a straight line at intensity 1 — k with absorption lines 
at the wave-lengths of the chromospheric features. The relative brightness of 
the B star to the K star will be. 

a=(i—k)/k. (1) 


Intensity tracings of the B-type spectrum and of the B-type spectrum plus 
chromospheric features have been obtained by this method from the intensity 
tracings of the composite spectrum of 31 Cygni at various phases. The sub- 
traction of the tracings is done semi-automatically using the intensitometer of this 
Observatory in the manner described by K. O. Wright(z). The ratiok by which to 
reduce the K-type spectrum is found by trial. A series of tracings of the K-type 
spectrum reduced by a series of factors k is prepared; by superposing these 
tracings in turn on a tracing of the composite spectrum one may quickly determine 
which value of k yields a result such that the B-type spectrum in regions outside 
the chromospheric lines falls along a line of constant intensity 1—k. ‘The tracing 
of the B-type spectrum will vary in a random manner about the constant value 
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1—k which represents the continuum of the B star, for the irregularities due to 
grain will be somewhat increased over those for a single tracing. The mean 
tracing from two grating plates obtained during totality was used as the standard 
intensity tracing of the K-type spectrum. The resultant intensity tracings of the 
B star alone, the B star plus chromospheric features and the K star alone are shown 
in Figs. 1 and2. The result for the B star alone (phase +258 days) is the mean of 
two plates, while those for the B star plus chromospheric features are each from 
one plate. In spite of the irregularities due to grain and to small differences in 
placing the continuum on each original tracing, it is clear that the main features of 
the growth in complexity of the chromospheric spectrum are portrayed. 
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Fic. 1.—The chromospheric spectrum of 31 Cygni at various phases during ingress in the 
atmospheric eclipse. The intensity is expressed in terms of the continuum of the B star. The 
phase is recorded in days from second contact, 1951 August 12. The spectrum of the K star has 
been subtracted, leaving only the spectrum of the B star as it is modified by passage through the 
outer layers of the atmosphere of the K star. 


If during the atmospheric eclipse because of obscuration or occultation the B 
star appears fainter with respect to the K star than when the stars are far apart, the 
empirically determined factor k will increase tok’. ‘The apparent relative bright- 
ness of the B star to the K star will then be 


= fa=(1—R')/k’. (2) 


Hence one will be able to determine the loss of light of the B star relative to the 
mean intensity of the K star in the wave-length range 3915 <A<3950A. The 
results are given in Table II. The phase during ingress is counted in days from 
second contact, which is assumed to have occurred (3) on 1951 August 12 (U.T.). 
During egress the phase is measured in days from third contact, 1951 October 13 
(U.T.). The second column of Table II gives the phase, the third gives the 
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empirically determined ratio k or k’, the fourth gives « or «’, while the fifth gives f, 
the factor by which the light of the B star has been reduced. The uncertainty in 
the determination of the factor k is about +0°03. This results in an uncertainty of 
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Fic. 2.—The chromospheric spectrum of 31 Cygni just before totality and the spectrum of the 
K star during totality. 
continuum of the B star. The intensity of the K-type spectrum is expressed in terms of the 
apparent continuum of the K star in the spectral range AA 3915-3950. 


Date 
=. 
July 31°3 
Aug. 2°3 
Aug. 4°3 
Aug. 63 
Aug. 7°3 
Aug. 9°3 
Aug. 10°3 
Aug. 12°8* 
Oct. 15°2 
1952 June 25'9* 


The intensity of the chromospheric spectrum is expressed in terms of the 


TaBLe II 
Light curve of the B star at AA 3915-3950 
Phase k 


—12 days 0°30 
0°30 
0°30 
0°33 
0°36 
0-40 
060 
1°00 
0°40 
0°30 


* The mean date of two exposures. 


some 20 to 25 percentina. However, even sucha large uncertainty does not mask 
the result that the light received from the B star begins to decrease some time near 
phase —6 days, and that the loss of light increases rapidly with the approach of 


second contact. 
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The equivalent widths of the chromospheric Fe 1 lines AA 3920, 3922, 3927 and 
3930 were measured at each phase, and are presented in Table III. The mean 
equivalent width of the four lines is given in the last column. At phases —3, —2 
and +3 days this mean is for AA 3920, 3922 and 3927 only, since at these phases a 
moderately strong chromospheric line of Ti 1 blends with Fe 1 43930. These 
mean equivalent widths are plotted against time in Fig. 3, which shows a general 
similarity to comparable curves for ¢ Aurigae (4). 


TaBLe III 
Equivalent widths of chromospheric Fe I lines 


Date Phase A3920 A3922 A3927. A3930 Mean 


U.T. Days A A A A A 
July 31°3 —12 0°04 0°03 0°05 0°02 0°04 
Aug. 2°3 —10 08 ‘Il 13 “15 12 
Aug. 43 — 8 08 13 16 “15 13 
Aug. 6°3 — ‘IO "15 ; 18 16 
Aug. 7°3 _ ‘29 ‘29 . ‘35 "31 
Aug. 9°3 - "51 ‘57 : "55 
Aug. 10°3 — “51 56 : 55 
Oct. 15:2 4 “44 "49 ; 54 
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Fic. 3.—The growth of the chromospheric Fe I lines of 31 Cygni. The mean equivalent width 
of the lines 4A3920, 3922, 3927 and 3920 is plotted against the phase expressed as days before 
second contact. 


Discussion.—The growth of the chromospheric spectrum as the B star 
approaches the photosphere of the K star is illustrated in Figs. 1 and 2. At 
phase —12 days the Fe 1 lines AA 3920, 3922, 3927 and 3930 are barely visible 
although the K line of Ca 11 is quite strong. However, two days later the Fe 1 lines 
are quite distinct. These lines increase in strength as the light from the B star 
traverses deeper and deeper layers of the atmosphere of the K star. By phase —3 
days many other chromospheric lines have appeared. The strongest of these 
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lines, identified from Miss Davis’ list of lines in the spectrum of § Pegasi (5) and: 
from the Revised Multiplet Table (1), are listedin Table IV. They arise from energy 
levels with excitation potentials from 0 to 3 volts, and they strengthen greatly as 
the line of sight approaches the photosphere of the K star. Their relative 
intensities are roughly as in the K-type star except for Fe 1 43918, which is weak 
relative to Cr 143919. Also the relative intensities of Ti 143915 and Cr 13916 
appear to be reversed in the chromosphere to what they areinthe K star. However, 
in view of the observational uncertainties of the intensity tracings shown in Figs. 1 
and 2, it is not certain that this apparent reversal is significant. 


TABLE IV 


Chromospheric lines in the region 3915 < A < 3933 A 


—* Element 


Multiplet pe Phase from second contact 
sow E.P. é 
number at which first seen 
eV 

391588 Til 
3916°24 Cri 
3917°18 Fe I 
3919°16 Cri 
3920°26 Fel 
392102 Cri 
3921°42 Til 
3921 °80 Zri 
3922°91 Fel 
3924°53 Til 
3925 °24 Vi 
3927 °92 Fel 
392864 Cr 
3929°88 Til 
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393366 Call 
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The light curve of the B star derived from the spectra as set forth in Table II 
gives the brightness of the B star at various phases relative to the mean brightness 
of the K star at A3915 and A3g950. Although these measurements are inferior in 
accuracy to those that may be obtained with modern photometric equipment, they 
serve to indicate that the light from the B star diminishes by an appreciable amount 
before occultation of the disk of the Bstar by the photosphere of the K star canoccur. 
Miss Vinter Hansen (6) has studied the radial-velocity variations of 31 Cygni and 
from her work it appears that the B star is moving at a speed of about 35 km/sec 
with respect to the K star. General considerations suggest that the spectral type 
of the B star is about B3 V or Bs V and that its diameter probably does not exceed 
1osolar radii. Since a star moving at a speed of 35 km/sec takes 2-3 days to traverse a 
distance of 10 solar radii, it follows that the occultation of the disk of the B star will 
probably not commence prior to phase —3 days. However, our observations 
indicate that by this time approximately one-third of the light of the B star had been 
extinguished. It is suggested that at about phase — 6 days the line of sight begins 
to pass through a more dense region of the atmosphere of the K star, and that this 
region causes an extinction of the light of the B star. This extinction is not due to 
the occurrence of chromospheric absorption lines, for the present results are based 
on a comparison of the flux from the B star with that from the K star in spectral 





No. 5, 1954 chromospheric lines during the 1951 eclipse of 31 Cygni 565 


regions free from chromospheric line absorptions. It must be caused by some 
process varying little with wave-length in the region 3915 <A<3950A. It is 
assumed that no significant intrinsic fluctuations in the light of the K star (or of the 
B star) occurred during the period of observation. Unpublished observations by 
C. M. Huffer* and observations by F. Bradshaw Wood (7) indicate that the light 
of the K star was constant during totality. Also unpublished observations by 
E. F. Carpenter* indicate that during egress there was extinction in the ultra-violet 
similar to that found here for ingress. Wood (7) also presents photometric data 
for egress, suggesting that atmospheric extinction plays an important part in the 
eclipse. 

The conclusion that before first contact the B star may be partially obscured 
by continuous absorption in the atmosphere of the K star is supported by observa- 
tions of ¢Aurigae. Fracastoro (8), from a spectrophotometric study of low- 
dispersion plates of ¢ Aurigae obtained during the 1947-48 eclipse, finds for 
¢ Aurigae that the light of the B star has been dimmed by some 25 per cent at A 3933 
when first contact occurs. Roach and Wood (9), from a detailed study of photo- 
metric data on the eclipse of ¢ Aurigae, also conclude that the B star is dimmed, 
particularly at shorter wave-lengths, by chromospheric absorption in the atmo- 
sphere of the K star. They conclude that the dimming is chiefly due to the growth 
of chromospheric absorption lines, although they cannot rule out the possibility 
that it is due to continuous absorption in the atmosphere of the K star. Photo- 
metric observations of the partial phases of an eclipsing system such as 31 Cygni, 
made with narrow-band filters which do not pass wave-length regions containing 
strong chromospheric features whose intensities would change greatly with phase, 
would be of great value in determining accurately the amount of extinction at each 
wave-length and its overall wave-length dependence. Only a limited interpreta- 
tion can be made from data such as those discussed by Roach and Wood because of 
the relatively long spectral range covered by each filter. 

It may be readily shown that the relative brightness of the B and the K star may 
be estimated from the residual intensities in the totality spectrum and in the 
composite spectrum at phases far from eclipse. We have 

_ Rex— Rx 


. nag Rex’ 


(3) 


where Rp is the residual intensity in the composite spectrum and R, is the residual 
intensity in the totality spectrum. Both intensities are measured relative to the 
apparent continuum. Equation (3) may be written in terms of equivalent widths 
(W): om, 
a ie (4) 
It is assumed that the intensity in the B-type spectrum is independent of wave- 
length and that the absorption lines appearing are due solely to the K star. 

From the apparent residual intensities at the centres of the lines AA 3920, 3922, 
3927, 3930 and 3944, the mean value of « is 1-7. From the equivalent widths of 
regions centred on these lines and about 1 A widex=1-9. The result found above 
from the subtraction of the intensity tracings isa=2-3. The value 2-3 is felt to be 
more correct than the other values, for it is found by a method in which the 
distortions due to instrumental broadening are not important. The effects of 


* Private communication to A, McKellar, 
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instrumental broadening on the present spectra are sufficient to cause the 
differences in« noted here. The use of equations (3) and (4) to estimate « is not to 
be recommended except for spectra insignificantly distorted by the spectrograph. 
Theoretically equation (4) should not be affected by instrumental broadening, 
but practically it is difficult to be sure of integrating the whole broadened profile 
when W is measured in a many-lined spectrum such as the present one. 

The variation of the mean equivalent widths of the Fe1 lines AA3920, 3922, 
3927 and 3930 with time during ingress is shown in Fig. 3. It is significant that a 
rapid increase in strength occurred after phase — 6 days when the light of the B star 
first became dimmed. The lines do not change appreciably in equivalent width 
after phase —3 days. ‘This is probably because the lines have become saturated 
and the curve of growth has become more or less flat. Reference to studies of 
K-type spectra by K. O. Wright (10) shows that the curve of growth begins to 
flatten out when the equivalent width is about o-1 A. This growth in strength of 
the Fe 1 chromospheric lines is similar in many respects to that of the Ca 11 K 
line (11). 

If it be assumed that the Fe 1 lines fall on the linear part of the curve of growth, 
as seems likely when the equivalent widths of these lines are small, then the equiva- 
lent width will be proportional to the number of atoms in the line of sight. If the 
density of neutral Fe atoms is p (r) at a distance r from the centre of the K star, then 
the number of atoms in an element of length ds along a column 1 cm? in cross- 
section which passes within a minimum distance p from the centre of the star is 


dN =p(r) ds =p(r)p cosec® ¢ dd, (5) 


where ¢ is the angle between the radius r passing through the element of volume 
and the line of sight. The total number of atoms in the line of sight is 


N(p)=P | p(r) cosec*¢ a, (6) 


Pi 


where ¢, and 4, are the limiting values of ¢. If R be the outer radius of the zone 
containing neutral Fe 1 atoms, that is, the outer radius of the stellar atmosphere, 
then 


N(p)=2p rf p(r) cosec® ¢ dd. 
~ sin-' p/R 
For purposes of illustration let us assume that 


p(r) = C/r? = C/p? cosec? ¢. 
Then 
N(p)= == 


| -= —sin“! F- 
Pp sin-' p/R 


The equivalent widths of the aa absorption lines formed along lines 
of sight with minimum distances p, and Pp, will then be in the ratio 
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Let the unit of distance be the daily motion of the Bstar. Then at any date the 
quantity p is the number of days from mid-eclipse. It will be assumed that the 
K-type atmosphere extends to twice the radius of the photosphere, as seems likely 
from observations of the Ca 11 chromospheric lines. Then R=62. The ratios 
of equivalent widths, W(p)/W (39), where W(39) is the equivalent width at phase 
—8 days, computed from equation (10), and the observed ratios are given in 
Table V. The computed ratios agree reasonably with the observed ratios in the 
first part of the table but between phases — 6 days and — 2 days the observed ratios 
are much greater than the computed ratios. Since allowance for a possible 
departure of the curve of growth from linearity would further reduce the computed 
ratios, it is clear that after phase — 6 days the line of sight must pass through a region 
in the atmosphere of the K star where the density increases much more rapidly 
than the assumed law p~1/r?._ We conclude that the break in the curve relating 
strength of the chromospheric lines to phase represents a real discontinuity in the 
density of the outer atmosphere of the K star. ‘The density in the region traversed 
by the line of sight between phases — 6 and o days is significantly greater than that 
of the outer region of the atmosphere. 


TABLE V 
The mean relative strengths of the Fe 1 lines at various phases 


Computed Observed 


Phase p W(p)/W(39) W(p)/W(39) 


Days 
—12 43 082 o'3I 
—10 4I o'9I 0-92 
8 39 1°00 1‘0o 
37 IIo 1°23 
36 1°16 2°39 
34 1°28 4°23 
33 1°34 4°23 


It has been pointed out that the line profiles shown in Figs. 1 and 2 are distorted 
to some extent by the spectrograph. The true profiles of the Fe 1 chromospheric 
lines are somewhat deeper and a little narrower than these profiles. However, 
a measure of the mean half-width of these profiles will give an upper limit to the 
half-width of the line absorption coefficient. The mean half-width at all phases 
corresponds to a Doppler width of 2okm/sec. Since the thermal motion of Fe 
atoms at 7'= 3500 deg will give a Doppler width of 1 km/sec to the line absorption 
coefficient, it is clear that there is some turbulent broadening of the profiles of the 
chromospheric Fe 1 lines. This conclusion is interesting in view of the results 
for Ca 11 at higher atmospheric levels during egress in late October, November and 
December 1951. McKellar, Odgers, Aller and McLaughlin (12) have shown ina 
preliminary study of the K line of Ca 11 that part of the material in the outermost 
regions of the atmosphere of the K star was in rapid motion a good deal of the time 
during egress, although the effect was scarcely, if at all, detectable for the Ca 11 
lines during ingress. ‘The amount of turbulent broadening of the Fe 1 chromo- 
spheric lines of 31 Cygni is a little larger than the results for  Aurigae. For this © 
star O. C. Wilson and H. Abt (13) have found by curve-of-growth methods that 
the turbulent velocity is about 8 km/sec for neutral atoms and about 12 km/sec for 
ionized atoms, ‘The general experience in the analysis of stellar spectra has been 
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to find that the turbulent velocity derived from curves of growth is equal to or less 
than that derived from the shapes of the line profiles (14). 
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TWO-COLOUR PHOTOELECTRIC PHOTOMETRY OF 
GO CYGNI 


M. W. Ovenden 
(Communicated by the Director of the Cambridge Observatories) 
(Received 1954 August 23) 


Summary 


Photoelectric light-curves in the blue and yellow have been obtained for 
GO Cygni from observations with the Cambridge pulse-counting photometer 
during the period 1950 August to 1951 October. The paper reports, 
analyses and discusses these observations. Both light-curves show strong 
asymmetrical terms in the variation of light between eclipses; the observations 
show that the processes giving rise to these terms suffer intrinsic variations. 
The observed differences of magnitude at maxima are compared with 
O’Connell’s statistical conclusions. Elements for the binary are derived by 
Russell’s method, Pearce’s spectrographic observations being utilized. The 
comparison of the light-curves in the two colours gives evidence that the 
observed secondary spectrum (type Ao) is due to reflection, the spectral 
type of the secondary if unaffected by reflection being later than F8. The 
assumption that the radiation of the secondary arises from a point close to the 
primary star serves to reconcile the deduced masses with the general mass— 
spectral type relationship. 





1. Introduction.—GO Cygni (BD + 34°4095) was chosen for observation with a 
view to testing, under observing conditions at Cambridge, the suitability of the 
pulse-counting photoelectric photometer for the photometry of eclipsing variables. 
‘The sensitivity of the photometer set an upper limit about 11™ to the magnitude of 
the variable at primary minimum, and the behaviour of the Cambridge sky 
suggested that the star’s period should be less than 3 days, but not close to a multiple 
of aday. ‘lo examine the value of photometry with the given equipment in the 
problem of eclipsing binaries, it was decided to take a star for which satisfactory 
visual and photographic observations were available, in addition to spectrographic 
observations. GO Cygni was chosen because Pierce’s visual observations (1) 
showed that the star was a critical case in which more extensive and accurate 
observations might have been expected to provide a much more accurate know- 
ledge of the elements of the system. In the event, this promise has not been 
fulfilled, owing to an intrinsic variability of the light-curve, but the photoelectric 
observations have enabled an analysis of these variations to be made. 

The variability of GO Cygni was discovered photographically by Schneller (2) 
in 1928. Kukarkin (3) showed that the star was a short-period eclipsing binary, 
and Szezyrbak (4) gave the first elements of the light-curve, considering the system 
to be of the Algol type. Kukarkin’s visual observations from 1929 to 1931 (5) 
(which were analysed by Kopal (6)) showed variation of the 8-Lyrae type, with 
equal maxima and no displacement of the secondary minimum. Kordylewski’s 
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visual observations (7) also show equal maxima, but give the phase of the secondary 
minimum as 0-509. A series of visual observations by Beyer (8) gave a strong 
eccentricity, but the observations are few. A set of photographic observations by 
Iwanowska and Dziewulski (9) gave no eccentricity and symmetrical minima, but 
the maximum following the primary eclipse about o™-1 brighter than the maximum 
following the secondary minimum. Observations by Payne-Gaposchkin (10) 
a few years later gave a similar excess of the first maximum in the photoblue, but 
an excess of about the same amount of the second maximum in the photored, both 
colours showing strong asymmetry of both minima. 

The most complete series of observations, with discussions, are those of Liau 
(photographic) and Pierce (visual, polarizing photometer). Liau (11) finds the 
primary minimum asymmetric, in the sense that ascent is more rapid than decline, 
and the first maximum 0-075 brighter than the second maximum. _ Pierce (1, 12) 
finds no asymmetry of minima, and equal maxima. Neither Pierce nor Liau 
found any evidence of eccentricity. 

The only published spectrographic observations are those of Pearce (13), from 
which radial velocities were obtained for both components. The spectral types of 
the components were given as Bgn and Aon, but, as will be seen later, the radiation 
of the secondary is due almost entirely to reflection. 

2. Experimental details —The present paper reports and discusses some 700 
individual photoelectric measures of GO Cygni made between 1950 August and 
1951 October with the pulse-counting photoelectric photometer designed by 
Yates (14) and mounted on the Huggins 15-inch visual refractor of the Cambridge 
Observatories. A graphical investigation of the resultant spectral transmissions 
of the two filters used (Wratten 5 yellow and Wratten 39 blue), the photomulti- 
plier (from a sensitivity curve for a typical photocell) and the objective (from data 
for a typical 15-inch visual object glass (1§)) gave maximum sensitivity at 5250A 
(yellow) and 4200 A (blue). The blue transmission curve is sensibly symmetrical, 
with 50 per cent transmission at 3840A and 4540A, but the yellow curve is 
markedly asymmctric, with 50 per cent transmission at 4920A and 5770A. 
These figures are only approximate, for in any case variations in sensitivity from 
one cell to another may occur. 

At the beginning of the programme, two comparison stars were chosen, viz. 
HD spectral 

type 

A 35° 4197 8:2 (7) A2 

E 34° 4098 7°9(2) K2 
but star A appeared to show some variability, and only the comparisons of the 
variable with star E have been used. ‘The correlation of the deviations of the 
observations of the variable from a mean curve with the phase of the variable, and 
the lack of any systematic effects of differential extinction give assurance that the 
use of only one comparison star has not caused any loss of accuracy in the measures. 

Observations were made in two colours on the comparison star, then in two 
colours on the variable, and then in two colours on the comparison star again. 
Sky readings were taken between every other pair of comparisons, except on nights 
of bright moonlight, when they were made between every pair of comparisons. 
A single cycle giving two comparisons in each colour took about 10 minutes, 
each exposure being of half a minute. ‘The one minute exposures more usual 
with this photometer proved too slow for the variable at certain phases. 


Star BD No, HD mag. 
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For each determination of the variable, the count on the sky and on the 
comparison star was found by linear interpolation of adjacent readings. After 
subtraction of the sky reading, the ratio of counts on the variable and the com- 
parison star was converted to a magnitude difference. Although the Cambridge 
programme of standard photometry provides a curve in each colour for the cor- 
rection of the observed magnitudes to the NPS system, according to colour 
index, no attempt has been made to reduce the original observations of the 
variable to the NPS system, since the correction curves are obtained from 
observations of supposedly single stars. Not only does the spectral distribution 
of the binary not correspond to that of a single star of the same mean colour 
index, but furthermore the relative contributions of the two components vary 
throughout the eclipses. This is a fundamental difficulty of eclipsing variable 
work, if the photometry is not monochromatic, and for this reason caution must 
be exercised in any attempts to combine in a single light-curve observations 
made with more than one photometer. These considerations do not affect the 
validity of the analysis of a light-curve made entirely in one photometric system. 
The whole analysis of GO Cygni has therefore been made in the photometric 
systems of the photometer, and the standard colour correction curves have been 
used only in the discussion of the colour of the binary when the contribution of 
each component can be estimated. 

For reference and comparison with other observations, a series of observations 
ot the comparison star E were made on two nights with a view to placing it upon 
the NPS scale. Four stars of Kapteyn Selected Area No. 81, stars BD 16°2466, 
13° 2643, 15°2561 and 16°2514, were compared with E. ‘This selected area 
was chosen because it could be observed at equal zenith distances with the field 
of the variable, and for the four chosen stars satisfactory magnitudes on the NPS 
scale are available from the standard photoelectric programme. Comparisons 
of E with BD 13° 2643 were made over a range of time including the instant of 
equal altitude, and the magnitude differences were plotted against the difference of 
the secants of the zenith distances of the two stars; this provides a correction 
curve against relative altitude which can be used for the comparisons of E against 
the other three stars of the selected area. The magnitude of the Selected Area 
stars on the NPS scale from the Cambridge standard photometry were taken as: 


Star Ipg Ipv 
BD 13° 2643 g°17 8-32 
BD 16° 2466 8-62 7°77 
BD 15° 2561 10°19 g*00 
BD 16° 2514 9°23 8-23 


‘he reduced magnitudes for each night, with the mean deviations for each night, 
and the adopted mean magnitudes for the comparison star E are: 


Date Ipg Mean 6 Ipv Mean 6 
m m m m 
1952 May 15-16 8-78 0:03 7°84 0°03 
1952 May 22-23 8-74 0-05 7°78 0°04 
Adopted 8-76 7°81 


3. Elements of the light-curve.—Light-time corrections were found from 
Prager’s tables (16), graphically interpolated to the nearest minute. 
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Pierce gave the elements of the light-curve as 
JD 2426 509-461 + 09-717 763 2 E. 


‘The best-determined minimum of the present series of observations, at corrected 
time JD 2 433 861-499, gave a period of 04-7177626+0%-0000005. An attempt 
was made to improve this value of the period by utilizing the whole of the light- 
curve. The observations were assembled on a period of 04-717 7626, and for 
each night the sum of the individual deviations along the time axis from a mean 
curve drawn freehand through the observations was divided by the total number 
of observations for that night, and plotted against the number of periods from an 
arbitrary zero. Each point was then weighted according to the number of 
observations and the quality of the night, and mean points found for three ranges 
of E, the number of periods from the arbitrary zero. The points so found lay 
close to a straight line, giving a correction of + 04-000 003 5 to the period, but the 
scatter of the individual points, due mainly to phases outside the primary minimum, 
gave no confidence in this correction, and the period of 04-717 762 6 was adopted. 
A comparison of this figure with previous values of the period by earlier observers 
does not give any definite evidence of a change in the period of the binary. 

4. Formation of normal points.—For forming normal points an objective 
criterion of the photometric quality of a given night is desirable. ‘The stability 
of the response of the pulse-counting photometer provides such a criterion in the 
internal agreement of the counts of the comparison star alone. Each observation 
of the comparison star was compared with the reading obtained by linear inter- 
polation of the adjacent counts of the comparison star, and the mean deviation 
in magnitudes, 7, of these comparisons found for each night. For a large number 
of such comparisons, with Gaussian errors, the empirical weight to be attached to 
observations on a given night would be proportional to 7~?, but this gives too 
strong a weighting in view of the uncertainty of on a given night due to the 
limited number of comparisons available. An empirical weight in the range 
1 to 5 was therefore adopted for each observation, the weight factor being 
approximately proportional to 7! for the night of observation. 

Normal points were made for ranges of about 15 minutes throughout the 
cycle, the normal values being taken as the arithmetic means of the weighted 
individual magnitude determinations. ‘The empirical weight, w,, of a normal 
point was taken to be the sum of the empirical weights of the observations 
contributing to the normal point. According to Kopal (17), the intrinsic weight, 
w;, of a normal point is given by +/w,; 0c dl/d(sin? @) for errors constant on an 
intensity scale, and by \/w; oc l-!. dl/d(sin? 0) for errors constant on a magnitude 
scale, / being the intensity of the binary in terms of the intensity at maximum. 
Observations by Yates on a standard laboratory source with the same photometer 
showed that the instrumental mean deviation is of the order of o™-005, whereas 
the mean deviations found in practice in these observations and the standard 
photometric programme are of the order of o™-oz. Thus the errors are due to 
variation of sky transparency, and would be approximately constant on a magnitude 
scale. An intrinsic weight w, oc /-*. (dl)*/(d(sin? @))? was therefore adopted, the 
final weight for each normal point being taken as w=w,. u;. 

The normal points and their weights are given in Tables | and II. Light- 
curves on a magnitude scale are given in Figs. 1 and 2. 
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"TABLE | 
Normal points (yellow curve) 
Vu, Vw Phase Var-Comp. 
P m 
6:0 0°5000 +0°734 
58 ‘5157 -0°723 
4°2 "5290 + O°747 
3°4 5447 O"705 
31 "5585 +0°698 
2°0 5738 +0696 
1'7 "5894 0°702 
13 6045 +0616 
o'9 6167 0628 
°°6 6312 +0610 
o"4 ‘6460 + 0°575 
06 6585 -+-0+562 
o°4 ‘6747 +0°588 
o'4 ‘6885 + 0°533 
°'4 "7047 +0'526 
o'2 "7192 +O°525 
o-2 7331 rO°§25 
o'2 7463 0518 
orl “7615 +0°503 
orl 7764 + 0°486 
0°3 ‘7847 + 0°537 
o'2 8165 0°532 
o-3 ‘8339 0°573 
o"7 ‘8491 + O'571 
O'4 8620 +0"590 
o'7 8776 +0648 
o'7 ‘8906 +0656 
o'8 “9066 FOr715 
o"9 "9195 0°751 
°'9 9358 +0821 
1'0 ‘9497 +0°904 
o'8 ‘9632 -o-958 
o'4 ‘9782 + 1°52 
o'2 0'9978 | 1°087 


No. of 


Var-Comp. ron 


° 
in 
a 


m 
+1056 
+ 1°046 
+0985 
+0°935 
+0825 
+o-764 

0'716 
+ 0°673 
+0618 
+ 0°603 
-0°564 

0°541 

-0°545 
+0508 

0°505 

7 O°507 
+ 0°500 
F0°495 
+0489 

0-500 

O°512 

oO°515 

-0°547 

0°527 

-0°544 
F0°595 
+0°605 
-0'576 
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0690 
0693 
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Fic. 1.—The light-curve from the normal points in the yellow system. Ordinates give the 
magnitude of GO Cygni minus the magnitude of the comparison star, on the system of the photometer. 
An approximation to international photovisual magnitudes, based on a colour correction corre- 


sponding to the mean colour index of the variable, may be obtained by taking the zero of the magnitude 
scale as +-7™-93. 
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Tasie II 
Normal points (blue curve) 


No. of 
obs. 
? m P m 
0:0068 +0°025 2°6 O'5150 —0°371 
O215 —0o-'oo! 2°4 "5284 —0°371 
"0334 — 0°032 20 "5422 —0°355 
‘0475 + —O°107 20 5565 —0°393 
0642 —0°255 2°4 ‘5734 —0°394 
‘0819 —0°325 2'2 5879 —0°405 
0954 —O'414 2°6 6020 —0'436 
‘1089 = —0°405 24 ‘6147, —0°434 
1273 0478 2°8 6317 —o-'476 
-1634 ~0°485 2°8 . 6466 —O'514 
‘1805 —0'527 2"0 " ‘6591 —o'501 
‘1952 0526 2°4 . 6744 —0*500 
"2093 0548 2°8 , 6889 —0o'4g91 
"2245 0°576 2°2 , “7041 —0'524 
°2375 —o°582 2°2 ‘ ‘7172 —~0°534 
‘2531 4 —0'558 2°8 , ‘7315 "547 
‘2684 -0°544 2°8 . ‘7448 —0°533 
2821 —0°536 2°4 ; “7605 —o'558 
2963 —o'561 2°8 ° “7760 —0'595 
"3223 0°535 3°3 ‘7856 =—-0°573 
"3398 —0'542 2°8 ‘ 8188 —0'556 
‘3590 —0o'502 8414 —0'541 
-3706 0-448 8641 —o'5il 
"3838 0°457 8772, — 0468 
‘3978 0487 ‘8918 = =—0°399 
"4135 0°478 “9080 —0368 
"4295 —0°484 9216 —0'283 
‘4427 0°425 ‘9347. —0'219 
"4571 —0°429 9485 —o'168 
‘4701 0°391 9631 —0:'076 
4856 0-366 ‘9779 +0020 
0'5001 —0°379 0'9925 +0058 
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Fic. 2.—The light-curve from the normal points in the blue system. Ordinates give the 
magnitude of GO Cygni minus the magnitude of the comparison star, on the system of the photometer. 
An approximation to international photographic magnitudes, based on a colour correction correspond- 
ing to the mean colour index of the variable, may be obtained by taking the zero of the magnitude 
scale as +8™-go. 
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5. Deviations of the individual observations from the mean curve, in magnitudes.— 
In addition to the supposedly random variation of sky transparency, four sources 
of deviation from the mean curve suggest themselves, viz. (i) lack of constancy 
of the comparison star, (ii) differential colour extinction due to the difference in 
colour of the variable and the comparison star, (iii) errors in timing of the observa- 
tions, and (iv) intrinsic variations of the light-curve of the variable. 

The possible effects of differential colour extinction may be investigated by 
making use of the stability of the response of the photometer from night to night 
over long periods. The absorption of the atmosphere is reflected both in the 
observed colour index of the comparison star (defined by the ratio of the counts 
in the two colours) and in the actual count in any one colour, e.g. the blue, which 
is the more sensitive to changes in atmospheric absorption. ‘The absorption 
being a logarithmic effect, a linear relationship should exist between the logarithm 
of the count of the comparison star in the blue, and the colour index of the com- 
parison star. Such a linear relationship was found to exist, with a mean deviation 
of about 2 per cent, for a given colour index. However, no correlation of the 
deviation of the variable from the mean light curve (5) with the logarithm of the 
adjacent count of the comparison star was apparent; differential colour extinction 
does not, therefore, contribute appreciably to the observed deviations. 

A cursory glance at the plots of the individual observations shows that the 
deviations from the mean curve show some systematic variations with phase of 
the binary motion in both colours. To examine the variation of deviation with 
phase it is necessary to allow for the quality of the night on which a given observa- 
tion was made. Although the magnitude of the mean deviation of the observations 
of the variable due to atmospheric effects will depend upon the photometric 
quality of the night, the ratio 5/n may be expected not to vary from this cause, 
5 being measured in magnitudes from the mean curve drawn freehand through 
the normal points. 

For the blue curve the mean vaiue of 8/y for equal ranges of |sin 20| was 
plotted against | sin 20 |, and a strong correlation was found (29). The variation 
of light between eclipses in the blue shows a sin 26 term with a large negative 
coefficient. Liau also found a variation of error with phase, although he made 
no detailed analysis. 

Treated in the same way, the yellow observations showed no correlation with 
|sin 20|. ‘The yellow curve has a term in sin 26 for the variation of light between 
eclipses which is of doubtful significance, but it has an important term in sin 6 
with a positive coefficient. Treated in a similar way for the variation with 
|sin @|, the yellow deviations showed a correlation in the opposite sense to that 
of the blue observations with |sin 20| (29). Since in each’ colour the major 
variation of intensity outside eclipses is of the cos 20 type, any changes in mean 
deviation due to changes in the measured magnitude difference between variable 
and comparison star would show the same correlation with phase in both colours. 
It may thus be concluded that the processes causing the sin @ term in the yellow 
and the sin 2@ term in the blue are subject to intrinsic variations of the order of 
a few hundredths of a magnitude. 

If the effects of these intrinsic variations were removed, 5 would be of the order 
of 7, but somewhat smaller than 7 since the interval of interpolation of the original 
observations for obtaining the magnitudes of the variable was less than the 
interval of interpolation between adjacent observations of the comparison star 
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used to obtain 7. ‘The intercept in the blue correlation line is 5/7=~0-6, and in 
the yellow correlation line is 5/n=o-8; this agrees with expectation if variations 
of the comparison star were making no significant contribution to the deviations. 

To investigate possible errors of timing, a correlation of 5/n with the slope of 
the blue-light curve (in magnitudes) was sought. None was found; errors of 
timing are thus negligible. 

6. The variation of the colour of the variable with phase.—Taking the colour 
index of the comparison star E as zero, the relative colour index of the variable 
was found from each pair of normal points (blue and yellow). The phase of the 
combined colour normal was taken to be the arithmetic mean of the phases in the 
two colours without any relative weighting, and the weight of each colour normal 
was taken to be the sum of the empirical weights of the individual observations 
in the two colours. The weighted colour normals were then combined into new 
colour normals for ranges of 0-05, these final colour normals being plotted in 
Fig. 3, which shows the variation of the colour of the variable wich phase. 
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Fic. 3.—The variation with phase of the colour of GO Cygni. Ordinates give the variation of 
colour index on the systems of the photometer, but the zero of the ordinate scale has been adjusted 
to give an approximation to international colour indices. 


7. The variation of light between eclipses——The magnitudes of the normal 
points were converted into intensity ratios. A provisional analysis of the variation 
between eclipses by Merrill’s graphical method (18) indicated that terms up to 
sin 20 should be included in the analysis in both colours. ‘The variation of light 
between eclipses was therefore represented in both colours by the series 

l= Ay+ A, cos 6+ A, cos 20+ Az sin 6+ A, sin 20. 

A difficulty frequently encountered with light-curves of the 8 Lyrae and W Ursae 
Majoris binaries is the finding directly from the light-curve of the duration of the 
eclipses. ‘The slope di/d(sin® @) is more sensitive to the moments of first and 
last contacts, and some suggestion of a discontinuity of slope in the plot of this 
quantity with phase in the blue appeared at +0":12; no discontinuity was 
detected in the similar yellow plot. Least-squares solutions for the 
A-coefficients front the weighted normal points of the blue light-curve were made 
for various assumed eclipse durations, the results being shown in Table III. 
Solutions 2, 3 and 4 are similar, and give last contact at a phase close to 0-125, 
but as there is a danger of including some points actually in the primary eclipse 
in solution 4, solution 3 was accepted as being the best possible representation 
of the variation of light between eclipses. The weighted yellow normals were 
then subjected to a least-squares analysis for the A-coefficients on the single 
assumption of an eclipse duration of +0"-150, and these coefficients are also 
shown in Table III. 

The major cause of variation between eclipses in both colours is ellipticity 
of components, shown by the cos 2term. The yellow curve shows some 
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Eclipse I 

Duration + o0?-200 

Blue 

Ay +0°965+0°043 
+0'045+0'013 
-0'018+ 0-048 
+0002 + 0:003 


Two-colour photoelectric photometry of GO Cygni 
Tase III 


Variation of intensity between eclipses 
Least-squares solutions for the A-coefficients 


2 
+ OP '175 


+0°901 + 0'009 
+-0'003 + 0°005 
—o'085 +0012 

0'004 + 0°005 


3 
+ oP 150 


+ 0916+ 0:006 
+-0°008 + 0004 


—0°069 + 0-009 
—0°OOI + 0°002 


4 
+ oP 125 


+ O°915 + 0°003 
+ O°'O1O+ 0°004 
-0'072 + 0°007 
-0:006 + 0:003 


5 
+ o?-100 


+-0'899 + 0°003 


0°014+ 0°004 


—0:098 + 0007 


0-001 + 0'003 


—-0°023 + 0:006 
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0°025+0:003 —0'023+0'002 -—o-010+0-002 —0-006+0-014 
| O°912+0°005 
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—0:076 + 0-008 
-++-0°012+ 0'002 
—0:'006 + 0-002 


Intensity 








N, 
N N 


















































O'O of fof fos.) fo 
Phase 


Fic. 4.—The rectified light-curve in the yellow system, in intensity units. 
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Fic. 5.—The rectified light-curve in the blue system, in intensity units. 
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reflection, but an amount much less than would be expected from the properties 
of the binary system. The reflection in the yellow may be partially disguised 
by a variation with cos @ in the opposite sense corresponding to the significant 
“negative reflection” present in the blue curve. Strong asymmetric terms 
are present in both colours. 

8. Rectification.—In the absence of any satisfactory theory for the causes of 
the asymmetric terms, and the “ negative reflection’’, no precise rectification 
of the light-curves is possible. As far as possible, the rectification followed the 
method laid down by Russell (1g). The reflection term in the yellow was 
treated as normal reflection, and the “negative reflection” term in the blue 
removed by subtraction. The sin @ term in the yellow was removed by sub- 
traction, following Dugan’s hypothesis of a difference in brightness of the 
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advancing and receding hemispheres(20). The terms in sin 26 and the ellipticity 
terms were removed by division. 

The rectified curves are given in Figs. 4 and 5. The secondary minimum in 
the blue is due almost entirely to ellipticity, and disappears on rectification. 
A small residual secondary minimum remains in the yellow curve. 

9. Computation of the elements of the system.—The absence of a significant 
secondary minimum to the rectified blue curve renders this curve virtually 
useless for the determination of the elements of the system, and the analysis was 
made with the yellow curve only. 

An analytical computation of the elements, following Kopal’s procedure (17), 
was attempted, but the iterations were found to diverge both for the assumption of 
primary transit and the assumption of primary occultation. The rectified curves 
do not, therefore, correspond exactly to the mutual eclipses of two spherical 
stars, and their failure to do so may be attributed to the presence of large terms 
which have no satisfactory theoretical basis. An approximate solution for the 
elements is all that is possible, and this was made on Russell’s simplified model, 
using the tables of Merrill (21) and Zessewitsch (22). 

The circumstances of the eclipse proved to be such that the light-curve is 
insensitive to the value of the ratio of the radii, k, and no determination of k was 
possible from the light-curves. A value of k was therefore found from the relative 
intensities of the two spectra and the relative depths of the two minima in the 
yellow curve. Pearce (quoted by Pierce (1)) gave L,=0-07+0-01, which 
yields k=0-67+0-05. For convenience in using the tables, a value of k=0-65 
was adopted, and curves for limb-darkening coefficient x =0-2, 0-4, 0-6, 0-8 and 
1-0 were computed. ‘The curve for complete darkening was considered to give 
the best representation of the observations. 

No determination of elements was possible from the blue curve owing to 
the absence of secondary minimum after rectification, but the blue curve was 
found to be consistent with the same assumptions of k and magnitude of eclipse 
as were found for the yellow curve. 

The remaining photometric elements were computed by Russell’s method. 
The factors NV for conversion from photometric to geometric ellipticities were 
taken from the Roche model (23), 


N=(15+x)(1+y)/(15 — 5x). 

The limb-darkening coefficients x were taken as 1-0 for the yellow (from 
observation) and o-8 for the blue (from Miinch and Chandrasekhar’s table (24) 
for an At star, the earliest type for which figures are given). ‘These limb-darkening 
coefficients are not entirely satisfactory. For a Bg star, electron scattering may 
make a significant contribution to the opacity. Also, according to the theory of 
radiative transfer for a grey atmosphere, the limb-darkening coefficient should 
increase with decreasing wave-length. The calculations of Miinch and 
Chandrasekhar for opacity due to neutral hydrogen atoms and negative hydrogen 
ions (24) show that the linear law of limb-darkening is a poor approximation, and 
particularly so in the photovisual region, for which they give no limb-darkening 
coefficients. ‘The geometric ellipticities were calculated with y=1 for want of 
a better approximation (18). 

The remaining elements may be computed from Pearce’s spectrographic 
observations. However, the radial velocity observations may be influenced by 
the fact that the secondary spectrum is due to reflection (see later) and will arise 
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from the hemisphere of the secondary which faces the primary. ‘The observed 
radial velocities will thus correspond, not to the centre of the secondary star, 
but to a point on the line of centres which is displaced towards the primary star. 
If this displacement is ga, where a=(a,+4,) is the separation of the two stars, 
then o<g<r,. With r,=0-333 as determined, elements have been derived on 
three assumptions as to the value of g, viz. 0-00, 0:25 and 0-33. The derived 
elements are given in Table IV. 
TABLE IV 
List of elements 
Photometric 
Ratio of radii, k =0°67 
Radius of larger star, 7, =0°513 
Radius of smaller star, r, =0'333 
Ellipticity, € =0'129 (yellow obsns.) 
=0°'134 (blue obsns.) 
=0'131 (adopted mean) 
Max. eclipsed area of smaller star, «0-761 
Inclination, i a 99° 
Eccentricity, e =0°'00 
Spectrographic 
I II III 
§=0'25 §=0°33 
Radius of orbit of larger star,a,(x10%°km) = 1°38 1°38 1°38 
Radius of orbit of smaller star, a, ( x 10° km) 1°62 2°62 3°12 
Separation, a (=a,+4) ( x 10° km) 3°00 4°00 4°50 
Radius of larger star, r, (< 10° km) 1°54 2°05 2°30 
Radius of smaller star, r. ( x 10° km) 0997 1°33 1°50 
p M, 
Ratio of masses, im 085 0°53 0°44 
. 1 
Mass of larger star, M, 1-120 3:20 4:80 
Mass of smaller star, M, 095 O 1'70 210 


10. The reflection effect.—The reflection will be confined to the secondary 
component, since this component contributes less than a tenth of the total light 
of the system. ‘The observed spectral type of the secondary, Ao, may well be due 
to reflection of the radiation of the primary. Eddington (25) has given a formula 
for the amount of radiation reflected at full phase from a point-primary, and this 
formula should give an order-of-magnitude estimate of the amount of reflected 
radiation in the present case. When applied, Eddington’s formula gives a 
luminosity due to reflection of 0-ogL,, where L, is the luminosity of the primary. 
This is of the same order of magnitude as the observed luminosity of the secondary 
spectrum, which may therefore be attributed primarily to reflection. It is not 
possible to deduce the amount of reflected radiation from the photometric reflec- 
tion effect in the yellow curve, for it may be masked by “ negative reflection” 
corresponding to the observed “ negative reflection” in the blue curve. 

On the assumption that the quality of the radiation from any surface is that 
of a black body at some effective temperature, an estimate of the average spectral 
type of the area of the secondary which undergoes eclipse may be obtained, 
provided that the relative response of the photometer in the two colours to radiation 
from stars of various spectral types is known. ‘The relative response may be 
found from the international colour indices and bolometric corrections on the 
assumption that the photometric systems of the photometer do not depart from 
the international systems, It is probable that this approximation would be 
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adequate for the present purpose, but a slightly better approximation can be 
obtained from the colour correction curves which convert the systems of the 
Cambridge standard photometry to the NPS systems, since the Cambridge 
standard photometry involved photoelectric observations with the same photometer 
as was used for the present investigation, in addition to photographic observations. 

With the relative response of the photometer to surfaces of various spectral 
types known, the nature of the reflection effect may be examined. ‘The relative 
brightnesses of the inside surfaces of the two stars may be computed for various 
assumptions of mean spectral type for the eclipsed area of the secondary and 
compared with the observed ratio of minima in the light-curves (uncorrected for 
reflection). In the yellow curve the observed ratio of depths of minima corre- 
spond toa mean spectral type of F8 for the eclipsed area of the secondary, taking 
the primary star to be of type Bo. 

In the blue curve it is possible to set only an upper limit of o-10 for the 
ratio of the depths of minima. However, this upper limit is precise enough to 
enable some conclusions to be drawn concerning the nature of the reflection 
effect. If the radiation of the primary is reflected without change of quality, 
then the colour of the eclipsed area of the secondary should correspond to spectral 
type Bg, this assumption giving a ratio of depths of minima in the blue of 0-17. 
If, on the other hand, the radiation of the primary is absorbed, and re-emitted at 
a lower effective temperature, then the colour of the eclipsed area of the secondary 
should correspond to spectral type F8, this assumption giving a ratio of 0-09. 
It may be concluded that the primary radiation is absorbed and re-emitted at 
a lower effective temperature. ‘This is in agreement with the conclusion of 
Pismis (26) from a statistical consideration of many eclipsing systems. 

Although the mean spectral type of the eclipsed area of the secondary is F8, 
Pearce’s spectrographic observations show that some of the eclipsed area has an 
effective temperature corresponding to spectral type Ao. ‘This suggests that 
(i) the effective temperature must change rapidly over the inward hemisphere of 
the secondary, (ii) the observed radiation from the secondary must arise from a 
region close to the primary star, and not from the centre of the secondary star, 
and (iii) the spectral type of the secondary if unaffected by reflection would be 
later than F8. 

11. The masses of the components.—These conclusions are confirmed by 
a comparison of the deduced masses of the components with the general mass— 
spectral type relationship. For this comparison the data given by Russell and 
Moore (27, Fig. 2) have been used. 

The mean mass for spectral type Bg is about 3-30 + 1-40. — Solution II 
of Table IV, corresponding to g=0-25, gives a mass of the primary close to the 
mean mass for this spectral type. Solution III is on the edge of the observed 
spread, but Solution I is outside the observed spread. (It should be noted that 
trigonometrical as well as spectroscopic binary data are available.) ‘The mass- 
ratio deduced from the ellipticity is 0-65, which agrees best with Solution II. 
‘The mass of the secondary component in Solution II corresponds to a spectral 
type about Go. The assumption that the secondary spectrum arises from a 
limited region close to the primary star thus serves to reconcile the deduced 
masses of both components with the general mass-spectral type relationship. 

12. The nature of the asymmetric terms.—Asymmetric terms of the form 
Ay sin @ and A, sin 26 are frequently found in the variation of light between 
eclipses with close eclipsing binaries, The sin 6 term is the most easily recognized 
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from the light-curve, since it produces a difference in the heights of the two 
maxima. O’Connell (28) has made a statistical analysis of the occurrence of 
differences in the heights of maxima of light-curves of the 8 Lyrae and W Ursae 
Majoris types. The difference of magnitude in the yellow curve of the present 
observations is taken to be +0-015, following O’Connell’s sign convention, 
where Am is taken to be positive if the maximum following the primary minimum 
is the brighter. (A value of Am taken directly from the light-curve is used in 
preference to one deduced from the size of the coefficient A,, since the expression 
of the variation of light between eclipses as a Fourier series is entirely arbitrary.) 
In Table V the deviations of the observed Am from O’Connell’s regression lines 
are given, together with the probable errors of Am given by O’Connell. 
TABLE V 
Difference of magnitude of maxima (yellow obsns.) compared with O’Connell’s statistics 
Am (obs) — Am (O’C) p.e.(O’C) 

Ellipticity, «€ —o-o16 +0016 

Separation, a (soln. IT) +0006 +0014 

Ratio of radii, k —0'040 +0016 

Mass of primary, M, (soln. IT) —0'OI7 +0:016 

Ratio of densities, © (soln. IT) —o'oll +0019 

2 

The agreement is satisfactory except for the ratio of the radii. The failure 
of the observed asymmetry to agree with the regression line for the ratio of the 
radii may be related to the conclusion arrived at above that the radiation of the 
secondary is from a limited area close to the primary star. The adopted value 
of k=0-65 was found from the relative intensities of the two spectra on the assump- 
tion that the secondary radiated uniformly; a failure of this assumption will 
cause the radius of the secondary to be underestimated. A ratio of about 0-80 
would suffice to bring the observed difference in the heights of the maxima into 
agreement with O’Connell’s statistics. 

O’Connell also examined the variation of Am with wave-length, and concluded 
that Am decreases with increasing wave-length, and may become negative for 
photovisual and photored observations. Of the stars used in arriving at this 
conclusion, GO Cygni is included with observations by Payne-Gaposchkin (10). 
The present observations are inconsistent with those of Payne-Gaposchkin, 
although it should be pointed out that Liau’s photographic observations do show 
an appreciable positive Am. The present observations are rather in agreement 
with Mergantaler’s conclusion that Am increases with increasing wave-length. 
In any case, the relationship may well be obscured by other asymmetries, and 
it is significant that one of the exceptions to his rule cited by O’Connell is RS CVn, 
where other asymmetries are marked. The possibility that the magnitudes of 
the asymmetries of GO Cygni vary with time is strongly suggested by acomparison 
of the various series of observations at widely separated epochs. 

The correlation of deviations from a mean curve with the asymmetric terms 
shows that the processes giving rise to these terms suffer intrinsic variations 
even during a single period. The significance of these observations has been 
discussed by the present author elsewhere (29). Following Mergantaler (30), 
the sin @ variation is attributed to absorption in a system of gaseous envelopes 
surrounding the binary, and the sin 2@ term is attributed to a similar cause. 
The asymmetry which increases with increasing wave-length (in this case the 
sin @ variation) is supposed due to absorption by H™- ions (cf. Mergantaler, 
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loc. cit.). If a model similar to that of Kuiper (31) for 8 Lyrae is adopted for 
GO Cygni the H~ absorption is important in the expanding shell system. For 
the close shell which constitutes a common envelope of the two stars, absorption 
by hydrogen atoms on the Balmer level will be the more important process. 
Since a few per cent of the radiation received in the blue system of the photoelectric 
photometer used in the present investigation has a wave-length less than that 
of the Balmer series limit, this may explain the importance of the sin 20 variation 
in the blue curve and its absence in the yellow curve. It would be of value to 
make photometric observations of a system like GO Cygni in two restricted 
wave-length regions, one on each side of the Balmer discontinuity. It may be 
remarked that the flare-like variation found by Huruhata for U Pegasi (32) shows 
a strong wave-length dependence of this type. 
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COLOUR-MAGNITUDE ARRAYS IN THE CLUSTERS 
TOMBAUGH 5, NGC 7789, NGC 1528 AND NGC 2682 


V. C. Reddish 
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Summary 


Magnitudes and colours of stars in four clusters have been obtained 
from simultaneous photographs on blue- and red-sensitive plates, using the 
twin refractor of the University of London Observatory at Mill Hill. The 
magnitudes have been reduced to the (B, V) system of Johnson and Morgan. 

The colour—magnitude arrays show the cluster Tombaugh 5 to be of 
similar type to Praesepe. NGC7789, previously classified as a galactic 
cluster, is found to be a globular cluster. Trumpler’s classification of 
NGC 1528, a normal galactic cluster, is confirmed. NGC 2682 is known 
to be peculiar. The author considers it to be a transition type, between the 
normal Population I and Population II types. It is tentatively classified 
as of Population I on the basis of the existence of a narrow Hertzsprung 
Gap, a basis which the author suggests should be used as the main criterion 
to differentiate between clusters of the two types of population. 

A continuous sequence of types of colour—magnitude arrays of clusters is 
arranged, in which the main variation along the sequence (apart from the 
brightness of the stars) is the width of the Hertzsprung Gap. This variation 
is probably a result of differences in the proportions of the stars’ masses 
contained within their convective cores. 





During the years 1953 and 1954, observations were made with the Radcliffe 
twin refractor of the University of London Observatory, Mill Hill, to determine 
colours and magnitudes of stars in four clusters. 

Because of difficulties in guiding for long periods near the pole, and the 
well-known disadvantages of using the stars in that region as a standard system, 
the (B, V) system of Johnson and Morgan (1) was adopted as the standard. 
The B—V colours in which the results are expressed were determined from 
colours (denoted B—R) obtained from the difference between blue and red 
magnitudes, the longer base-line of the latter resulting in smaller proportionate 
random errors. 

Equipment.—The 18-inch visual refractor attached to the 24-inch photographic 
refractor was adapted for photography in 1953 September. It thus became 
possible to photograph simultaneously on red-sensitive and blue-sensitive plates. 
Guiding was done by means of an eyepiece attached at the side of the plate-holder 
on the 24-inch refractor. Photographs were taken on Ilford Astronomical Zenith 
(and in some cases on Kodak Oa-O) plates on the 24-inch refractor, which was 
stopped down to 18 inches to make the images of comparable density to those 
on the red-sensitive plates. Red-sensitive Kodak Oa-F plates were used, 
together with a Chance OR2z red glass filter, on the 18-inch “visual” refractor. 

Exposure and development.—All exposures were of one hour’s duration, the 
limiting visual magnitude being 15-2. Plates were developed for three minutes, 
in the pairs as exposed, in Kodak D 19b developer at 18 deg. C. 
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Measurement and reduction.—No photometer being available, diameter 
measurements of the images were used to obtain magnitudes, and a Hilger 
‘Spectrum Plate’”’ measuring machine was adapted for the purpose. All plates 
were measured twice—in the north-south and east-west directions. Magnitudes 
were estimated in the case of images which were too faint for diameter measure- 
ments. ‘The colour-magnitude diagrams therefore extend to the limiting 
magnitude of the equipment. Satisfactory relations were found between 
magnitude and square root of image diameter. 

Extinction.—It was not possible to make satisfactory measurements of 
atmospheric extinction and estimated extinction coefficients were employed. 

The standard system 

Visual magnitudes.—All visual magnitudes were standardized by comparison 
with plates of one-hour exposure on the clusters Praesepe and IC 4665, for which 
magnitudes on the (B, V) system have been given by Johnson (2, 3). A linear 
relation between red magnitude and square root of image diameter was found 
and is represented by the equation 


R= —2:231/d,+ 22-74. (1) 
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Fic. 1.—Observed minus standard visual magnitudes, against visual magnitude. 
Each dot represents a magnitude determined from a single plate. 
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Fic. 2.—Observed minus standard visual magnitudes, against colour index. 
Each dot represents a magnitude determined from a single plate. 


The relation between the colours was then found to be 
B—V=o0-70(B—R) (2) 

and hence the visual magnitudes are given by 

V = —2:23+/d,+0-42(B—V)+22-74, (3) 
where d, is in units of o-o1 mm, and corrections for extinction have been made. 
Values of “observed minus standard” visual magnitudes are plotted against 
magnitude and colour in Figs. 1 and 2. There is a slight non-linearity in the 
plot against colour, of the order of a few hundredths of a magnitude. T he use 
of non-linear relations to eliminate this would have greatly increased the labour 
of reduction, and was not considered worth while in view of the very small 
increase in accuracy which would result. 
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Blue magnitudes.—Because of a misunderstanding regarding supplies of 
photographic plates, no Kodak Oa-O plates were available to continue the 
programme after photographs had been obtained on the clusters Tombaugh 5 
and NGC 7789. At that time, no exposures had been made on the standard 
(B, V) system, and previous work had to be adapted for the reductions, as 
explained in the next paragraph. The remainder of the programme was carried 
out using Ilford Astronomical Zenith plates. 
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Fic. 3.—Square root of image diameter against blue magnitude (corrected for extinction). 
The units of d are o-o1 mm. In addition to the curve, equation (5) is used for magnitude 
determinations. 
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Fic. 4.—Observed minus standard blue magnitudes, against blue magnitude. 
Each dot represents a magnitude determined from a single plate. 
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B-V 
Fic. 5.—Observed minus standard blue magnitudes, against colour index. 
Each dot represents a magnitude determined from a single plate. 


(a) Clusters Tombaugh 5 and NGC 7789.—The blue magnitudes were obtained 
from photographs on Kodak Oa-—O plates, using an equation derived from a 
simple relation between square root of image diameter and magnitude. This 
relation, based on the (P, V) system and obtained from photographs of the 
North Polar Sequence, also on Kodak Oa—O plates, was found in earlier work (4). 
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The equation, corrected for extinction, is 

B= —1:25/(d—10)+0:03(B— V)+ 19°17, (4) 
where d is in units of oor mm. The equation holds down to d=o-15 mm 
(about B=16-5). The details of the manner in which this equation was 
obtained are too long for publication in this paper, but because some ap roxi- 
mations had to be made in its derivation it cannot be expected to give magnitudes 
as reliable as those in the following section (5). 

(6) Clusters NGC 1528 and NGC 2682.—The blue magnitudes were obtained 
from photographs on Ilford Astronomical Zenith plates, and standardized by 
comparison with photographs on similar plates of the clusters Praesepe and 
1C 4665. The relation between square root of image diameter and blue magnitude 
was found to be non-linear; it is shown in Fig. 3, corrected for extinction. The 
symbol B, is used for the blue magnitudes because there is an additional 
correction for colour given by the equation 

B=B,—0-14(B—V)+0-05. (5) 
The method of reduction was to read off B, from Fig. 3 and then apply equation (5) 
to get B. Values of “observed minus standard” blue magnitudes are plotted 
against magnitude and colour in Figs. 4 and 5. 

Errors.—The internal probable errors of the magnitudes and colours deter- 
mined from one pair of plates (red and blue) have been found by repeated 
measures to be as follows: 


V 10m 12m I 4m 15m 
Probable error in V + 0-04 + 0-06 +0™-06 + 0M-07 
Probable error in B—V +0™-07 +0™-08 +o™-08 +o™-10 


‘Two pairs of plates were obtained for each of the clusters NGC 7789, NGC 1528 
and NGC 2682, and three pairs for the cluster Tombaugh 5. 

The main probable sources of systematic error are: 

(1) Uncertainty in the accuracy of equation (4). This affects only the clusters 
Tombaugh 5 and NGC 7789. 

(2) The use of a linear relation—equation (3)—for reducing the visual magni- 
tudes ; the true relation probably differs from linearity by a few hundredths 
of a magnitude. 

(3) The use of estimated extinction coefficients. This probably makes the 
zero point of the colours uncertain by about a tenth of a magnitude. 


The clusters 


Tombaugh 5.—R.A. 3" 42™-9 Dec. + 58° 54’ sles 
Gal. Long. 112° Gal. Lat. +4°:8 94 


This small faint cluster was discovered by Tombaugh in 1941 (5). The 
colour—magnitude array is given in Fig. 6. It appears to be a normal galactic 
cluster of Population I, with a few yellow giants. The six brightest main-sequence 
stars are separated from the rest of the main sequence by a gap of about 0-4 
magnitude. One of them, at V=11™-76, B—V = +0™-07, must really be very 
blue if it is a cluster star, because comparison of the main sequence with that 
of Praesepe (2) and the position of the yellow giants indicate a colour excess of 
about 0:35 magnitude. This would give for the star a true colour B— V = —o™-28. 
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There is a very red star at V = 11™-76, B— V=2"™-35; correction for colour-excess 
gives the true B—V=2™-oo. It is almost certainly a cluster member—if it were 
more distant than the cluster and its red colour were due to interstellar absorption 
behind the cluster it would be much fainter than it is. 

If the absolute visual magnitude of the yellow giants is +0-3 (cf. Praesepe 
—a similar type of cluster), the distance modulus is m—M=12-30. Total 
absorption A,,,=3E,_, (6). Hence the total absorption for this cluster 
A,i,= 1™-05 and the distance modulus corrected for absorption is 11™-25. The 
corresponding distance is 1800 parsecs. 

NGC 7789.—R.A. 23" 52™0 Dec. +56° 10’ hi on 

Gal. Long. 83° Gal. Lat. —5°-2 {*9 

This cluster was listed by Shapley (7) as a galactic cluster, and by Trumpler (8) 
as a galactic cluster of type 2—3a, signifying probably more yellow giants than 
main-sequence stars. 

The colour—magnitude array is given in Fig. 7. It is typical of the bright red 
giant branch of globular cluster arrays, showing that the cluster is of Population II 
and that the previous classifications as a galactic cluster are incorrect. The cluster 
is unusually open for a globular cluster—indeed it is only for the faintest stars 
measured that clustering becomes obvious at all. It is probable that the stars 
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Fic. 6.—Colour—magnitude array of Tombaugh 5. Fic. 7.—Colour—magnitude array 
of NGC 7789. 
near the centre could be resolved down to much fainter magnitudes than have 
been reached in these observations. Comparison with the colour-magnitude 
arrays of various globular clusters obtained by Arp (9) (with allowance for the 
difference in the colour systems) suggests that NGC 7789 has a colour excess 
of about o™-2; hence the total absorption 4,,,=o™-6. Assuming that the 
brightest stars on the giant branch have an absolute visual magnitude of — 3-5 
(from which it follows that the horizontal branch, usually including the RR Lyrae 
variables, is just below the limiting magnitude of the diagram), the distance 
modulus is m—M=14™"-9. The distance modulus corrected for absorption is 
14™-3. The corresponding distance is 7200 parsecs, and hence the cluster is 
about 650 parsecs south of the galactic plane. The luminosity function for the 
giant branch is given in Fig. 10. It has been shown by Sandage (10) that the 





588 V. C. Reddish, Colour-magnitude arrays in the clusters Vol. 114 


luminosity function for stars between M,,,= +1 and M,,,= +4 in the globular 
cluster M3 is closely represented by the formula 

logy) $(M) =0-4M,, + const., 
which would be expected from some evolutionary theories. Fig. 10 shows that 
the formula closely represents the luminosity function for the giant branch in 
'NGC 7789 for stars brighter than M,,,, = —0°5. 
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Fic. 8.—Colour—magnitude array of NGC 1528. Fic. 9.—Colour—magnitude array 
of NGC 2682 (= M67). 


NGC 1528.—R.A. 4" 7™-6 Dec. + 50° 59’ hr = 
Gal. Long. 120° Gal. Lat. + 1°-2 9 

The colour—magnitude array of Fig. 8 confirms Trumpler’s classification (8) 
as a galactic cluster of type 1—2b~a, signifying that most of the brightest stars 
are on the main sequence and are of spectral types B to A. 

A colour excess of 1™-o is indicated by the position of the main sequence, 
giving a total absorption A,,,=3™-0. Comparison of the main sequence with 
that of Praesepe suggests that the absolute magnitude zero corresponds to 
V=11™1. Hence the distance modulus corrected for absorption is 8™-1 and 
the distance 420 parsecs. 

The star at V=g™-82, B— V =o™-54, must be very blue if it is really a cluster 
member. Correction for colour excess gives a true colour B—V = —o™-46 at 
M,= — 13. 

NGC 2682.—R.A. 8" 45™-8 Dec. + 12° 11’ ha on 

(=M67) Gal. Long. 184° Gal. Lat. +33°-4{ °9 

Listed by Shapley (7) as a galactic cluster, and by Trumpler (8) as a galactic 
cluster of type 2—3a, NGC 2682 was investigated by Ebbighausen (11) who 
pointed out that the Hertzsprung—Russell diagram bore some similarity to those 
of globular clusters. Further work on the cluster has been done by Vander- 
linden (12) and by Becker and Stock (13) who classified it (on the basis of the 
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colours of the brightest clusterstars)asa Population I globular cluster. The present 
measures extend about a magnitude fainter than those of the latter investigations. 

The colour—magnitude diagram, Fig. 9, is similar to those of globular clusters, 
but there are striking differences. These can be seen most easily by comparing 
a schematic diagram for NGC 2682 with one of a typical globular cluster, M3 (14). 
Since the absolute magnitudes of the stars in NGC 2682 are not known, the 
ditferences depend on how the clusters are compared. 
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Fic. 10.—The luminosity function for the giant branch of the globular cluster NGC 7789, 
full line. The dashed line gives the curve to be expected if the number of stars per unit luminosity 
range is inversely proportional to the luminosity. 


(1) If the stars in the region V=11™-5, B—~V = + 0™-3 (Fig. g) are assumed 
to be the same as those on the horizontal branch of M3 about M,=o-0, 
the comparison is as in Fig. 11. In this case the main sequence of 
NGC 2682 has already been reached by these observations, but the 
‘“‘turn-off” of this main sequence takes place at about M,= +2°5, 
compared with M,= +3-6 in the case of M3. Also, the vertical giant 
branch in NGC 2682 would appear to be o™-7 redder than the main 
sequence at the “turn-off” point, the corresponding figure in M3 being 
about o™-2. 

(2) If the brightest red giants in NGC 2682 are supposed to be equal in 
luminosity to the brightest in M3, the comparison is as in Fig. 12. Here, 
the “horizontal” branch in NGC 2682 is nearly two magnitudes brighter 
than in M3, and if the “‘turn-down” at V = 14™ (Fig. g) is supposed to 
be the start of the vertical giant branch, it is both fainter and bluer than 
the vertical branch in M3. 

The first comparison appears more likely to be near the truth as regards 
comparing equivalent sequences. The author considers, however, that the cluster 
may represent a stage of transition from Population I to Population II. Despite 
the strong resemblance to colour-magnitude arrays of Population II, there are 
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characteristics typical of Population I—a narrow Hertzsprung Gap at V = 13™-6, 
the “bunching” of the red giants bordering the gap (at V=13™-3) and the 
comparatively small number of red giants brighter than the bunching. 

It is possible to arrange the colour—magnitude diagram of NGC 2682, 
together with diagrams of the clusters M3 (14), NGC 752 (15), Praesepe (2), 
M38 (16) (reduced to the (B, V) system (4)) and the cluster in Perseus (17) 
into a sequence in which the main variable along the sequence (apart from the 
brightness of the stars) is the width of the Hertzsprung Gap (Fig. 13). NGC 2682 
has been placed in what appears to be its most natural position in the sequence. 
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2682 2682 
Vv V 
1OF 1OF 
lif M3 Vif ai 
4My=O as 
12h "ce 
13F 13F M3 
¢ My=O0 
14F 14 fg rx 
ISP Sr NS 
0-90 +1.0 B-V 0.0 +10 py 
Fic. 11.—Comparison of M3 (full line) Fic. 12.—Comparison of M3 (full line) 
and NGC 2682 (shaded). and NGC 2682 (shaded). 


See the text. See the text. 


Comparison of the colours of the stars in NGC 2682 with those to be expected 
from the spectral types published by Ebbighausen (11) shows good agreement 
and indicates that there is no interstellar reddening of the cluster. This is not 
unexpected in view of the cluster’s high galactic latitude. If the zero point of 
the absolute magnitudes indicated by Fig. 13 is correct, the distance modulus 
is 12™-o and the corresponding distance is 2500 parsecs. 

“he decreasing width of the Hertzsprung Gap along the sequence of diagrams 
in Fig. 13 may be expected from theoretical considerations. Along a sequence 
of stars of decreasing masses the proportion of the mass contained within the 
convective core decreases until finally it vanishes in stars of about the solar mass. 
This is because electron scattering becomes relatively less important in the 
opacity and the proton—proton reaction becomes relatively more important in 
the energy generation as the masses decrease. Both these changes act to reduce 
the proportion of mass contained within the convective core (18, 19, 20). It 
has been suggested by the author (21) that when the hydrogen in the core has 
been consumed in the nuclear reactions, evolution to the red-giant stage will be 
rapid and will thus cause a gap in the sequence of the stars on a colour-magnitude 
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array—the Hertzsprung Gap. The red-giant models of Sandage and Schwarz- 
schild (22) show that as the proportion of mass contained within the hydrogen- 
exhausted core increases, the radius of the star increases, and its surface temperature 
falls. It is probable, therefore, that the width of the Hertzsprung Gap depends 
on the proportions of mass contained within the convective cores of the brightest 
main-sequence stars. Briefly, this is because the larger the core the more it must 
contract after hydrogen exhaustion—and hence the more the envelope expands 








+4 








' = 

0.0 +1.0 

Fic. 13.—Comparison of various colour—magnitude arrays, 
showing the narrowing of the Hertzsprung Gap. 


B-V 


Cluster in Perseus. 
M38. 

Praesepe. 

NGC 752. 

NGC 2682. 

M3, outline only. 


Only those portions of the arrays bordering the Hertzsprung Gaps are given. 


and the redder the star becomes—before the temperature at the surface of the 
core rises high enough for the carbon—nitrogen cycle to act there and slow down 
the rapid evolution. Thus the width of a Hertzsprung Gap on the colour- 
magnitude array of a cluster may depend on the masses of the brightest main- 
sequence stars in the cluster. Since the masses of the stars decrease with decreasing 
brightness, the width of the Hertzsprung Gap would, according to this theory, 
decrease with decreasing brightness as it does in Fig. 13. This explanation is 
of course tentative and requires further investigation. 
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The classification of clusters according to Population must be based on 
some easily recognizable criterion. The author suggests that the presence or 
absence of a Hertzsprung Gap is the best criterion to use. It suffers from the 
obvious disadvantage that the ability to detect narrow gaps depends on the 
accuracy of the colour measurements, but this difficulty will probably apply to 
any criterion based on the use of colour-magnitude arrays. One of its principal 
advantages is that it is closely connected with a structural characteristic of the 
stars—the presence or absence of convective cores—which is itself connected 
with the age of the stars. 

On the basis of this criterion, NGC 2682 would be classified as Population I. 

Addendum.—Since this paper was presented, a colour—magnitude array of 
NGC 2682 determined by Johnson (23) has been published. Fig. 9 appears to 
agree well with Johnson’s diagram except for stars fainter than V =14", where 
the colours of this paper are slightly bluer. 


The author acknowledges with pleasure the help and encouragement he 
received from Professor C. W. Allen during this investigation. 


University of London Observatory, 
Mill Hill Park, 
London, N.W.7: 
1954 September 9. 
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STAR-STREAMING FROM THE PROPER MOTIONS OF THE 
A-TYPE STARS OF THE BOSS GENERAL CATALOGUE 


T. R. Tannahill 
(Communicated by the Director, University Observatory, Glasgow) 
(Received 1954 November 4) 


Summary 


1. The proper motions of 8676 stars of spectral types B8—A3 in the Boss 
General Catalogue have been analysed, according to the two-streams 
hypothesis, by Eddington’s method. 

2. The galactic coordinates of the vertex are (365°-5, —o''5). 

3. The equatorial coordinates of the solar apex are (264°°8, +-31°°5). 

4. The ratio of the number of stars in the drifts is 1°54. 

5. Comparison of the results with those from later types reveals large and 
irregular variations in the constants of star-streaming which are not to be 
accounted for by accidental error. It is suggested that inadequacy of the 
two-drift hypothesis, or non-homogeneity of the groupings of the data, may 
be the explanation of these. 





1. Ina previous paper (1) the proper motions of 22 374 stars of spectral types 
As-M in the General Catalogue of 33 342 Stars by Benjamin Boss (2) have been 
analysed on Eddington’s two-drift hypothesis in three spectral groups and four 
zones of galactic latitude. A later paper (3) analyses the motions of 1536 O- and 
B-type stars of the same catalogue. ‘The present paper analyses the motions of 
8 676 stars in the spectral range B8—A3, and completes the determination, by this 
method, of the constants of the star streams, and the directions of the solar motion 
and the vertex of star-streaming, from the data of this catalogue. 

2. The treatment of the observational data follows almost exactly that 
described (1) for the stars of types A5s—M, and reference may be made to this 
paper for the system of assigning weights to the proper motions, for the systematic 
corrections applied and for the grouping into regions. Stars belonging to known 
clusters have been rejected—z2o of the Pleiades group, 46 of the Ursa Major 
group and g of the Taurus cluster. As before, diametrically opposite regions 
have been combined, and (in cases where the precaution seemed necessary) 
corrections have been applied to the position angles of the proper motions to allow 
for the convergency of the meridians in regions of high declination. 

The results of the analysis are given in Table VI.* In this table, kV, andhV, 
are the projections of the space-velocities of the drifts on the tangent-plane at 
the centre of the region, 6, and @, the position angles of the directions towards 
the drift apices and N, and N, the numbers of stars allocated to the drifts. As 
with the As5—M stars, the regions round the Drift I apex offer difficulty in analysis, 
particularly in the determination of the Drift II constants. A further difficulty 

* Tables VI to IX appear at the end of the paper. 
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in this investigation is the higher concentration of the stars in the galactic plane, 
the regions in the higher galactic latitudes being very sparsely populated—in 
Zone IV, comprising the six regions round the galactic pole, the total weight of 
stars does not exceed 100 in any region. Analysis has, however, been possible 
in all regions. 

3. Solutions for the drift apices have been found for three groupings of the 
data—a solution for all the stars involved in the discussion, and solutions for low 
galactic latitude (obtained by combining the regions of Zones I and II of the 
previous paper) and for high galactic latitude (Zones III and IV). ‘The mean 
galactic latitude for the low latitude group is 13°-7, and for the high latitude 
group 51°°6. ‘The results are given in Table VII, in which X, Y and Z are the 
linear components of the space-velocity of a drift relative to the Sun, referred to 
the usual equatorial axes, /W is the resultant space-velocity and A and D are the 
equatorial coordinates of the apices. The decrease of AW, with increasing 
galactic latitude, already noted in the previous paper (1) on the A5—M stars, is 
quite apparent. In order to establish this variation without ambiguity, a further 
solution for the A5—M stars has been made, in galactic latitude groups similar 
to those used for the A-type stars (namely, combining Zones I and II and Zones III 
and IV). The results, given in Table I, show that this variation with galactic 
latitude holds for limited spectral ranges. It is not, therefore, a reflection of 


TABLE I 
Variation of hW, with galactic latitude 
B8-A3 As5-F5 F8-G5 Ko-M p.e.* 
Zones I, II 1°742 1°934 1-756 1°491 +0°04 
Zones III, IV 1°518 1661 1°416 1°299 +0°04 


any similar variation in spectral type. A comparison of the results from the 
spectral groups of the previous paper with those of the present paper (‘Table II) 
shows the surprising feature that the decrease of AW, with advancing spectral 
type, which is so clearly marked in the A5—M types, is not maintained in the A-stars 
of the present paper. A systematic change in the Drift II apex is now, however, 
apparent, the right ascension decreasing, the declination increasing, with advancing 
spectral type. In noting these apparently systematic trends it is desirable to 
emphasize that there occur also large and apparently irregular variations in some 


Tas_e II 
Variation of drift constants with spectral type 

B8-A3 As-F5 F8-Gs5 Ko-M p.e. 
A, g1°4 gi°'1 92°'9 96°-o +1" 
D, —16°-6 —8°-9 —10°°3 —10°8 : 
AW, 1°628 1°770 17564 1°385 +0°03 
A; 309°°4 301°°o 273°°5 265°°8 + 6° 
Dez ~—gg" 8 —63°°4 e's —74°"4 sd 
hW, 0695 0°735 0889 0°770 +0-03 


* The probable errors quoted in Tables I to IV are approximate values roughly applicable to 
all the quantities in the same line. The computed probable errors are given in the main tables. 
of this and of the previous paper. 
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of the quantities determined—variations which, judged by the probable errors 
involved, are statistically significant, and yet indicate no systematic trend. ‘The 
variation of hW, with spectral type (Table II) is one such case, the value of D, 
for the A-type group another. (It should be mentioned that the probable errors 
are derived from the residuals obtained from the separate regions, and measure 
the degree of accordance of the regional results with the results of the final 
solutions). ‘Two possible reasons for these curious discrepancies may be 
tentatively suggested: firstly, that the two-drift hypothesis is merely a first 
approximation to the representation of proper motion data, and that the data 
are now sufficient tc indicate the short-comings of the hypothesis; secondly, 
it is possible that the traditional grouping into spectral classes and zones of galactic 
latitude is not the grouping which best exhibits the variations in the constants 
involved; for example, a division of the data into giants and dwarfs rather than 
into spectral types, or a distance grouping of the material, might produce more 
homogeneous groups and more regular and systematic variations in the drift- 
constants derived therefrom. 

4. The theoretical values of hV and 6, calculated for each region from the 
constants of the final solutions, are given in Tables VIII a to VIII c under the 
headings AV,, and @,, with the residuals hV.—AV and 6-0. The residuals are, on 
the whole, rather higher than those obtained for the As—M stars. The values of 
A—the distance of the region from the apex—and of S—the distance of the 
calculated apex from the great circle defined in each region by the observed 
direction of the apex—are also shown. The values of f, are satisfactory, only 
four in each table exceeding 10°; those of 8, are rather high, usually exceeding 
10° and attaining nearly 40° in Table VIII a. A noticeable improvement occurs 
in Tables VIII B and VIII c, where the observed regional values are compared 
with the zonal constants. 

5. The linear components €, 7, ¢ of relative motion of the drifts, with the 
resultant 2, and the corresponding right ascension «, and declination 5, of the 
vertex of star-streaming are given in Table IX, together with the galactic co- 
ordinates /, and b, of the vertex nearest to the galactic centre. ‘The outstanding 
feature of the results is the high value of the longitude of the vertex, which exceeds 
that of the galactic centre by 31°. This value agrees substantially with the 
value 351°-5 obtained by Delhaye (4) from his analysis of A-type stars on the 
ellipsoidal hypothesis. ‘Table III shows the variation of the longitude of the 
vertex and of the relative speed © of the drifts with spectral type, the results for 
the later types being obtained from the previous paper (1). ‘The decrease in the 


Tas_e III 
Variation of l, and Q with spectral type 
B8-A3 As-F5 F8-—Gs5 Ko-M 


356°°5 348° 6 341°°7 342°°4 
1-887 2°075 I°QI2 1638 


longitude with advancing spectral type is clearly shown for the first three groups. 
The value of Q follows the same trend as hW, (‘Table II), and the same remarks 
are apposite. 
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The variation of 2 with galactic latitude is shown in Table IV for the same 
groups of stars, in limited spectral ranges, as are usedin Table I. It is somewhat 
perplexing to note that while the later types show a well-marked decrease with 
increasing galactic latitude (a result which was suggested by the results of the 
previous paper), the A-type stars reverse this tendency, and show an equally well- 
marked increase. It may be noted in Table IV that in the non-galactic zones 
the values of 2 suggest that this quantity decreases with advancing spectral type. 
This is not true of the zones of low galactic latitude, nor is it true of the values for 
the whole sky (Table III). The anomalies in Tables III and IV may therefore be 
summarized in the statement that, for the B8—A3 stars, the value of Q is abnormally 
low in regions of low galactic latitude. 


TABLE IV 
Variation of Q with galactic latitude 
B8—-A3 As-F5 F8-G5 Ko—-M p.e. 
Zones I, II 1775 2°209 2033 1°740 +0-06 
Zones III, IV 1987 1-981 1°838 1°578 +0°06 
TABLE V 
Solar motion with respect to spectral type 
B8-A3 As-F5 F8-G5 Ko-M 
Ag 264°°8 266°:o 272°'8 277° '8 
Dy 31°°5 +28°-1 37.4 +40°°3 
hU, 0973 0948 0883 0°770 
n,/ng 1°54 1°22 1°18 1°02 
(1°17) (1°04) (1°08) 


Table IX contains the results for the right ascension A, and declination D, 
of the solar apex, with the corresponding solar speed hU, and the ratio of the 
numbers of stars in the two drifts. The ratio m,/n, is remarkably high. This 
is the only determination of the ratio for the A-type stars alone. Eddington’s 
investigation (5) of the PGC stars (including the A-types) gave the value 1-48. 
Smart and Tannahill (6, 7) from the Cape photographic proper motions found the 
value 1-41 (including the A-stars) and the value 1-10 (excluding the A-stars). 
There is therefore some previous evidence for this high value of the ratio from the 
A-type stars. The variation of the elements of the solar motion with spectral 
type is shown in Table V. ‘Two of the elements show a steady trend. The 
right ascension of the apex increases with advancing spectral type, and the solar 
speed decreases. The ratio m,/n, is high for the early types, and reaches a mini- 
mum in the K-type. This feature is also shown by the Cape proper motions 
(6), the relevant values being shown in brackets in Table V. 

University Observatory, 

Glasgow, W.2: 
1954 October 28. 
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Analysis of the drift curves 
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TasBLe VIILA 
B8—A3 stars—whole sky 
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9 9-9 = Bs As 


+23°0 14°9 
152.99 +12°9 12°6 
30°4 +90°4 29°2 
38°6 +1386 8-2 
140°2 +50°2 9°3 
217799 —I2°I II°7 
229°1 —I10°9 g°2 
231°6 —18-4 14°0 49°6 
219°4 —30°6 30°6 93°1 
221°9 ++ I'°9 1°9 gi‘! 
218°3 — 1°7 1°7 88-9 


Region 
hV,-hV, 
10°7 §4°0 1°23 
2°1 716 1°44 
2'0 94°2 «=I°51 
91 79°9 1°49 
3°7 665 1°39 
77° §=6937°8 0°93 
3°7 §5°O 1°24 
2'0 65°4 1°38 
3°0 80°8 1°50 
6°1 56°3 1°26 
3°6 34°6 0°86 


138°8 
102°6 
150°8 
167°5 
167°8 

76°2 

57°6 


—o'17 163°0 
—0-'06 
+o°o1 
—o'ol 
+0°29 
+0°13 
+0°14 
—0'*42 

0°00 
—O'l4 
—o'l4 


NN 
Co oO 


w 
° 


136°1 
119°9 
148°6 
67°4 
44°0 
67°5 


+o'1o 
—0°23 
+0°22 
+0o°19 
—0o'24 
+0°09 


TABLE IX 


The vertex of star-streaming and the solar motion 
B8-A3 B8-A3, Zones I, II B8-A3, Zones III, IV 


+10°2 71 
— 08 o'7 
—99°O0 31°0 
226'2 —23°8 21°8 
242°0 — 3°0 2°‘! 
2216 + 66 6"1 


2°38 80:3 1°50 
2°9 75‘l 1°47 
6°1 Q2°*I 1°52 
6-2 581 1°29 
67 74°3. «1°46 
14 78:6 1°49 


= 


140°2 
149°2 
III‘O 


N 


tw WWW WwW Ww 
sue W& 


—0'285+0°045 
+1°862 0°044 
+O°IIO 0°047 


° 


aja" 1°*s 
—3°4 14 
1°887 0:047 
356°'5 1""4 
—o'-s  1°4 
264°°8 
+31°°5 
0°973 
1°54 


—0'220+ 0°049 
+1°762 0°057 
—0'004 0°057 
277°'1 2°°0 
+o°1 1°8 
1°775 0°O57 
358°-5 1°°8 
+2°6 1°°7 
266° -3 
+31°°5 
1°074 
1°40 


—0'478+0°079 
+1'908 0-061 
+0284 0-071 
284°-1 2°°3 
—8°-2 2°1 
1°987 0-062 
354°°6  2°1 
“74.22 
268° -6 
+26°-6 
0948 
2°04 


40% 








ORBITAL ELEMENTS OF FOUR VISUAL BINARIES 


W. D. Heintz 
(Communicated by R. v. d. R. Woolley) 
(Received 1954 September 17) 


Summary 


Provisional elements have been derived for = 3050 and OX 358, and the 
orbits of & 2 and = 228 have been recomputed in accordance with the new 
methods of W. Rabe. 





1. & 3050 (ADS 17149; 23"56™-9, +33°27'). New orbit 

This “‘ beautiful pair’? (37 Andromedae) has been measured very frequently 
since W. Struve’s times. ‘The distance decreased from 4” to 1”-5, but the 
observed arc includes neither periastron nor apastron, so that the computation 
meets with some difficulties, even in the application of Rabe’s method for short 
arcs (1). 

The position angles P are expanded in the following: power series: 

P(r) =219°'14 + 6°-587 + 0°-7307? + 0°-06297* + 0°-002174, 
where 1910-0 is taken as zero epoch (osculation) and 10 years as time unit, 
i.e. 7 =(t—1910)/10. Generally, Rabe does not use second and higher differential 
quotients of the coordinates, but if the value dp/dt is determined from the 
position-angle series, as in this case, it is advantageous to include the fourth 
power in the series, so that the second power coefficient can be fairly accurately 
derived. 

The curvature of the apparent arc as given by the most recent normal places 
indicates high eccentricity, but this seems to be inconsistent with the older 
measurements. ‘The computation, using Rabe’s “solution B”’, was repeated 
several times under various assumptions about e and periastron time 7. ‘The 
final solution depends greatly on the accuracy of the 1831 normal place; the 
1792°8 measurement 180° + of Herschel gives O-C= —2°-5. The formulae of 
Russell and Moore (2) yield a dynamical parallax of 0”-065 and a total mass of 
the system of 1-8 solar masses. Table I gives the representation of the 23 normal 
places and an ephemeris. The change of position angle due to precession is 
negligibly small. 

As the star is approaching periastron quickly now, it will be obvious in a 
few years whether the elements derived from the small arc approximate to the 
true orbit. But these elements will at least give a sufficient description of the 
motion of this interesting couple up to the present. 


2. &2(ADS 102; 08 06"-3, +79° 26’) 

Recent observations give deviations up to + 15° from van Arnam’s orbit (3). 
As the components were too close for the observers in the years 1865-85, the 
observed arc consists of two separate parts, and the motion in periastron (1888) 
can hardly be determined with sufficient accuracy. ‘The position angle ¢% of 
the tangential velocity at periastron time, which is important for Rabe’s 
improvement formulae, and which usually may be taken from the previous 
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orbit, had to be changed by nearly 7° in this case. Additionally to van Arnam’s 
17 smoothed normal places, there are 8 new ones, covering the time since 1930 
(Table II). They were not smoothed out, for in such rapid angular motion the 
interpolation process may introduce some error. The resulting elements show 
little difference from those of van Arnam. Although the companion has moved 
more than 300° in position angle since 1830, the direction of the minor axis of 
the apparent ellipse is not yet determined, so that, owing to the high inclination, 
the values of eccentricity and period are still uncertain. From Russell’s 
formulae, the orbit gives a parallax 0”-o0g and a total mass 5-50. 

Before finishing the computations, the writer was informed of a recent orbit 
published by van Biesbroeck, where the period is increased to 419 years. It 
will be seen in a few years whether the recent observations support this 
considerable change. 

3. 2228 (ADS 1709; 2" 10-8, + 47° 15’) 

Previous orbits of this pair showed that the older measurements could be 
hardly represented by the areal constant derived from the more recent observa- 
tions. Rabe (1914) gave full credit to the position angles, while Kuiper’s orbit 
(1930) satisfied the distances. The recent data now allow a solution which 
represents both coordinates fairly well. It should be a useful approximation 
to definitive elements as almost the full orbit is observed. Kuiper’s orbit (4) 
was denoted “‘assez bon”’ in Baize’s new catalogue (§), but gives residuals + 6° 
recently. Using these elements as initial values, and adding 8 normal places 
(Table III), an improvement was easily obtained by means of Rabe’s formulae 
of ‘‘ variation of apsides’’. Period and eccentricity seem to be settled now, but 
an orbit of definitive accuracy will not be available before the companion has 
passed apastron, for the observations at the time after its last apastron passage 
are scanty and discordant. ‘The dynamical parallax comes out to be 0-023, 
giving a total mass of 3:00. Other parallax determinations, referred to by 
Baize (5), yielded values 0”-026 to 0-047. 

4. OX 358 (ADS 11483; 18"33™-6, +16°56’). New orbit 

‘The numerous measurements, from which 14 normal places were formed 
(‘lable LV), cover 45° in position angle only, but with a distinct change in distance. 
The observed are probably contains apastron; so 193c:0 was chosen as 
osculation and apastron epoch, and the position-angle series is (time unit =25 
years): 

P(7) = 182°-95 —9°:137 —0°-4037? + 0°-0567* + 0° 0107". 
As in the previous case, the observations seem to disagree with the law of areas 
somewhat ; the distances were mainly represented. The companion will continue 
its slow motion towards the main star for a long time. The mass of the system, 
3°3 solar masses, results from a dynamical parallax of o”:o21. 


The writer wishes to thank Professor W. Rabe, Munich, for suggesting that 
computation of orbits of these binaries was worth undertaking, and also for 
communicating his unpublished measurements up to 1954:0. 


Commonwealth Observatory, 
Mt. Stromlo, 
Canberra, A.C.T.: 
1954 September 7. 








Normal 
t 


1831 °82 
41°66 
54°44 
66°58 
74°11 
82°83 
9096 
96°94 

1904°17 
09°85 


13°90 
17°65 
21°19 
23°66 
26°22 
30°19 
35°05 
38°20 
40°66 
42°48 
45°98 


49°30 
52°24 
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TABLE [ 
x 3050 
P 3209 e 0°62 i+ 73°:06 w 65°*12 
T 2044°0 a3" °686 NR 353° °68 pe 1125 
places : 
(n) a 1 O-C Po Po O-C 
(11) 190°'2 19080 —o-60 4-02 3-962 4 0:06 
(5) 192°9 19306 -0°16 3°90 3°834 07 
(15) 196°3 196-28 +002 3°59 3°620 — 03 
(16) 200°0 199°76 +oO'24 3°22 3°371 — 15 
(21) 202'1 202°19 —0'og 3°04 3°206 — 17 
(21) 2058 205 °35 +0°45 2°07 2998 . > Of 
(29) 208°7 208°75 —0°05 2°78 2°794 — ol 
(32) 212°0 211-60 +0'40 2°67 2°641 + 03 
(26) 215°2 21548 —o:'28 2°49 2°454 + o4 
(37) 219°3 219°00 +0°30 2°32 2°308 + of: 
(51) 221°9 221°78 +o-12 2°19 2206 — o2 
(33) 224°7 224°58 +0°12 2°13 2‘111 + o2 
(35) 227°6 227°50 +o-'10 2:00 2'027 — 03 
(33) 229°2 229°65 —0°45 1°93 1'971 > a 
(43) 231°6 232°03 —0°43 1°99 I‘gIo + 08 
(48) 235°9 236°03 —0O'l3 1°92 1°824 + 10 
(43) 241°3 241°57 —0'27 1°84 1°726 II 
(52) 244'8 245°24 —0'44 1°74 1°673 07 
(40) 248°4 248-37 +0'03 1°64 1633 o! 
(41) 251°2 250°81 0°39 1°59 1-607 — o2 
(33) 256°0 255 °69 +031 _— 1561 Ol 
(33) 260°3 260°56 —0'26 1°45 1°528 — 08 
(37) 265°0 265 °05 —0'05 1°45 1509 — 06 
Ephemeris : 
a p ?P p 
1954°0 267°8 1°50 1968-0 289°7 1°52 
56°0 270'9 1°49 70'O 292°8 1°54 
58-0 274°! 1°48 72°0 295°7 1°56 
60-0 277°2 1°48 74°0 298-6 1°58 
62-0 280°4 1°49 76-0 301°4 160 
64°0 283 °6 1°49 78-0 304°0 1°63 
66-0 286-7 1°50 80-0 3066 1°66 
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TABLE II 


x2 
P 3009 € 0°575 i+113°°2 Ww 330°°1 
T 1888-2 ao”-686 2 169°-9 fA 1200 
Normal places : 

t (n) P, P, O-C Po 
1829°00 343°4 34391 —o'51 °°75 
1834°00 341°9 342°20 —0°30 
1839-00 340°2 340°22 —0o'02 
1844°00 337°6 337°87 —0'27 
1849°00 335°4 33494 +0°46 
1854°00 332°4 331°13 +1°27 
1859-00 325°8 32580 ° 
1864°00 316°4 317°93 —1°53 


1889-00 181°5 180°57 +0°93 
189400 168°5 169°68 —1'18 
1899-00 160°5 160°49 +0o-°o1 
190400 151°3 151°30 ° 

1909°00 140°8 141°13 —0°33 
191400 130°3 129°29 +o-7I 
1919°00 116-0 115°50 +0°'50 
1924°00 99°3 100°31 —1‘o1 
1929'00 81-0 85-02 —4'02 


1931°47 77°7 77°95 —0'25 
1935°12 76°0 68-41 +7°59 
1937°46 64°6 62°94 + 1°66 
1939°34 54°3 58-92 —462 
1941°38 59°0 54°94 +4:06 
1944°81 49°03 +1°37 
1948-44 444 43°70 +0°70 
1952°22 38-9 39°00 —o'1o 


Ephemeris : 
P 


° 
1954°0 37°0 
56-0 35°0 
58-0 33°1 
60:0 31°4 
62:0 29°8 
64°0 28°3 
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TasBLe III 
= 228 
P 1449-7 e 0°26 i+ 63°°4 
T 1898-2 ao”-908 2 97°°6 
(Equinox 1900) 
Normal places : 
t (n) P,.(1900) P, O-C 
1831 °46 (7) 262°4 262-81 —o'4! 
43°55 (4) 273°9 272°37 + 1°53 
56°54 (7) 284°1 282°73 + 1°37 
64°29 (10) 288-9 29061 —1°7I 
71-02 (8) 300°5 300°51 —o°ol 
76-12 (16) 312°6 312°57 +0°03 
81-51 (9) 336°0 335°36 +064 
84°52 (3) 356°2 356°13 +0°07 
88-78 (14) 251 31°03 —5°93 
89-97 (10) 40°7 39°85 +085 
92°25 (7) §2°3 54°04 —1°74 
95°18 (11) 66-0 67°57 —1°'57 
97°31 (20) 75° 74°98 +002 
99°81 (14) 84°3 81-96 +234 
IQOI ‘13 (13) 87-6 85°16 +2°44 
03°25 (10) go'l 89°74 +0°36 
06°54 (19) 94°2 96:06 — 1°86 
100! (22) 104°4 102°20 +2:°20 
13°14 (17) 107°8 107 ‘66 +o°14 
14°89 (15) III‘! 110°79 +0o°31 
17°34 (22) 114°4 115°37 —0'97 
20°72 (20) 120°4 122-10 —1°70 
22°62 (9) 123°9 126-69 —2°79 
25°23 (12) 131°3 133°38 —2:08 
29°90 (6) I50°1 147°89 +2:21 
1930°24 (21) 149°6 149°06 +0°54 
34°15 (23) 163°7 16401 —0'31 
37°47 (24) 177°2 178-48 —1'28 
39°41 (28) 187°5 186-96 +0°56 
42°20 (29) 200°0 198 -66 + 1°34 
46°64 (21) 214°7 214°82 —o'l2 
49°61 (19) 224°9 22368 +1°22 
52°51 (26) 231°1 230°92 +018 
Ephemeris (equinox of date) : 
P p 
1954°0 234°4 0°65 1968-0 
56-0 2384 0°69 70°0 
58:0 241°9 0°73 720 
60-0 245°1 0°77 74°0 
62-0 248-0 0°80 76:0 
64:0 250-6 084 78:0 
66-0 253°0 0°87 80-0 


1°14 
1°00 
0°93 
0°72 
0°59 
0°38 
0°26 
0°40 
0°39 


0°34 
0°43 
0‘60 


066 
0°63 
0-60 
0°64 
0-72 


0-69 


062 
0°61 
o'61 
0°64 
0°46 


0°56 
0°53 
0°45 
0°44 
0°47 
o's5I 
0°60 





255°3 
257°4 
259°4 
261-2 
263°0 
264°7 
266-3 





Pe 


1°02 
1°08 
1°00 
0°87 
0°69 
0°54 
0°39 
0°34 
0°34 
0°36 


0°40 
0°48 
0°53 
0°59 
0°62 
0°66 
0°70 
o'71 
o-7I 
0°70 


0°68 
0°64 
0°61 
0°57 
0-50 


0°50 
0°46 
0°45 
0°45 
0°47 
0°52 
0°57 
0°62 
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t+ 


++ 


ee 


p 


o'9I 
0°94 
0°97 
0°99 
IOI 
1°03 
1°05 
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TaBLe IV 


OX 358 
P 292) e 0°48 i+ 134°°23 
T 1784:0 a1”-358 R 16°-63 
(Equinox 1900) 


Normal places : 
t (n) P.,(1900) P. O-C 
1845 °30 (6) 220°4 219-78 +0-62 
65°27 (14) 207°5 208 02 —0'52 
74°05 (22) 202°3 203 °93 —163 
81°47 (45) 200°9 200°78 +012 
86-92 (53) 198 ‘6 198°58 +0-02 
96°73 (50) 194°8 194°84 —0°04 


1905°78 (73) 190°9 1QI'52 —o-62 
15°01 (52) 188-1 188-20 —o'1o 
20°45 (52) 185°9 186-26 —o-36 
25°07 (77) 184°8 184°59 | +0o'21 
30°28 (58) 183°2 182-70 +050 
38-04 (82) 180-0 17983 +0'17 
44°36 (53) 178°0 177°41 +0°59 
51°03 (54) 174'8 174°76 +0°04 


Ephemeris (equinox of date) : 
P 


173°2 
172°4 
171°5 
170°6 
169°7 
168-8 
167°9 











MEAN DAILY AREAS AND HELIOGRAPHIC LATITUDES 
OF SUNSPOTS IN THE YEARS 1951 AND 1952 


Royal Greenwich Observatory 
(Communicated by the Astronomer Royal) 


(Received 1954 July 27) 


The following results are in continuation of those given in M.N., 114, 134, 
1954, and are derived from the measurement at Herstmonceux of photographs 
taken at the Royal Observatories of Greenwich and the Cape and at the 
Kodaikanal Observatory, India. 

Table I gives the mean daily areas of umbrae, whole spots and faculae for 
each synodic rotation of the Sun included in the years 1951 and 1952; the 
means for the two years as a whole are included in Table II, which summarizes 
the yearly values since the last sunspot minimum of 1944. 











TaBLe | 
Mean daily areas 
Projected* Corrected for foreshorteningt 
No. Rotation Days A~— oe “ S 
of commenced photo- Whole Whole 
aitiiilisi UT. graphed Umbrae ape Faculae Umbrae epots Faculae 
1950-52 

wer Den o74 86 80 538 749 51 341 gIl 
1302 Jan. 6:07 28 225 1344 842 154 gil 1024 
1303 Feb. 2-41 27 117 533 944 89 556 1173 
1304 Mar. 1°75 27 233 1170 935 172 888 1179 
1305 Mar. 29:06 28 337 2030 1081 237 1449 1314 
1306 §=6Apr. 25°33 27 650 3882 1176 479 2885 1457 
1307. May 22°56 27 607 3678 1387 450 2724 1687 
1308 June 18-76 27 292 1853 1449 213 1385 1712 
1309 6 July 15:96 27 181 1038 1497 130 742 1797 
1310 ©. Aug. 1218 28 118 669 1222 go 520 1445 
1311 Sept. 8-43 27 236 1386 1282 165 953 1508 
1312 Oct. 5-70 27 179 1073 1256 132 807 1455 
1313. Nov. 1°99 27 148 885 1035 III 664 1196 
1314 Nov. 29°30 8 28 121 845 882 86 614 1084 
1315 Dec. 26°63 27 148 840 949 99 566 1128 
1316 Jan. 22:96 28 76 475 749 51 322 880 
1317. Feb. 19:30 27 32 165 519 25 132 633 
1318 Mar. 17°63 27 67 402 540 46 280 670 
1319 Apr. 13°92 27 58 338 660 55 329 794 
1320 ©May 11°16 28 48 254 599 36 194 704 
1321 June 7°37 27 136 727 652 95 522 773 
1322 July 457 27 136 861 757 100 639 932 
1323. July 31°78 927 126 7°9 641 94 543 772 
1324 Aug. 28-o1 27 116 664 602 go 510 735 
1325 Sept.24:28 27 7° 390 476 49 280 559 
1326 §=©Oct. 21°56 8328 77 456 439 56 318 520 
1327. Nov. 17°87. 27 138 799 352 100 589 421 


* Expressed in millionths of the Sun’s disk. 
+ Expressed in millionths of the Sun’s hemisphere. 
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Table III gives for each rotation in the years 1951 and 1952 the mean daily 
area of the whole spots (corrected for foreshortening) and the mean heliographic 
latitude of the spotted areas for both northern and southern hemispheres. 

The mean heliographic latitude of the entire spotted area and the mean 
distance from the equator of all spots are also tabulated. ‘The mean values for 
the years 1951 and 1952 are included in Table IV. 


TABLE II 
Mean daily areas 


“~ 


No. of days Projected* Corrected for foreshorteningt 





dt A A 








_ o ‘ 


Photo- Without Umbrae Whole Faculae Umbrae Whole Faculae 
graphed _ spots spots spots 
366 157 30 160 284 23 126 344 
365 14 102 560 774 78 429 940 
365 ° 389 2381 1794 291 1817 2188 
365 558 3559 2326 405 2637 2894 
366 419 2618 1849 314 1977 2331 
365 485 2873 2140 356 2129 2597 
305 274 1638 1421 203 1222 1754 
365 259 1552 1145 188 1136 1379 
366 24 96 548 594 7° 404 7il 
* Expressed in millionths of the Sun’s disk. 
+ Expressed in millionths of the Sun’s hemisphere. 


Tables II and IV are in continuation of similar tables in Monthly Notices : 
for the years 1874 to 1888, 49, 381, 1889: 1889 to 1902, 63, 465, 1903: IgOI 
to 1914, 76, 402, 1916: 1913 to 1924, 85, 1007, 1925: 1923 to 1933, 94, 870, 
1934: and 1933 to 1945, I1@, 501, 1950. 

The rotations in Tables I and III are numbered in continuation of 
Carrington’s series. 

The chief features of the record for the years 1951 and 19§2 are as follows :— 

(1) Both years are of the descending phase of the sunspot cycle, but in 
1951 the decreasing sunspot activity is largely obscured by the superposition 
of a considerable peak of activity included in Rotations 1306 and 1307 (see 
Table I). In the former rotation, the mean daily area of sunspots was 2885 
millionths of the Sun’s hemisphere and in the latter 2724 millionths. A giant 
sunspot (a late-comer after the four great spots of 1946-47) with its central 
meridian passage on May 16-0 contributed 1757 units to the mean daily area 
for Rotation 1306. 

Another notable spot, the predecessor of that in May, with its C.M.P. on 
April 18-7, had a mean area during its disk passage of 2064 millionths. A third 
spot crossing the central meridian on June 18-4 was of nearly the same area, 
viz. 2100 millionths. 

The mean daily area of faculae for each rotation in 1951 rose fairly steadily 
to maximum values in Rotations 1307-09. 

In the following year (1952) there were no large excursions in the sunspot or 
facula values which for the year as a whole gave greatly reduced values as 
compared with the year 1951, viz. 644 per cent for spots and 484 per cent for 
faculae (see Table II). 
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Tas_e III 
Spots north Spots south 
of the equator of the equator Mean 
— ~ . ¢ A . Mean distance 
Mean Mean latitude from 
. , Mean > Mean : : 
No. Rotation dailv helio- dail helio- of entire equator 
of commenced ns graphic toe graphic _—_ spotted of all 
rotation Ur. —_ latitude area “latitude area spots 
1950-52 ° ° ° ‘ 
1301 Dec. 9°74 268 14°94 73 12°10 9°13 14°33 
1302 Jan. 6:07 671 9°17 240 10°56 +3°98 9°54 
1303 Feb. 2-41 482 10°31 74 11°33 +7°43 10°44 
1304 Mar. 1°75 753 10°65 134 9°07 +766 10°41 
1305 Mar. 29:06 1183 12°41 266 9°79 +834 11°93 
1306 Apr. 25°33 2350 13°14 535 12°90 +8°31 1310 
1307 May 22°56 1660 13°37 1064 13°37 +2:°92 13°37 
1308 June 18-76 302 11°68 1084 10°97 —6:04 I1'I3 
1309 July 15-96 239 7°82 503 8-32 —3°11 8-16 
1310 Aug. 12°18 236 9°52 284 10°96 —1°65 10°31 
1311 Sept. 8-43 729 10°64 223 6°25 |+-6°68 9°62 
1312 Oct. 5°70 234 10°52 574 9°03 —3°37 9°46 
1313 Nov. 1°99 282 13°62 382 9°75 +O°17 11°39 
1314 Nov. 29°30 549 8°59 65 9°79 + 6-64 8-72 
1315 Dec. 26°63 543 5°66 23 12°61 +4°90 5°95 
1316 Jan. 22-96 35 14°36 287 11°13 —8-38 11°48 
1317 Feb. 19-30 20 15°45 112 11°56 —7°38 12°16 
1318 Mar. 17°63 164 5°37 116 5°66 +o-81 5°49 
1319 Apr. 13°92 75 10°55 254 7°53 —3°42 8-22 
1320 May 11°16 10 10°48 184 6:28 — 5°43 6-49 
1321 June 7:37 60 15°67 462 11°00 —7°93 11°54 
1322 July 4°57 120 10°40 519 5°07 —2:17 6:07 
1323 July 31-78 335 10°50 207 9°47 + 2°87 101! 
1324 Aug. 28-o1 209 5°81 301 6-29 —1°34 6-10 
1325 Sept. 24:28 221 9°57 59 6°34 +6-22 8-89 
1326 Oct. 21°56 278 5°16 40 13°38 2°80 6-21 
1327 Nov. 17°87 528 6-09 61 10°30 +440 6°53 
‘Tas_e IV 
Spots north Spots south 
of the equator of the equator Mean 
- a ,_  — “a . Mean distance 
Mean Mean latitude from 
: Days Mean helio- Mean helio- of entire equator 
Year : daily ; daily 
y aily 
photographed et graphic ie graphic spotted of all 
latitude latitude area spots 
° °o ° ° 
1944 366 42 19°00 83 22°81 —8-70 21°53 
Old Cycle 3 4°17 6 7°61 —I1'l4 6°72 
1944 3664 28, Cycle 35 22:18 77 24°02 —9g ‘63 23°45 
1945 365 121 20°13 309 «=. 20:26 —8-92 20°22 
1946 365 1127 20°74 690 ~=—:11879 T§°73 20°00 
1947 365 992 16°58 1645 17°86 —491 17°38 
1948 366 936 14°92 1041 13°53 —o:06 14°19 
1949 365 1178 14°30 gst «12°14 +249 13°33 
1950 3836. 365 771 13°62 452 13°05 3°77 13°41 
1951 = 365 730 11-62 406 = 10°78 +362 11°32 


1952 366 195 7-60 208 8-37 —o-64 8-00 
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The ratio of mean corrected areas of faculae/sunspots for 1951 is 1:21 and 
for 1952, 1:76. (See comments, Monthly Notices, 112, 572, 1952.) 

For other characteristics of the years 1951 and 1952 see Monthly Notices, 
112, 332, 1952, and 113, 382, 1953. 

(2) There were two days without spots in 1951 and 24 days in 1952. 

(3) The mean weighted latitude of all spots in 1951 (‘Table IV) crept 2°-1 
nearer the equator as compared with 1950. A further decrease of 3°-3 followed 
in 1952. ‘The latitudes of the three giant spots of 1951 in order of date were 
+12°, +13° and —12°. 

The northern hemisphere in 1951 had much the greater mean spotted area 
but this was largely due to the great spot in May which itself contributed 
130 units to the mean for the year. 

(4) The number and distribution, northern and southern hemispheres, of 
spot groups of (a) two days’ duration or longer; (b) one day only, are as 
follows : 

1951 1952 
(a) (b) (a) (6) 
Northern spots 142 39 56 26 
Southern spots 115 29 73 20 
Total 257 68 129 46 


The range in latitude for the longer-lived spots in 1951 was from 1° to 25°. 
In 1952 the range was from 1° to 23°. The latitudes of the 1-day spots ranged 
from 2° to 26° in 1951 and from 2° to 22° in 1952. 


An appended table gives the mean daily areas of sunspots and of faculae 
{corrected for foreshortening and expressed in millionths of the Sun’s hemi- 
sphere) for each calendar month of the years 1951 and 1952. 


1951 1952 
Month Spots Faculae Spots Faculae 
Jan. 746 1005 464 1034 
Feb. 627 1173 305 774 
Mar. 845 1210 142 639 
Apr. 1600 1347 354 782 
May 2899 1422 257 703 
June 2477 1628 393 800 
July 780 1942 654 906 
Aug. 664 1534 689 768 
Sept. 960 1548 307 
Oct. 757 1427 339 518 
Nov. 666 1205 455 
Dec. 602 1091 477 


Royal Greenwich Observatory, 
Herstmonceux Castle, 
Sussex : 


1954 July 15. 








STELLAR PARALLAXES DETERMINED PHOTOGRAPHICALLY 
AT THE CAPE OBSERVATORY (NINETEENTH LIST) 


Royal Observatory, Cape of Good Hope 


(Communicated by H.M. Astronomer) 
(Received 1954 October 4) 


This list is in the same form as the previous lists beginning in M.N., 89, 
497, 1929. 

The units in the last five columns are 0”-001. Two values of the proper 
motion in R.A. are given. The first (as also the proper motion in Dec.) 
is derived from meridian observations or other sources; the second is the 
relative value obtained from the parallax plates. The two values of probable 
error entered in the last column refer respectively to the parallax and to the 
photographic proper motion in R.A. 

Nearly one-half of the stars in the present list are faint stars with a proper 
motion exceeding 0”-5 annually, given by Luyten in Vol. III, No. 1 of the 
Publications of the Astronomical Observatory of the University of Minnesota. 
Approximately one-third are stars selected from the General Catalogue or from 
the Cape Astrographic Zone as having an annual proper motion in excess of 
o”-2. A few stars of considerably smaller proper motion—relatively bright 
stars of early and late spectral types—are also included, but none of these has 
a noteworthy parallax. The remainder of the list includes stars of interest 
derived from various sources, particularly Vyssotsky, Ap. 7., 97, 381 and 
104, 234; while three faint stars, Nos. 1729, 1752 and 1792 appear on account 
of their relative proximity to Nos. 1730, 1753 and 1791 respectively, and the 
consequent original uncertainty as to the identity of the proper-motion star in 
each case. 

Where possible, the magnitudes given are taken from the Henry Draper 
Catalogue; otherwise these are generally as given by the authority quoted for 
the proper motion. 

Nine of the stars have parallaxes exceeding 0”-1. Of these, six are stars 
given by Luyten. The other three are BD + 4°4048 (No. 1773) and BD + 1°4774 
(No. 1799)—two of Vyssotsky’s dwarf M stars—and CoD —32°16135 
(No. 1778), a star with only moderately large proper motion. 


. Mag. Proper Motion P.E. 
mm oe - Vis. Type RA»Dec. Unit Parallax PE. 
— <n Phot. (1900°0) R.A. Dec. Weight 
hm °s 
1701 LPM 6 ve m © 03 53 + @o +130 tar +83 it 8&8 
—27° 16 12'0 —27° 410 + 692 + 4 
1702 GC 452 6:44 Go © 17 44 + 388 +66 +27 +48 +09 
9 Ceti 7°00 —12° 46’-0 + 407 x's 
1703 LPM 26 aes m © 27 36 + 560 -—100 +24 + 21 +7 
—63° 51 II‘! — 63° 38'°3 + 506 = 
1704 LPM 29 9°4 Gs © 30 21 + 920 —s80 +426 +59 i9 
—64° 51 10°! —64° 13'°7 + 868 = 5 
1705 GC 721 8-58 Gs © 31 16 + 372 —142 +23 +74 =+%+7 
— 49° 138 9°2 —wars +3 2 
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No. Name and — . 


D.M. No. Phot. 


Proper Motion P.E. 
Unit Parallax 
Weight 


R.A., Dec, 


Type (1900°0) 


R.A. Dec. 


GC 865 
B Ceti 


2°24 Ko 
3°24 

GC 922 5°62 

—23° 293 6°9 

LPM 53 

—10° 216 

GC 1296 

— 48° 259 

LPM 81 

— 68° 77 


GC 7127 
—15° 1126 


LPM 248 
— 44° 3045 
GC 9264 
—61° 754 


GC 9280 
—67° 686 


LPM 257 
—63° 703 
FS 3446 
—44° 3764 
GC 11306 


—1° 2005 


GC 11398 
—68° 771 
GC 11837 
—6° 2664 


GC 11906 Mb 
— 16° 2541 . 


GC 12307 
—4° 2490 
GC 12554 
—52° 1904 
GC 12996 
—43° 5262 


GC 13043 
—21° 2802 


GC 13168 
—10° 2857 


LPM 322 
— 38° 5760 


hm =°5s 
© 38 34 
— 18° 32’*1 


© 41 14 
—23° 04’'1 


© 56 57 
—10° 24'7 


© 59 56 
—48° 28-6 


1 38 20 
—68° 10-2 


5 37 22 
—15° 40°°6 


6 54 56 
—44° og 2 
6 58 50 
—61° 12’'0 


7 00 O1 
—67° 46'°7 


7 12 15 
—63° 112 


7 46 34 
—44° 209 
8 12 48 
—O1° 29'°3 


8 19 07 
—68° 40’-9 


8 32 57 
—o6 27°°S 


8 35 18 
—16° 56'°7 


8 49 22 
—05° 03'°3 
9 OI 23 
— 52° 066 


Q 21 22 
—43° 328 


9 22 44 
—21° §4°°3 


9 27 36 
—10° 44'°6 


9 34 37 
— 38° 55'*4 


231 + 40 
221 


198 
184 


200 
213 


7O 
72 





H- 


H- 


Ib He 


L 


H- | 


H- 


tL 


HH HEHE HEE 


He He 


H- H- 


H- H- 


He He 
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‘ Mag. Proper Motion P.E. 
No. a = Vis. Type RA»Dece Unit Parallax PLE. 
ie. Be, Phot. (1900) = RAL = Dec. Weight 
hm °3 
1727. LPM 336 site m 10 06 54 + 480 —640 +22 +118 +7 
12°2 —0o2° 10’°0 + 528 ~ % 
1728 LPM 341 —- k 10 II 53 — 400 —600 +20 +106 + 6 
12°4 —11° 26:8 — 388 + 3 
1729 N.p. No. 1730 " od 10 18 44 aa ~~» <6 =— 89 - 9 
12 -—wrs + % 3 
1730p LPM 345 me m 10 18 48 + 370 —440 +17 + 58 + 6 
11°5 —59° 39°00 = + «(357 £3 
1730f LPM 346 abe m 10 18 50 + 370 —440 +29 + 57 +12 
13°5 —59° 389 = + «(357 + 6 
Mean re 
t 357 = 2 
1731 LPM 361 int g5 10 47 33 — 240 —750 +21 +45 +6 
—1° 2457 11°O —a" 35. — 378 t 4 
1732 LPM 366 ran K2 II 02 28 + 350 —360 +26 + 14 +7 
— 30° 8970 11*4 — 30° 16-2 + 331 + 5§ 
1733 LPM 368 dis ss II 04 44 — 790 —470 +23 +106 47 
—23° 9765 11°8 —24° 02’°6 — 769 +4 
1734 LPM 386 ees m tr 22 §2 + 510 —820 +32 +45 +9 
13°0 — 08° 36'°7 + 565 + 5 
1735 GC 15886 5°63 Ma II 30 24 — 9 — 6 +332 + § +41 
— 46° 7199 Ya -— a3. ~- =? 
1736 ZC 9655 8+3 Ma _ sit 30 47 — 54 “st 2 Rise. ae 
— 46° 7203 9°6 —-~ 9 —- +7 
1737 LPM 398 — m II 41 35 + 760 —780 +24 +79 =t+ 7 
12°6 —13° 26-0 + 9727 : 5 
1738 LPM 402 I1'9 M4+ 11 48 10 — 150 —520 £17 + 79 x= ¥ 
13°8 —06° 483 — 157 +2 
1739 Ross 925 _ oi I2 02 09 — 490 —220 +22 + 10 + 6 
— 4° 3208 9°4 —05° 10'°3 — 246 t 2 
1740 GC 16584 2°88 B3p 12 03 10 — 37 — 20 + 23 + a3 : 
6 Centauri 3°05 —50° 099 — 24 + § 
1741 Ross 931 9°91 Ko 12 07 18 — 250 —260 +28 + 36 - 9 
—5° 3450 10°91 —ef 4y°F  — B83 + 4 
1742 GC 16724 3°08 B3 12 og 50 — ae r25 ° 
8 Crucis 2°91 —58° 116 — 13 + § 
1743 LPM 430 aia G5 12 23 11 — 500 — 30 £26 + 34 + 7 
—16° 3469 10°5 —16° 21'°5 — 544 . + 23 
1744 Ross 974 és as 12 56 53 — 10 —570 +17. + 75 1. 6 
12 —o1° 32':8 - 34 - 2 
1745 LPM 475 Sd m 13 20 19 — 510 —100 +23 +116 + 6 


—27° 9225 12°0 —27° §1'"1 — 493 aS 
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Mag. 


ap | ; ak 
pu aOeneeeT s Proper Motion P.E. 
is. 


A R.A., Dec. Boge 
. Type : pcviisaiailactiaiaiiae Unit Parallax 
D.M. No. Phot. (1900°0) RA. Dec Wess 


hm =°s 
LPM 477 ee < 13 22 26 — 480 —120 
—27° 9236 - —37 46%. — 478 


LPM 478 me I 13 23 11 + 160 
Ross 486 ° —oO1° 50°"! + 161 


FS 11167 dive ain 13 29 33 — 358 
CAZ —47° 951 , —46° 11'-2 341 


FS 11346 ii 3 13 39 03 223 
— 47° 8579 , —47° 34/0 - 209 


GC 18666 ° 13 43 39 46 
— 33° 9358 4 —3a Sy *t + a 


LPM 509 aa 14 03 41 - 470 
—30° 11195 ' — ae. ae 


S.p. No. 1753 wal one 14 13 51 
—06" 50°'2 


LPM 519 ii 14 13 5 
Wolf 534 , — 06° 500 


LPM 521 5 14 14 50 
Ci 20 851 , —o9° 088 


ZC 12226 ° 35 14 15 04 
— 42° 9252 , —43° 01°7 


ZC 12336 oe ; 14 22 35 
— 46° 9347 - —46° 17'°8 


LPM 530 eit 14 26 50 
—27° 9894 : —27° 48'°6 


LPM 554 —_ < 15 03 42 
—19° 350 


LPM 557 om 15 07 25 
—21- 35°° 


GC 20537 "95 rs I5 II 32 
—7° 3992 5.p. , —0o7" 54" 


—7° 3992 N.f. , j 15 11 


Mean 
LPM 569 
— 26° 10870 


LPM 595 
—70° 2196 


LPM 601 


LPM 612 
—14° 4454 





Stellar parallaxes determined at 


Name and 
D.M. No. 


GC 23708 
« Arae 


FS 17263 
—40° 11804 


— 1° 3474 


LPM 676 


FS 18061 
—51° 11572 


GC 25488 
— 43° 12699 


Ross 142 
Ci 20 1111 


ZC 18044 
—50 12397 


Wolf 1055 
4° 4048 


FS 18978 
— 46° 13364 


FS 18979 
— 46° 13365 


Mean 


LPM 724 


LPM 736 


GC 28425 


— 37° 13741 


— 32° 16135 


ZC 19181 
- 46° 13835 


GC 29782 
~43° 14464 


FS 19554 
CAZ —51 


FS 19595 
CAZ —48 


LPM 806 


Mag. 
Vis. 
Phot. 


2°97 
2°80 


10°o 


88 


Type 


R.A., Dec. 
(1900°0) 
> = 

17 24 07 

—49° 47'°8 


17 38 33 
— 40° 45°°3 


18 14 39 
—o1° 58':8 


18 18 27 
—0o5° 11°°5 


18 19 45 
—sr oF 3 
18 32 24 
—43° 16°*3 
18 45 O1 
+ 02° 59° 
19 08 22 
—50° 21°°7 
19 12 03 
to5° 018 
19 56 57 
— 46° 04°°4 


19 56 57 
— 46° 040 


20 O1 19 
—19° 41'°6 


20 15 58 
—o6° 44°°3 
20 20 24 
— 37° 43'°6 
20 35 36 
—32° 46'°9 


20 59 54 
— 46° 40'°3 


21 Il 54 
—43° 45°°3 


21 27 59 
—50° 00’'! 
21 35 07 


— 47° 28':9 


21 55 43 
—7O° 24°°3 


the Cape Observatory 


Proper Motion 


R.A. 


32 
28 


231 
206 


26 
8 


Dec. 


77 


P.E. 


Unit Parallax 


Weight 
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ae Mag. Proper Motion P.E. 
No. Name and Vis. Type RA.Dec, —sUnit Parallax 
Phot. (1900°0) R.A. Dec. Weight 


D.M. No. 


hm 8 
1784 LPM 808 ee 21 56 22 + 920 + 50 +13 + 78 
—19° 58’-0 + 924 


1785 LPM 817 oii 22 09 07 80 -—620 +19 
Wolf 1332 “ —og° 13'°4 106 


1786 LPM 831 nem 22 20 36 —600 +21 
—21° 19'°7 


1787 FS 19954 “a 22 25 27 +24 
CAZ —51° 35 , —51° 441 


1788 LPM 837 22 32 56 +31 
~— 359 


1789 LPM 840 ate 22 42 32 
—37° 18'°6 


1790 GC 31913 ‘ 22 47 43 
— 49° 13988 ? —49° 07'°7 


1791 LPM 847 pes 22 55 31 
—55° 020 


1792 5.f. No. 1791 ee on 22 55 32 
—55° 02'°2 


1793 LPM 854 ee 23 04 48 
—26° 16501 . — 26° 28’5 


1794p LPM 859 ° 23 08 52 
—9° 6149 “ —og° 28’*1 


1794f LPM 860 , 23 08 53 
—9g° 6150 ° —og° 28°5 


Mean 


1795 FS 20274 poe ies 23 11 15 
CAZ — 42° 96 ° —41° 197 


1796 —2° 5958 eas 23 21 55 
—o1” 50°1 


1797 LPM 875 sail 23 29 03 
Ci 20 1433 . —oo° 207 


1798p LPM 881 ia 23 37 50 
—24° 17814 . —24° 394 


1798f LPM 882 oat 23 37 53 
—24° 40'8 


Mean 
GC 33053 
+1 4774 


LPM 895 
—6° 6318 
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Vol. 
Notes 
This star was originally selected as a reference star for 8B Tucanae (No. 1401). 


This star has common proper motion with No. 4. 
A difficult double, of which the mass was measured. The observations extend over 18 
years. 


This is CPD —21° 5912. 
This is ADS 9564. 


This star, which is itself a close double, has common proper motion with No. 1189. 
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